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The Absorption of Light by Alkali Metals 
By-Pa N, BUTCHER 


Department of Mathematics, Imperial College, London 


Communicated by H. Fones; MS. received 6th April 1951 


ABSTRACT. The theory of the absorption of light by photo-excitation of the valency 
electrons in cubic metals is developed. The matrix elements vanish when calculated with 
respect to the Wigner—Seitz approximate wave functions for the occupied states in the first 
band, and free electron wave functions for the states in the higher bands. Calculations of 
the contribution to the absorbing powers of the alkali metals due to photo-excitation are 
made, using the approximation of nearly free electrons to obtain all the unperturbed wave 
functions. The calculated values for sodium, potassium and rubidium agree quite well with 
those observed if the first non-vanishing Fourier coefficient of the self-consistent potential 
field has the values 0-323, 0-305 and 0-350 ev. respectively. 


§1. INTRODUCTION 

BSORPTION of visible and ultra-violet light by metals is due to two processes, 
Af scattering of the oscillating electrons by the vibrating lattice, and 
secondly photo-excitation of the valency electrons. ‘The former process 
has been treated by considering the damped oscillations of free electrons (see, 
for instance, Seitz 1940). This theory predicts for the alkali metals absorbing 
powers which are only about one-tenth of the observed values. Hence in these 
metals it appears that the absorption is primarily due to photo-excitation. We 
shall, therefore, make a detailed investigation of this effect. In previous treat- 
ments (Wilson 1935, Fan 1945) it has been assumed that the metal is isotropic with 
regard to absorption. We shall present the theory so as to verify this assumption. 
We shall show that the matrix elements vanish when calculated with respect to the 
Wigner—Seitz approximate wave functions for all the occupied states in the first 
band, and free electron wave functions for the states in the higher bands. The 
contribution to the absorbing powers of the alkali metals due to photo-excitation 
will be calculated, using the approximation of nearly free electrons to obtain all the 

unperturbed wave functions. 


§2. THE UNPERTURBED ELECTRONIC STATES 
The self-consistent field V(r) of a cubic lattice has the periodicity and symmetry 
of the lattice. ‘The unperturbed wave functions are the solutions of the equation 


SP re Wa (1) 
om Fee SF eereasines dent sac 
having the usual time dependence, which are triply periodic with period L, where 
L?=V is the volume of the metal. We shall normalize these wave functions, so 
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that they form an orthonormal set in the cube of side L. It is convenient to use the 
reduced zone scheme to identify the electronic states. ‘The energy levels fall into 
a series of bands which are distinguished by associating with each a positive integer 
l. The states within each band are distinguished by the wave vector k of the 
corresponding wave function. This wave vector may be restricted to take all 
values 27n/L, where n is any triple of integers, which lie within the first Brillouin 
zone. ‘Thus we may speak of the state (J, k) with energy E,(k) and wave function 
Y(k, r, 2) =,(k, r) exp {—iw,(k)t}, where E(k) =fw,(k). The space part of the 
wave function has the Bloch form %,(k,r) =e™*'u,(k,r), where u,(k,r) has the 
periodicity of the lattice. As an immediate consequence of the symmetry of 
V(r), we find that u,(k, r) and £,(k) have the respective symmetries expressed by 
the relations 


u(k, Rx) =u,(Rk,r), ...... (2) and E(k) = 5,RR) te eee (3) 


where R is any one of the operations of the cubic group. 

We shall be concerned ultimately with the alkali metals which have one valency 
electron per atom. At the absolute zero these electrons fill completely the lower 
half of the energy levels in the first band. Hence the number of electrons per 
unit volume of the metal per unit volume of the occupied region of k space is 
1/47°. ‘The modifications necessary in the following theory, when the metal is 
not monovalent, are trivial. When the frequency of the light is not near the 
absorption limit, it is sufficient to consider the metal at the absolute zero. ‘This 
is because the electrons in the tail of the Fermi distribution, which exists at 
ordinary temperatures, absorb only a small fraction of the total light energy 
absorbed. 


$3: THEORY OF THE RATE OF EXCITATION 


We will suppose that an electron in the state (1,k) is perturbed by a wave of 
visible or ultra-violet light, which has scalar potential zero and vector potential 
A=4(ae*”+a*e'). The vector a can be taken as constant throughout the 
metal, because we can always choose the fundamental cube sufficiently small to 
justify this assumption, while still sufficiently large to contain a large number of 
electrons. 

It is well known that the electron can only make transitions to states in the 
higher bands with the same wave vector k, i.e. to states (/,k). We will suppose that 
the perturbation is switched on at time ¢=0, and denote the probability of tran- 
sition to the state (/,k) in time t by P(/,k,t). P(l,k,t) may be calculated by 
applying Dirac’s method of the variation of constants to the perturbed wave 
equation. Using the periodicity of the functions u,(k, r), the result may be ex- 
pressed in the form 


PALS /oca\e 1 —cos [w,(k) — w,(k) — wt 
P(t) = 5 (=) lac K| Vik = a[oxk\ son ee ; 


where (1,k| V|1,k) = | uy*(Kk, 2) Vuj(kyt)dV. eae (4) 
‘ 


We have neglected the terms of P(/, k, t), which involve w,(k) — w,(k) + w, because 
they contribute nothing to the integrals in § 4. 
The symmetry property (2) of the functions u,(k, r) leads to the useful result 


(i, Rk V[1,RkY=(LK[RV|1,K)...  wervns (5) 
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The isotropy of a metal with a cubic lattice with regard to absorption depends on 
this relation, which allows us to simplify considerably the final expression for the 
absorbing power. 

The probable rate of excitation at time ¢ to the state (/,k) is d[P(/,k, t)]/dt. 
This can be written in a very convenient form by putting w,(k) — w,(k) — w = ©,(k). 


Thus we find 

d Pear een\" , | sin Q,(k)t 

aPkt)= 5 (=) la. (k|V|1,k)| | aa i: 
Since the perturbation is due to radiation, it is possible to choose ¢ large enough 
to make the last factor of d[P(J,k, t)]/dt approximately a 6 function, without 
invalidating the variation of constants method by making the transition pro- 
babilities too large. Hence for large t we have approximately 


d 7 2 
ah k j= 5 (=) ever< 2, ke AV) Lk) |20/ OKT eas (6) 


§4. THE INTRINSIC ABSORBING POWER 
Summing (6) over all the occupied states in the first band, and then over all 
1#1, we find that the total rate of excitation of photo-electrons per unit volume of 
the metal is 1 


n= 3(5) y | la- (1,k|V]1,k)P38[Q,(k)] dk... (7) 


We choose dk to be a cylindrical volume element between two adjacent surfaces of 
the family ©,(k)=constant, with base dS, on Q,(k)=constant and length 
ds = dQ,| V,.Q,(k) |, where V,=(0/0k,, 0/0k,, 0/0k,). We can then carry out 
immediately the integration with respect to Q, and obtain 


Awl (oe Ja. (,k|V|1,k)[? 
n= 2 (saz) > ean or] Sama slew (8) 


The summation is over those bands, except the first, for which the surface Q,(k) =0 
passes through the occupied region of k space. ‘The integration is over the part of 
the surface (,(k) =0 within the occupied region, which we denote by S,. Using (3) 
we see that this surface can be divided into 48 regions, each of which can be 
obtained from a single fundamental region, by means of one of the 48 operations 
-of the cubic group. Let J, denote the integral over the whole surface, Sj, the 
fundamental region of S,and RS, the region generated from it by the operation R. 
We have 
Ja. </, a ja. Zk| V1, k)>) ei Ais (Lyle | Va)iL, kee 
= i Sapa (re hype Se) 


=| la. </, Rk| V | 1, Rk) |? 
ae | Vi,02,(k) | ze 
Using (3) and (5) this becomes 


dS. 


1 
=| ap hala: GRIRVIL EDP dSy. 
meg @la arid aciaie 
Writing out the summation explicitly, we find that all the cross products 
between different components of a cancel, leaving 


1,=16|aP | [<2 kV [1,k)>P 


LS AR 9 
Sie [Vx 22,(k) |i 4 ( ) 
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Finally, substituting (9) in (8) we have 


2g LR apg i CE IE 
n=2(——) lal =| aon 2 diet oe (10) 


One quantum of light energy is absorbed for each electron excited. Hence 
the mean rate of absorption of electromagnetic energy per unit volume of the metal 
is 


R, =nhw =2ihw (=) lal, | ee FESO RE (11) 
Since the scalar potential of the light wave is zero, the electric field is E= —Ale 
Hence the mean square electric vector is 
F2=%|.al?/2c2a «a) sean) ahaa (12) 
Eliminating |a|? between (11) and (12) we have 
Sane Ee V|1,k 
R,=— =(=) ES | See Reet (13) 


We now see that R, is proportional to E’, 1.e. the metal is isotropic with regard to 
absorption due to photo-excitation. 

If the frequency of the light is sufficiently small, none of the surfaces Q,(k) =0 
will pass through the occupied region of k space, and there will be no absorption. 
The frequency at which absorption begins, i.e. when Q.(k) =0 touches the surface 
of the occupied region and transitions to the second band commence, is called the 
absorption limit, and denoted by vy. 

In classical electromagnetic theory the intrinsic absorbing power of a homo- 
geneous isotropic medium o(w) is defined by R= =o(w)E?, where R is the mean rate 
of absorption of electromagnetic energy per unit volume. R may be separated 
into a small part R, due to scattering by the vibrating lattice, and the part R, due to 
photo-excitation. We define an intrinsic absorbing power for each of these by 
the relations 

Reo wR a ele (14) and Rese (eo) Et, kh ae (15) 


so that o(w)=o,(w)+o,(w). Comparing (13) and (15) we obtain the final 
expression for the absorbing power due to photo-excitation 


2e\? Kiklv Lip 
aes (=) Al a TN Se (16) 


§5. APPROXIMATIONS TO uw ‘k, r) AND E, (k) 


Henceforth we shall be concerned only with the alkali metals sodium, 
potassium, rubidium and caesium. The Wigner—Seitz treatment would 
provide approximations to u,(k,r) and £,(k) for all the occupied states of these 
metals. In this treatment u,(k, r) is expanded as a power series in k within each 
atomic sphere. Wigner and Seitz obtained the approximation 


u(k,r)=f,([r])+k-rf((rl), |r] <7, 


where 7 is the radius of the atomic sphere. Since u,(k, r) must have the periodi- 
city of the lattice, we find f,’(79) =0, f,(79) =0. Accurate approximations to the 
wave functions of the states in the higher bands are difficult to obtain. It is 
reasonable to suppose that these wave functions are those of free electrons, so that 
u,(k,r) is independent of r for />1. ‘The integral in (4) is 3V/4x7r,° times the 
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integral (of the same periodic integrand) over an atomic sphere. It is easily seen 
that the latter integral vanishes when the above approximations to u,(k, r) are used. 
Thus, in this case, there is no absorption due to photo-excitation. 

These considerations suggest that the relevant properties of the wave functions 
of the occupied states are their nearly free electron characteristics. We shall 
therefore obtain all the unperturbed wave functions by using the approximation of 
nearly free electrons. It is simpler to use the extended zone scheme, instead of the 
reduced one, to distinguish the electronic states in this approximation. ‘The 
various bands of the reduced zone scheme are each allotted to a separate zone of k 
space, by adding appropriate vectors of the form 2rn/a to the wave vectors 
associated with the statesin the band. The electronic states are then completely 
distinguished by the associated wave vector k and the band number / becomes 
unnecessary. ‘Thus we may speak of the state (k) with energy E(k). It is 
convenient to adjust the wave vector which distinguishes the electronic states so 
that E(k) tends to the free electron value (%k)?/2m when V(r) tends to zero. 

The periodicity of V(r) allows us to expand it as a triple Fourier series 
V(r) ==, V, exp {—71(27/a)n.r}, where ais the lattice constant. Asaconsequence 
of the cubic symmetry of V(r) we find that Vp,=V, for all R. It is therefore 
convenient to denote every one of the equal coefficients Vz, by the single symbol 
V,. Moreover, the structure factor for body-centred cubic lattices is 
1+cos7(n,+n,+n,), so that V,, is zero if n,+n,+n, is an odd integer (Mott and 
Jones 1936). Provided k does not approach any zone boundary, the correspond- 
ing normalized wave function has the space part 

o(k, r) =e '"u(k, r) 
pet 2nma* V ; 
= al Sa Gace oat ex? (—im/a)m. x} |, Ney: (17) 
to the first order of V,, (Wilson 1936). 

We still have to consider E(k). It is well known that for sodium this has very 
nearly the free electron value with only small energy gaps. There is little satis- 
factory information about E(k) for potassium, rubidium and caesium. We will 
therefore assume that E(k) is given approximately by the free electron formula, 
with only small energy gaps, for these metals as well as for sodium. ‘Thus, 
neglecting the energy gaps, we have 

(keto kie(hik)e 2m 1") eer Ie Fi, 28 (18) 
As the alkali metals are monovalent and have body-centred cubic lattices, if E(k) is 
given by (18), the occupied region of k space is the Fermi sphere with radius 
Ry (O77)! fa. 

In order to calculate o,(w) we shall neglect the effects of the zone boundaries, 
and use (17) and (18) for all k. 


§6. CALCULATION OF a, (w) 
The selection rule for body-centred cubic lattices in the notation of the extended 
zone scheme is | 
k’ =k+ (27/a)n, een oS CVEM, INtCVEL,. s.5->- «:- (19) 
where k and k’ are the wave vectors of the initial and final states respectively. 
Writing (16) in the notation of thefextended zone scheme we have 
h (2e \? | <k + (27/a)n| V|k)/* 
—)} 2% dS cy cece ee 
WwW (=) ¢ lg [Vx Qu(k) hk dS nfs ey) 
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where 
{k+(2z7/a)n| V[k)= | een) Bo EK (27/2) Bs f) Ve(K, t) eV aes (21) 
JV 
_ and O.(k) =k (27/a)n] = w(k) on 9 Gr eem (22) 


Consider the totality of the surfaces of Q,(k)=0, for which n,+n,+7, is an even 
integer and |n| has the same value x. The summation in (20) is over all values of n 
for which this composite surface passes through the occupied region of k space. 
The integration is over the fundamental region S.,, of the part of the composite 
surface within the occupied region. 

Substituting (18) in (22) we find 


On(k)=(jmaz)(an, K-70) —@,.) (Ae eas (23) 

which gives after differentiating 
VO(k) = (Hid) nw OP ee (24) 
Substituting (17) in (21) we find, to the first order of V,,, 

—imavV,,, 

For the alkali metals and light frequencies in the visible and ultra-violet, the — 
only planes of the set (2,,(k) =0 which pass through the Fermi sphere are those for 
which n=1/2. Hence the sum in (20) reduces to the single term with n= 1/2. 


Let wy and w, be the angular frequencies of the light at which S/2- enters and 
leaves the Fermi sphere respectively. We find immediately 

Wy =(1/2h|ma*)(/20—ak,,) and = w,=(+/2h/ma?)(4/27 + ak,,). 
The absorption limit is vy =w,/27. Inserting the values a=4-24, 5-25, 5-62 
and 6:05 x 10-$cm. for sodium, potassium, rubidium and caesium respectively, 
we find the absorption limits vy =5-11, 3-23, 2-94 and 2-37 x 1014 sec“! respectively. 
It appears that the absorption limits of sodium and potassium given by Fan (1945) 


in his Table 1 should be interchanged. If wy»<w<w,, the fundamental region 
S\/op iS a quadrant of a circle with radius 
1 =[Rp»? — (wmal/2h — 4/21] a)?}1? = (ma/1/2h)[(w — &9)(w, — w) }#?. 
Substituting (24) and (25) in (20) we find that the integrand is constant over S \/2¢ 
so that the integral can be immediately evaluated. The final result may be 
written in the form (Wilson 1936, equation (197)) 
4./2maen* |_, W — W_)(w,—@ 
NCD I rare rere fae 
In the Figure we show oa,(w) calculated from (26) for sodium, potassium, 
rubidium and caesium by assuming the values |V \/.| =0-323, 0-305, 0-350 and 
0-522 ev. respectively. According to the approximation of nearly free electrons 
2|V ./2| is the energy gap across the surface of the first Brillouin zone. This 
interpretation is perhaps not relevant here because | V,/.| arises in (26) from the 
wave functions of states which are not near the planes bounding the first zone. 
The circles are the observed values of the absorbing power less the contribution 
due to scattering by the vibrating lattice. The former was calculated from the 
observed values of the refractive index 7 and the extinction coefficient Ky, obtained 
by Ives and Briggs (1936, 1937), by using the well known expression 


o(w) =nKyw/2z. 
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o,(w) was calculated from the classical formula (Seitz 1940) 


e{n 2 
o,(w) = ant ae pets (28) 

where p is the ordinary resistivity and my) =2/a is the number of electrons per unit 
volume of the metal. p was given the empirical values for the bulk material, 1.e. 
4:35, 6:62, 12-0 and 19:0x 10-*ohmcm. for sodium, potassium, rubidium and 
caesium respectively (Seitz 1940, Table III). The magnitude of o,(w) is deter- 
mined by |V',y,| and.the above values of the latter were chosen so as to obtain the 
best fit between the theoretical curves and the points which were determined from 
experimental data. 

The question arises as to whether (27) will give the intrinsic absorbing power 
of the perfect metal, with which we are concerned, or that of a distorted and 
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Absorbing powers of the alkali metals due to photo-excitation. (a) Sodium. (6) Potassium. 
(c) Rubidium. (d) Caesium. 


contaminated surface layer. We shall suppose that the former is the case for two 
reasons. Firstly, Ives and Briggs took great care to obtain uncontaminated 
specular surfaces without polishing or cutting. Secondly, in the case of sodium, 
we may verify from (28) and the Figure that o,(w) would be greater than the 
observed value of o(w) for low frequencies if p were increased by more than 
3-63 x 10-¢ohmcem. Such an increase would not alter the position of the points 
which were determined from the experimental data for higher frequencies 
sufficiently to affect the appropriate value of |V./.| very much. Further, in the 
case of potassium, rubidium and caesium, such an increase of p would not seriously 
alter the position of the points which were determined from experimental data 

The predicted absorption edge of sodium coincides with that observed as 
indicated by the circles. Unfortunately, the optical constants of the other alkali 
metals have not been observed for sufficiently long wavelengths to indicate the 
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position of their inner photoelectric absorption limits. The agreement between 
the theoretical curves and the points determined from experimental data is fair for 
sodium, potassium and rubidium, but there is little agreement for caesium. Since 
the assumed value of |V'./2| is about the same for each of the former three metals, it 
appears that the differences between their absorbing powers are due largely to the 
differences between their lattice constants. Further, it appears that the nearly 
free electron model is inadequate for the discussion of the absorption: of light by 
caesium. 
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ABSTRACT. A method previously used for ‘ distorting’ the wave functions of a simple 
harmonic oscillator to fit the potential curve for a particular electronic state has been 
modified, giving wave functions for the vibrational states which are close approximations to 
those of Morse, with comparatively little labour. The transition probabilities obtained 
from these functions are compared with experimental values for the first negative bands of 
N,*, the Swan bands of C,, the CN violet bands, and the a-system of BO. 


oy INTRODUCTION 

NUMBER of attempts have been made to compare the relative intensities. 

of the bands in a band-system, or the transition probabilities deduced 

from these, with the corresponding values predicted theoretically, making 

various assumptions as to the form of the vibrational wave functions %,, 4, and 
taking the transition probability as proportional to [[, #.dr]?. Since many 
diatomic molecules have vibrational energy terms of a form approximating 
closely to w(v+4)—xw,(v+4)?, with the usual notation, the Morse potential 
function and the corresponding wave functions are likely to give values for the 
transition probabilities in good agreement with those obtained from intensity 
measurements, and this has been found to be so where the test has been made, 
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But the length of the calculations involved, which must be made independently 
for every electronic state considered, prevents extensive use of this method, 
limiting it to a few bands of low vibrational quantum numbers. In some cases 
the probabilities have been calculated assuming simple harmonic vibrations. 
For molecules which are symmetrical, or nearly so, this can be done by the method 
of Hutchisson (1930), and for any molecule the wave functions for such motion 
can be plotted from the tables given by Gaydon and Pearse (1939) and by Pillow 
(1949, 1950), and the integration carried out numerically. Here, as might be 
expected, there is fair agreement with experiment where intense bands occur at 
low vibrational levels for both states, that is, where the equilibrium internuclear 
distance is nearly the same for the two electronic states; but the departures 
from agreement become rapidly more marked as the quantum numbers increase 
above | or 2. 

Between these methods lie various ways of distorting the simple harmonic 
wave functions to fit potential curves derived from band measurements, such as 
those used by Gaydon and Pearse (1939) and Pillow (1949, 1950). Both of these 
methods are open to the objection that the central zero or turning point of each 
wave function is fixed at the equilibrium internuclear distance 7,. ‘This is not 
so for the Morse wave function, and it is almost certainly not true of the vibrations 
actually occurring. An attempt is made here to improve the ‘ distortion’ process 
in this respect. So far, the method has been applied to systems for which calcula- 
tions had already been made using Morse or simple harmonic functions, so that 
comparison could be made with these as well as with experiment. It will be seen 
that the new method gives results agreeing with experiment at least as well as 
the Morse function over this limited range, and rather better than the other 
methods used, and that the form of the wave function obtained is a close 
approximation to that of Morse. 


S27 METHOD OF DISTORTING WAVE KUNCTIONS 


Method A. 

In Figure 1 the dotted curve APE represents the potential function for simple 
harmonic vibrations, U=Da?é?, and the full curve BPF represents the Morse 
potential function U=D(1—e-“)?, where D=w,2/4xw,, a? =87?ucxw,/h, 
&=7—7,. For any chosen value of U, write é;,=AC=CE, &,=CF, & =BC, 
taking the numerical value of each length, without regard to the direction in 
which it is measured. ‘The line PQ represents the locus of the central zeros 
or turning points of the Morse wave functions, and it can be shown that on this 
locus the abscissa increases by x/a for each vibrational level, provided that 
(vx)? is small compared with vx. (x is generally a fraction of 0-01.) PQ is 
therefore nearly, but not quite, a straight line. 

Now U=Da%é,?2, and £,, —&, are the roots of U=D(1—e-“)? and can be 
written —a@€=In[1 + (U/D)"?] =In(1 + a€,), so that €, = —2-303a 1 logy (1 —a@g)), 
€,=2:303a-1 logs (1+.2£,). 

Write CD=z. Then dz/dv=x/a; but dU/dv=w, so that dz/AU=x/aw,, 
and z=0 for U=0, thus z=xU/aw, = U/4Da=aé,?. Hence in Figure 1 


DF = —2-303a logyy (1 — a&) — 4a), 
BD= 2:303a~1logy9(1 + 4&9) + fa€,?. 
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The simple harmonic wave functions are tabulated already for ¢ integral values 
of &)/d where d?= h/472 cuw,. The plot of such a function is shown in Figure 2. 
In plotting the distorted function shown in Figure 3 the ‘centre’ D is placed at 
approximately the Morse ‘centre’ for the corresponding level, i.e. at 
r=r,+(v+1)x/a. Now lay out the abscissae along the axis so that if, for instance, 
the point €)/d=1-5 originally lay at a distance CE (Figure 1) from 0, it now lies 


Figure 2. 


Figure 3. 


at DF, while &)/d= —15 lies at BD behind 0. On these abscissae plot the 
canned as for the simple harmonic wave function, and so obtain the distorted 
graph. é 

.The previous distortion method consisted essentially in the use of CF, BC; 
where now DF, BD are used, and in practice these distances were read off from 
the graph of the potential curve instead of being calculated.) 
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Method B. 

Another method of distortion tried, but not found very successful, is as 
follows: For each vibrational level, a value 7, can be obtained from the rotational 
constant B, for that state, using the relation B,=A/8r2 cur,?. 7, might then be 
considered as a ‘mean’ internuclear distance for that level, and as the ‘centre’ 
for the wave function. A locus similar to PQ in Figure 1 was plotted using 
these values. The distortion was then carried out as described above, except 
that the distances were read off from a potential curve graph instead of being 
calculated. In the-case tried, that of N,+, the slope of PQ was considerably 
greater when found in this way, and PQ roughly bisected such lines as BF. 
There is much doubt as to whether the function of the internuclear distance 
which is averaged by this means is the appropriate one for the purpose, and 
comparison of the results with experiment did not seem to justify investigating 
the method further. 


§3. DIRECT COMPARISON OF WAVE FUNCTIONS 


It was considered interesting to compare the actual form of the wave functions 
obtained from the Morse expression with those obtained by the present method 
of distortion (Methea A). This has been done, using first a state in which the 
anharmonicity of the vibrations is large, and the approximations involved in 
the distortion method are therefore considerable. Dr. K. E. Shuler, of the 
Johns Hopkins University, kindly furnished the figures from which he had 
already plotted the Morse functions for the v’=0, 1, 2 levels of the excited 
*>+ state for OH, and from these the Morse wave functions have been plotted, 
together with the distorted form. It can be seen that even here the agreement 
is close (Figure 4). 

we 
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Figure 4. Wave functions for OH excited Figure 5. Wave functions for N,+ ground 
state :—Morse function; - - - ‘ distorted ’ state :—Morse function; - - = ‘ distorted’ 
function. function. 

A more favourable state is the ?&,+ ground state of N.t. The Morse functions 
for v"=0, 1, 2 have been calculated, and plotted together with the distorted 
functions (Figure 5), and the diagrams show how very close is the agreement. 
Comparison has also been made in the same way for the excited state of N,*, 
and the agreement here is equally close*. 


* From figures kindly supplied by Mr. D. R. Bates. 
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§4. COMPARISON OF RESULTS FOR TRANSITION 
PROBABILITIES OR INTENSITIES 

In Tabie 1, a number of results for the first negative N+ system are shown. 
The experimental values are estimates of transition probabilities made from 
intensity measurements due to Herzberg (1928), and these have already been 
compared by Bates (1949) with values calculated from Morse functions, and by 
Stephenson (1951) with those from simple harmonic functions. Calculations 
have now been made by the previous approximate method, and by both the 
modified methods of distortion described here. Since Bates’ calculations were 
made, Herzberg (1950) has suggested amended vibrational constants for the 
excited state, and as the results are very sensitive to these constants the values 
for transition probabilities have been calculated from both sets, so that com- 
parison can be made both with previous calculations and with experiment. 
All are expressed on such a scale that the sum of probabilities in one complete 
progression is as nearly as possible equal to 1. 

In a recent paper Turner and Nicholls (1951) have calculated the same 
probabilities by a distortion method, and for completeness their results are 
included in the Table. They do not quote the values used for the molecular 
constants, neither do they make clear whether the method used is that of Gaydon 
and Pearse (1939) or that of Pillow (1949). These are not quite the same, and 
would probably give somewhat different results. 


Table 1. First Negative N,+ Bands. ‘Transition Probabilities 


B2X,,* (1) we’ =2419-84, Rais I). Tete OURS 
(2) we =2419-84, KO — Doe 19N ey = MNO TS 
Renee ays — PAYNE), NO Onl. (ay = OMAR Ne 
Experiment Morse Distorted functions Seyldls 
Ui aN (H) (B) IND) AN@AY Oil (AD ea IND). (iB) functions 
(S) 
0,0 3914 0-69 0-65 0°65 0-67.) 0765) 10270 0:65 0-69 
0,1 4278 0:26 0:26 0-23 0:23 Se O=li/e 022 0:37 0:23 
0,2 4709 0-04 0:07 0:08 0:08 0:04 0-05 0-13 — 
i OMS OS2, 0:23 0-30 O22 95025 e ee OLS ()e3533 0-30 0-29 
1,1 3884 0:27 O22 OL2 6505255 0-24 027 0-19 0-26 
1,2 4236 0-38 0-29 O24 0260-18 O77, 0-38 0-41 
1,3 4652 0-11 — —~ 0-14 — 0-07 — == 
2,0 3308 0:07 0-05 0:05 0:04 0:04 0:06 0-04 0-02 
2,1 3564 0-39 0-41 0:38 0-42 0:38 0:43 0-34 0-26 
PoP, SI5ss 0-07 0-04 0:08 0:07 0:08 0:08 0-03 — 
2,3 4199 0-25 — OOF — 0-24 — = 
2,4 4600 0-18 — — 0-20 _ 0-09 — == 
3,0 3078 0:00 0-00 0:00 0-00 0:00 0-00 0-00 0-00 
Boil SA, 0-08 0-11 0-11 0-11 0-08 0-15 0-08 = 
ye Siew, 0:36 0-41 0-45 0-42 0:46 0-44 0-31 0-29 
BS) 5855 0-03 — — 0-03 _- 0-01 — == 
3,4 4167 0:20 — — 0-14 = 0-17 — = 
3,5 4554 0:27 — — 0-17 ~ 0-08 = = 


(H), Herzberg; (B), Bates; A, B, methods A, B; (S), Stephenson; (T & N), Turner 
and Nicholls. 
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Table 2 shows a similar comparison made for the C, Swan bands. The Morse 
function has not been applied to these bands, but two pairs of workers, McKellar 
and Buscombe (1948) and Patel and 'Tawde (1950) have made estimates of 
transition probabilities from simple harmonic wave functions. ‘The experimental 
values quoted here were obtained by King (1948) and by Johnson and 'Tawde 
(1932) from their intensity measurements. In each case the temperature of the 
source was estimated by the workers concerned, and used in reducing the values 
to probabilities. Comparison should therefore be valid for ratios in one 
v” progression, but there is room for a constant factor of error in passing from 
one progression to another. 


Table 2. C, Swan Bands. Transition Probabilities 


AIT, De —17/88°22" XWe’ = 16°44, Te =1-266 A. 
ye SIDE, We” = 1641-35, Oe — AO ie SUSI A, 
Distorted Simple harmonic 

DO" eK Experimental values functions functions 

(K) | (J & T) A Old (McK & B) (P & T) 
ORO a5 1651 0-670-540-4958 0"6150:53 7 0°6705 0°77. 0-83 0-72 0-73 
OS OS DEO Om O>2 G0 27022 2a 02285 0-29 ee O-205 OI 0:25 0-24 
0,2 6191 0:03 0-23 0-24 0-17 0-19 0:04 0-05 0-02 0-04 0-03 
Ones e245 0-228 0270-206 Ont n0:-405 0223" 50:38 0-23 0-23 
Lele Sit 98 0-35. 10-2950'20)10:2450°2470:22. 5 0-38) | 0-40 0-31 0:34 
ieee SOO 568022300255 0-198 0-250-200-2802 1 0377 0:36 
tees 20055 02002298 0rs7 Oss a0 a Ostd O04 0-09 0:07 
2,0 4382 0:03 0:02 0-01 0:03 0:02 0:02 0:02 0-07 0:02 0-04 
Dele -5 502558 0-36050:3080-159 0-455 9 0-38 0-55 0:30 0-30 
D2 S097 Osx WxNer O20 Hoy, WHO) Cable) 9 Ries 8278) 0-10 0:12 
DS Sey Osis Wess Ooms Weiley Weve Oongr -— {iavey (oat 0-39 0:39 
2,4, 6059 0-13 0-26 0:33 0-42 0-53 0-20 0-18 0-08 0:14 0-13 
3,1 4371 0:07 0:04 0:04 0:08 0:02 0:09 0-08 0-17 0-11 0-10 
3,2 4697 0-48-0°32 0-33 0-41 0:18 0:50 0:44 0:56 0-28 0-28 
3,3 5070 — 001 — — 0:04 — OHO; Oeil 0-02 — 
3,4—0501 40-53 0°26 0°27 0:24 0°24 0:24 0:22- 0-20 0:36 0:36 
4,2 4365 0°17 0-09 0-04 0-23 0:02 0-13 0-13 a= 0-15 0-15 
4,3 


4684 0:63 0:37 0:36 0:50 0:20 0:44 0-46 — 0:24 0:24 


(K), King; (J & T), Johnson and Tawde; A, method A; (McK & B), McKellar and 
Buscombe; (P & T), Patel and Tawde. 


Table 3 gives the same comparison for the CN violet bands, the experimental 
values being due to Ornstein and Brinkman (1931) and to Tawde (1936). The 
calculations from simple harmonic wave functions were made by McKellar 
and Buscombe (1948) and by Patel and Tawde (1950). 

Table 4 shows values for the «-system of BO. Elliott (1933) set out to compare 
the results of calculations from both Morse and simple harmonic functions 
with those of experiments made by himself on this system, which he chose 
because there is very little overlapping of bands. His experimental values for 
the intensities are the result of integrating over the whole of each band. For 
the results of his calculations he quotes intensities at infinite temperature instead 
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of transition probabilities, and since he has not estimated the temperature of his 
source the comparison with experiment can be made only within one v” 
progression. 


Table 3. CN Violet Bands. ‘Transition Probabilities 


B2+ We =2164-15 KOei— 20/2255 te 1150 AL 
repay we” =2068-7 Oooo (ie NON DIN 
v,v’ 2 Experimental value Distorted functions S. H. functions 
(O & B) (T) A Old (McK & B) 
0,0 3883 0:91 0-84 0:91 0-91 0-90 
0,1 4216 0-05 0-13 0-09 0-03 0-09 
0,2 4606 0-00 0-02 0-02 0-00 0-00 
LO 30 0-15 0-09 0-09 0-13 0-08 
i orf | 0:79 0-68 0:77 0-81 0-73 
1,2 4197 0-09 0-21 0-12 0-06 0-15 
1,3 4578 0-01 0:05 0-02 0-00 0-01 
PN BESS) 0:26 O17 0-16 0-19 0-14 
2,2 3861 0-70 0-52 0-71 0-71 0-60 
2,3 4180 0-11 0:27 0-14 0-08 O22 
2,4 4553 0-01 0-05 0-06 0-00 0-01 
3,2 3584 0:22 0-13 0:23 0-31 0:19 
355 Ss 0-70 0:52 0-65 0:58 0-49 
3,4 4167 0-14 0-38 0-15 0-11 (0:27 
Art esosS 0-18 — 0-24 0-41 0-23 
4,5. 4158 0-17 0:47 0-15 0-09 0:31 
De Oma oD 0-23 —- 0-15 0-11 0-51 


(O & B), Ornstein and Brinkman; (T), Tawde; A, Method A; (McK & B), McKellar 
and Buscombe. 


Table 4. «-Bands of BO. Intensities 


ACen Oe —1125071k BYES = NDCC, Te =1-348 a. 
EDI @e 1885-45 Noe alesis Te OO 
UU" A Experiment Morse Distorted Sain 
(E) (E) functions functions 
(E) 
0,0 4248 40 46 30 21 40 
0,1 4613 96 102 88 69 WS 
0,2 5040 84 94 WS 81 80 
0,3 5548 64 80 55 66 61 
0,4 6160 34 41 Not’calc. 31 37 
1,0 4036 88 92 73 85 110 
1 143652100 100 100 100 100 
1,2 4744 49 49 39 we 42 
2,0 3847 — 158 82 158 152 
2,1 4143 45 85 41 7 36 
2,2 4485 0 0 v 0 0 
IRS A bateys, 26 56 17 sil LZ 
2.4.0 “S3u 32 81 Not calc. 24 13 


(E), Elliott. 
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The scale here is so chosen that the intensity is 100 for the (1,1) band. We 
should expect therefore that in the v’ =0 progression the calculated values would 
be lower than the experimental, and that in the v’ =2 progression they would be 
higher, and this is roughly true. 


§5. CONCLUSION 

The ‘distortion’ method of plotting wave functions requires the preliminary 
tabulation of the simple harmonic functions; this has already been done for 
levels up to v=16, .and the values are applicable to any molecule. The 
calculation required to distort one wave function takes about the same time as 
that for the Morse function of zero quantum number; and when this has been 
done any other function in the range of levels covered by the table can be 
plotted in a few minutes, whereas each Morse function.requires a separate 
evaluation, the complexity of the working increasing with v. The results here 
suggest that at least for fairly low values of vw the additional labour is not justified. 
Any loss of accuracy involved in the numerical integration is probably within the 
limits of other errors in the method. 

While the older method of distortion could cause large departures from 
normalization, the functions obtained by the new method are more nearly 
normalized, at least for low values of v. 

So far, the modified method of distortion has been applied only to low 
vibrational levels, with the deliberate purpose of comparison with Morse and 
simple harmonic functions, and in each case a Morse potential form has been 
assumed. The method has still to be applied, with any necessary modifications, 
to higher levels, and to electronic states for which the dissociation energy is 
known to differ widely from the Morse value, w,?/4xw,. This is being considered 
for some of the oxygen systems, which had previously been investigated by the 
earlier distortion process. 

It must be noted that none of the methods discussed here takes any account 
of the interaction of rotational and vibrational energies. ‘The method of 
estimating the intensity of a band experimentally is likely therefore to affect to 
some extent the agreement between theoretical and experimental results. ‘The 
comparisons are strictly valid only if made between corresponding lines in the 
bands, or for the integrated intensities over whole bands, which is possible only 
where there is very little overlapping of bands. In addition, the effect of 
temperature must be accurately calculated before the comparison can be made. 
The experimental results collected here show that the values obtained by different 
workers can differ at least as widely from each other as from the calculated values, 


ACKNOWLEDGMENTS 

The author’s thanks are due to Dr. R. W. B. Pearse for his suggestions and 
help, to Mr. D. R. Bates and Dr. K. E. Shuler for the use of their calculated 
values, and to Dr. G. Stephenson for discussion of his own methods and results 
before publication. 

Note added in proof. ‘This newer method is an improvement on older ones, 
particularly for low vibrational levels, but still has the defect of giving excessive 
values to the wave functions for large internuclear distances, thus making the 
functions far from orthonormal when moderately high vibrational levels are 
considered, and limiting the usefulness of the method. 
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A further modification, which gives promise of greatly improved results in 
these respects, at the cost of a slight addition to the calculation, is indicated shortly 
here, using the same notation as above. There the distortion is accomplished by 
substituting abscissae BD, DF for AC, CE, both taken at the same energy level U. 
But on the simple harmonic potential curve the vibrational levels correspond to 
U=wa,(v+4), while on the Morse curve, still used as the most useful approxi- 
mation to the actual curve for the state, the corresponding level is 
U' = 0,(0+ 4) —xp02,(0+ 4). | 

The distortion here suggested is the substitution of B’D’, D’F’ at the level U’ 
for AC, CE at level U. A calculation similar to that used before then gives 


223 
BID! = = log {1 + afo(1 —Ja%dg?)!9} + fady?(1 — 20%") (1 — 16), 
"7303 agi a 
D'F’ = — =~ log {1 ~ af(1 — 2aGy?)!} — bago%(1 — 4a°E,2)/(1 — 40°E,%) 


where &,=AC=CE as before. This rather more elaborate formula does not 
increase to any serious extent the work involved. 

Results so far obtained indicate that the functions obtained in this way are more 
satisfactory than the other approximations so far tried. 
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ABSTRACT. Some theorems concerning the free energy of the electric double layer of a 
colloidal particle in a dilute sol are developed by making use of the Poisson—Boltzmann 
equation. It is shown that this free energy can be derived by employing the fictitious 
process of charging the excess ions in the double layer and that the surface density of ions 
may be an arbitrary function of the degree of charging. By integrating with respect to 
the charging parameter at constant surface density of ions we obtain an alternative form 
for the free energy which is expressed in terms of the electrostatic energy of the double layer 
and the excess osmotic pressure in the diffuse layer. The latter expression for the energy 
can also be derived by attributing to each volume element in the diffuse layer an energy 
term due to the thermal motion of the ions and an (electrostatic) field energy term. 


Si INTRODUCTION 


N a series of previous papers (Levine 1946, 1948, 1950, 1951 a,b, hereafter 
] to be referred to as I-V), the author derived a number of theorems relating 

to the free energy of the double layers of colloidal particles. We have found 
since that there exist other theorems of a similar nature, and it is the purpose of 
this paper to present these in a unified form. We shall consider the case of a 
single colloidal particle of general shape immersed in a very large (infinite) volume 
of electrolyte and make use of the Poisson—Boltzmann equation. The extension 
to more than one particle is fairly obvious. Also, we shall neglect the contribution 
to the free energy of our system from the Coulomb interaction of the ions in the 
interior of the dispersion medium (see V). 

Three main results will be derived in this paper. The first, which is given 
by equation (12), is a generalization of a theorem obtained by Casimir (1948). 
Equation (12) states that the use of a Lippmann equation in deriving the free 
energy of the double layer is equivalent to imagining that all the ions in the 
colloidal system are subject to a charging process which is carried out in such 
a way that the density of surface ions varies quite arbitrarily with the degree of 
charging. The second result was first obtained by Verwey and Overbeek 
(1948, p. 78) in the. particular case of two parallel plates. The electrical free 
energy of the double layer can be determined by charging under the condition 
that the density of surface ions remains fixed. It is possible to transform this 
expression into one which involves the electrostatic energy of the double layer 
and the integral of the (excess) osmotic pressure over the volume of the diffuse 
outer layer; this theorem is expressed by equation (14). The third result is an 
alternative derivation of the theorem just referred to, using an entirely different 
approach, which throws further light on the nature of the approximations 
introduced in the Poisson—Boltzmann equation. Associated with each volume 
element in the diffuse outer layer we have an osmotic term and an electrical term 
(namely the ordinary field energy in electrostatics); these are the last two terms 
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on the right of equation (26). A similar expression has been used by Debye (1927) 
in his theory of ‘ salting-in’ effect, and also appears in the work of Frumkin (1943) 
(mentioned by Frenkel 1946), Derjaguin (1940) and Coolidge and Juda (1946). 

The question arises as to the value of the various theorems put forward. 
As yet, we have not applied equation (12) to any practical problem, and it may be 
regarded as further evidence for the self-consistency of the general theory. 
However, the other two results are of immediate use. ‘The equation (14) 
(or rather the equivalent first relation in (21)) has been applied in a forthcoming 
paper by Levine and Suddaby to determine the free energy of two parallel plates 
immersed in a 1-2 valency type of electrolyte. Also, work is now in progress on 
the numerical evaluation of the free energy of two spherical particles, and here 
again the right-hand side of (14) is found to be a suitable form for the energy. 
Finally, the expression (16) for the free energy differs from the other forms 
obtained, both in this paper and the earlier ones, in that it is equally applicable 
to either concentrated or dilute colloidal systems; at present it is being used by 
the author to study concentrated sols. ‘The equivalence of the various charging 
processes proposed by Derjaguin (1940), Verwey and Overbeek (1948) and the 
author is based on the assumption that the volume of electrolyte is infinite, and 
a number of difficulties, which remain to be resolved, are apparently encountered 
when extending these processes to finite volumes. 

Various authors have attempted to improve upon the Poisson—Boltzmann 
equation when applied to the electric double layer. Among the types of 
deviations there are: (i) fluctuation terms (Casimir 1944, Loeb 1951), (11) effect 
of finite size of hydrated ions (Grimley and Mott 1947, Grimley 1950), 
(ii) dielectric saturation * (Grahame 1950), (iv) electrostrictive effects (Mueller 
1935) and (v) image forces (Loeb 1951). These extensions are undoubtedly 
desirable but, in the opinion of the author, a satisfactory analysis of the corrections 
to the Poisson—Boltzmann equation entails the simultaneous treatment of all 
the above effects, particularly since they partially compensate one another. 
The Poisson—Boltzmann equation can be considered a useful first approximation 
and the basis from which a more exact equation can be developed. It is likely 
that considerably more experimentation would be required on the stability 
properties of colloidal systems, accompanied by a detailed application to special 
cases of the general theorems developed here and in the earlier papers, before 
we can establish the precise limitations of the Poisson—Boltzmann equation. 
This point of view receives support from the fact that at colloidal surfaces we may 
have such complications as the formation of the Stern layer of counter ions 
(Levine 1951c) and the distribution of charge in the crystal lattice of the 
colloidal particle (Grimley and Mott 1947, Grimley 1950). . Finally, it should be 
mentioned that equation (12) and the left-hand member of (14) may still be valid 
with the ‘corrected’ Poisson—Boltzmann equation, whereas this is not true of 
either the right-hand member of (14) or of the expression (26) for the free energy. 


§2. FIRST FORM FOR ELECTRICAL FREE ENERGY 


It is assumed that the surface charge on the colloidal particle arises from the 
adsorption of 7 ions of type 1, charge e,, and that v is the density of such ions at 
position S on the surface, which has an area S. The volume density of ions of 
species 7 (¢=1, ..., s) in the interior of the solution is assumed constant 

* More recently Booth (1951) and Conway et al. (1951). 
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(independent of ) and denoted by 7°. Then, according to the analysis in the 
earlier papers (see, in particular, III and IV), the part of the free energy of our 
system which depends on 7 is written as 


UJ AVAL a) ee ees (1) 
F(n)=2 te EQ,n) d=, | BG ila edn tiga init eons (2) 
E(A, n) = LAemb(A, n) +4 is bp dv = = le [eradur|? dujeey ee (3) 
Hn) == [wraS= | dr(SdS, ees (4) 


where A denotes the degree of charging, f,(”) is the free energy (depending on z) 
of our system in the completely discharged state A=0, —Ay,(1) is the change in 
the electrical self-energy of an ion of type 1 when it is transferred from the 
solution to the particle surface, F’,(n) is the so-called electrical part of the free 
energy associated with the double layer and is defined by the first equation in (2), 
the quantity E(A,z) in the integrand being the average electrostatic energy of 
the double layer at degree A. We imagine that the integration with respect to A is 
carried out subject to the condition of a fixed distribution of surface ions. As in 
III and IV, we write this distribution as v=1(S)=nv)(S), where v9(S) is subject 


to the normalization condition | vo(S)dS = 1, but is otherwise an arbitrary function 
Jig 


of position S. In the surface integral in (4) % is the potential at position S (and 
stage A) on the particle surface. In the volume integrals in (3), % and p are the 
corresponding potential and charge density respectively at position r in the 
diffuse outer layer; these integrals are taken over the volume V of the outer layer, 
which can be assumed to extend to infinity, where the potential % vanishes. 
aba, n) is the average potential (at stage A) over the surface, the weight function 
being »,(S). D is the dielectric constant of the dispersion medium. 

At stage A and position r the volume density of ions is given by the 
Boltzmann distribution law 


Beet A, 2S)” exp (ACA RL ler ons eins mmasraroiee (5) 
where Fk is Boltzmann’s constant and T the absolute temperature, and the 
Poisson—Boltzmann equation for the potential % is expressed as 


Ve arp DENN ee (6) 
f(A) = - 2 i WOOD Wea ACH Rule — 8S © bates (7) 


Inside the particle the potential, which we denote by y", satisfies Laplace’s 

equation since there is no charge. At the surface we have the boundary conditions 
os oyp" 
ip! ts ET) Ugseiead ar e* 
p=, D an, D an, ANC ae I On aaa (8) 
where n, denotes the outward normal at the particle surface into the region of 
integration and D! is the dielectric constant of the colloidal particle. It should 
be observed that the second equation in (2), which was derived in III, is valid 
for an arbitrary function f(A), provided that the latter vanishes at infinity, 
8 

where # is zero; we note that & n°¢,=0. 


t=1 
51-2 
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§3. GENERALIZATION OF CASIMIR EQUIVALENCE THEOREM 


Instead of charging under the condition of fixed density of adsorbed ions, 
we shall imagine that this density varies quite arbitrarily with A. For simplicity 
we first consider the case of uniform density of surface ions (v)(S) = 1) which would 
apply to a sphere, or to a section of a cylinder of infinite length or of an infinite 
plate. Then the function »=n(A) will represent the variation of the number of 
surface ions with the stage A. Also (A,7) will denote the (uniform) potential 
at the surface, a function of A and ; this potential now replaces ¥(A,7) in the 
equations of the preceding section. 

We begin by equating the first expression in (3) for E(A,7) to an alternative 
formula which was obtained in III (equation (30)), namely, 


EQ,n)=Pe, | Lia Iede eee (9) 


where the subscript A indicates partial differentiation, 1.e. subject to fixed n. 
Since the relations (3) and (9) are identities in A and 2, we can choose n=n(A), 
which, therefore, becomes the upper limit of integration in (9). We ee derive 
now that 


1, vlA, n(2)) nd 
Je pi herds = =e), n(A, Ws) dys + Ae, ls den dneaaae (10) 
where the function m =n(A, #) is the inverse to F=*(A, nm). But 
n(2) d via, nay) 
b,(A, n) dn = ary n(A, w)dyb + n(A)b’[A, n(A)], 2. ee (11) 


after integrating by parts; the prime denotes total differentiation with respect 
to A. If we substitute (11) into the right-hand side of (10) and then integrate (10) 
with respect to A, we obtain 


vil, n(1)] 
~ 0 


n(A)] 


n(1, pd +e, lim A Ie n(A, Wi) deb 


dn 
S -| ON iE [pln —nay do — ey (h An(A)p'[A, n(A)]drA. «2... (12) 


The first term on the left of (12) can be obtained by integrating the well-known 
Lippmann equation from electrocapillarity theory. It occurs in the expression 
U(n) for the free energy in the following way. The number of surface ions is 
not arbitrary but is determined by the condition of thermodynamic equilibrium, 
which is that U(m) is a minimum, 1.e. U’(n) =0. If% is the equilibrium value of n, 
we readily derive from (1) and (2) 


U(%) =fy(") —f,'(R) — ie Les er (13) 


If we identify 7 with n(1) then the integral term on the right of (13) is given by (12). 
Now the resulting form for U(n) can also be obtained from thermodynamic, 
statistical reasoning by imagining the charging process carried out along the 
path =n(A), the end point being A=1, m(1)=%. In this thermodynamic proof, 
the details of which will not be given here, the potential of the mean force at an 
adsorbed ion is identified with the corresponding mean potential, and the same 
potential function is assumed for the different types of ions in the diffuse layer. 
In view of this alternative derivation of (12) it might appear that we have shown 
the Poisson—Boltzmann equation to be a good approximation yielding reliable 
estimates of the free energy of the double layer. It should be pointed out, 
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however, that (12), like the second relation in (2), is based on the use of (6) and is 
valid for any form of f(A’). The interesting feature is that as we choose different 
functions for f(A’), the two members of the equation (12) will have different 
values, but they will always be equal to one another. ‘Thus the relation (12) is 
less specific than the Poisson—Boltzmann equation, and it seems that no 
conclusions concerning the reliability of the latter can be reached in this way. 

If we can choose a path m=n,(A) such that (i) the surface potential 
BA, n(A)] =%, a constant independent of A, and (ii) the second term on the left 
of (12) vanishes, then (12) reduces to the theorem of equivalence proved by 
Casimir (1948). ‘These two conditions are satisfied in the case of an infinite 
plate-like particle or a cylindrical particle of infinite length, even though in the 
latter case m,(0) is infinite. However, for a spherical particle (and probably for 
any particle of finite extent) m)(0) is again infinite, and the second term on the left 
of (12) no longer vanishes (see III). 

The result (12) can be generalized to apply to a particle of general shape, the 
surface of which is at a uniform potential, say % (A), which is an arbitrary function 
of A. The analysis is quite similar to that given at the end of §4 in III, where we 
considered a more particular case, namely (A) =, a constant independent of A. 
We first choose an arbitrary path of integration n=n(A). ‘Then at each stage A 
the surface ions are so distributed over the particle that the latter is at a uniform 
potential %(A), the value of which depends therefore on n(A). In this way we 
obtain (12), provided the function yA, n(A)] is replaced by (A). 


$4. SECOND FORM FOR ELECTRICAL FREE ENERGY 


The following theorem will be proved: Let Ap(r) denote the difference between 
the ideal osmotic pressure of the ions at position r in the diffuse layer and that in 
the interior of the solution ; then the electrical free energy F’,() can be expressed as 


es | v1, n)dn = — \, Ap(r)dv—E(1,n) +enW(l, n). .....- (14) 


Whereas the left-hand member is valid for arbitrary f(A), the right-hand 
member is-derived by using the Boltzmann form (7) for f(A). If we substitute 
the first equation 1n (3) into (2), then we have two terms in the resulting expression 
for F.(n). ‘The first term reads if 


es I dn | peas, ae (15) 


and, the second term can be written as 


[12] doden 2 fin srncnnees 2 fa), wre 


To obtain (16), we write %=(A‘s),—A%s,, where the subscript A again denotes 
partial differentiation, and make use of the first equation in (6). It is understood 
that the integration in (15) and (16) is carried out subject to fixed surface 
distribution of ions, so that v,=0. On introducing the formula (7) for (Ad), 
interchanging the order of integration and making use of (5), the first term on 
the right of (16) reduces to 


kT | @l, )—nPjdo=— | Apleydo, oes. (17) 
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It is assumed here that A’—+0 as A+0, which is true for a single particle immersed 
in an infinite volume of electrolyte. The second term in the right-hand member 
of (16) can be expressed as 


- {OP ies =| ace a 
= -2[. [pVv* Verdot > i an | ee si —, sa as. sree (18) 


The surface integral on the right is obtained by niger Green’s theorem, 
applied to the functions and %,. We note here that when A=0, % =0, or possibly 
y=constant, throughout the dispersion medium, in which case [V7b],9=0; 
the latter applies to an infinite plate. The second term on the right of (18) can 


be written as at 
he, | dd { (he Vase ee (19) 
0 S 


substituting the boundary conditions (8). To derive (19) we also make use of 
the relation v,=0 and of the result that the surface integral on the right of (18) 
vanishes when we replace the function / by %'. The latter follows on transforming 
the surface integral by means of Green’s theorem into a volume integral taken 
inside the particle itself and applying the conditions that both %' and %,' satisfy 
Laplace’s Saas If we now add (15) and (19), we obtain 


oe ja ay (ub),vdS =4e | Whew dS=demp(l,m), (20) 
integrating by parts swith, respect to A and putting v,=0. The free energy F’,(n) 
consists of three terms, namely, the quantities (17) and (20) and the volume 
integral on the right of (18). By making use of (3) it is readily verified that the 
sum of these three terms is identical with the right-hand member of (14). The 
formula (14) will be applied in later papers to determine the interaction of colloidal 
particles. 

If the particle is metallic, the preceding analysis must be modified, but since 
the details are quite similar to those described in paper V, we shall omit them. 
Although v, may no longer vanish, it can be verified that there are no additional 


terms in (19) and (20), since % is constant over the surface and n=| v dS is 


independent of A. It follows that for a metallic particle the potential (1, 7) in (14) 
is replaced by (1, 7). 

Verwey and Overbeek (1948) considered the case of a dielectric particle for 
which the potential is uniform over the surface. Thus, if the particle is an ionic 
crystal and an adsorbed ion is identical with one of the lattice ions, then under 
thermodynamic equilibrium no (isothermal, reversible) work is required to 
transfer such an ion from one part of the surface to another. If we now assume 
that both the non-electrical (chemical) energy of adsorption and the electrical 
self-energy change are constant over the surface of the particle (which is of general 
shape), it follows that the particle surface is at uniform potential. (We identify 
the electrical potential of the mean force at an adsorbed ion with the corresponding 
mean potential and ignore possible surface imperfections.) In these circumstances 
the charging process which yields the second relation in (2) is carried out under 
the condition of fixed surface density of ions v=nv,(S,n), where the function 
vo(S,7) is so chosen that the potential is uniform over the surface at the final 
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stage. We note that because of the non-linearity of the Poisson—Boltzmann 
equation v9(S, 2) actually varies with m; in the case of very small 2 we may apply 
the linear Debye—Hiickel equation and then »(S, 2) no longer depends on 7. 

There are differences between a metallic particle and a dielectric particle at 
uniform surface potential in the normal equilibrium state A=1. First, in the 
former (metallic) case the potential is uniform over the surface at any intermediate 
stage during the charging (owing to the presence of free electrons), even though 
the distribution of adsorbed ions is fixed, whereas in the latter case this is not true. 
Secondly, for a metallic particle the second relation in (2) is valid when (1, 7) 
is replaced by the uniform potential (1,7), the relative distribution of adsorbed 
ions vo(S) being quite arbitrary. But again this property does not apply to a 
dielectric particle of general shape, since it becomes necessary to replace v9(S) by 
the definite distribution v9(S, 7) in (4), in order to achieve a uniform potential 
#(1,) at each value of n. 

In the case of either a metallic particle or a dielectric particle of general shape 
under the conditions specified above, we may replace #(1, 2) in (14) by (1,7), 
which we denote simply by %. On integrating by parts the left-hand side of (14) 
we now obtain 


=e, | n(Wdb=— | Aple)do— EC, 


4 dx J Yo 
= | - | [pp],,dv—e, lim d[nQ Wa, Aa (21) 
0 r JV A4>0 0 


the second relation following from (12), where the subscript % denotes that the 
integration is carried out subject to the condition that the surface potential is 
uniform and given by yp, 1.e. it is independent of state A and position S. Fora 
metallic particle this condition is attained by suitable choice of the total number 
of ions as a function of A and y%. For a dielectric particle we would also need to 
specify the relative distribution of adsorbed ions as a function of A and yp. 

The equation (375) of Verwey and Overbeek (1948, p. 80), which applies to 
the case of two parallel plates, is equivalent to the second equation in (21), and 
it is possible to derive the latter by generalizing the proof of these authors. The 
first term on the right of this second equation is first written in the form (16), 
and we again obtain the formulae (17) and (18). Making use of the condition 
“sb, =0 on the surface of the particle, we interchange the order of integration in 
the second term on the right of (18) and then integrate by parts. This second 
term thus simplifies to 


a a, E |. .¢s = — 7 Nify + Fey lim An(A, Wo). .. «ess (22) 
Ns _\,=0 70 


The antes side follows on application of the boundary conditions (8), 
since %'=%, a constant. (In the case of a spherical particle the second term on 
the right of (22) does not vanish.) Now 

Ho lim An o)= lim J" nQ dat ivi (23) 


+70 
which can be shown from the theory of electrostatics as follows. When A=0, 
but g0, where q is the charge on the particle, we simply have a charged conductor 
immersed in a dielectric and qg is proportional to the surface potential, being 
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related by the coefficient of capacity. This coefficient is independent of the value 
of q, and hence (23) is valid. By substituting (22) for the second term on the 
right of (18) and making use of (23), we readily find that (16) becomes identical 
with the second relation in (21). 


§5. TWO PARALLEL PLATES 


The analysis which yields the expression on the right of (14) for F’.(m) must 
be modified to apply to sections of two parallel plates of infinite area. We shall 
consider the double layer on the inner face only of such a section, which for 
simplicity has unit area, so that we write =v in the following analysis. ‘The 
diffuse layer of this section extends from the plate face to the median plane; 
if R is the separation of the plates, the volume of the diffuse layer is R/2. If we 
charge at constant density v of surface ions, then, as already explained in the 
earlier papers (see, in particular, II), the potential over the whole of the double 
layer tends to a common infinite value as A>0, namely As—Ayj(A, v) =m, say, a 
finite quantity different from zero. For finite plates (i.e. in a real sol) m=0 and 
the origin of a non-vanishing m is therefore the incorrect use of the model of 
two infinite plates to represent two finite plates in the vicinity of A=0. As a 
consequence we have a spurious initial (or residual) electrical free energy in the 
completely discharged state, which can be expressed in one of two forms, namely, 


~e, | mdv=JkTR > np fexp (—em|kT)—1} 
¥ (0 t=! 


+imR & ne, exp (—em/kT)—4eyvm. ....... (24) 
i=1 

This can be shown in the following way. The formula (14) is always true 
for a particle of finite extent, provided that the Poisson—Boltzmann equation is 
used and the volume of electrolyte is infinite. We can also apply (14) to each of 
two identical finite parallel plates (with edges in line); the volume of electrolyte 
is divided into two by the median plane and the diffuse outer layer of each plate 
is assumed to extend over an (infinite) half-volume. Since no charging process 
is involved in the equation (14), the troublesome behaviour at A=0, described 
above, does not arise when we imagine the area of the two plates to increase 
indefinitely. Consequently (14) remains valid for each of the two unit sections 
taken from the inner faces of the infinite plates and (1,7) is now replaced by 
the uniform potential over each plate :J(1, v). But we can apply the first formula 
in (2) for F(z) =F .(v) to unit (inner) area of each infinite plate and by suitable 
transformations obtain two equivalent forms for F,(v). The first form, which 
was derived in III, is the sum of the left-hand members of (14) and (24), since 
%(1,”)=y(1,v); the second form is obtained by adding the right-hand sides 
of (14) and (24). The validity of the equation (24) follows immediately. It 
remains to indicate the derivation of the three terms on the right of (24). The 
first term must be added to (17) and accounts for a difference in the densities 
of ions between the plates and in the bulk of the electrolyte at\=0. The second 
term would occur in the right-hand member of (18) and is due to a non-vanishing 
value for {V*s=Adsf(As’) at A=0. Finally, the last term in (24) would appear 
on the extreme right of (20). 

The validity of (24) can be verified in another way. In an earlier paper 
(Levine and Suddaby i951) a series expansion for the potential in the double 
layer which applies to the case of a binary symmetrical electrolyte (7,°=7,°, 
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€,=—€,) and moderate potentials has been obtained by solving the Poisson— 
Boltzmann equation. Making use of this expansion, we find that 
el = 3 (BY mt +z (fx) met (25) 
o— 6 lpr) Mo + G0 lap) Mom eee ee e 


where m)=vkT/Re,n,°. The approximation m=my is obtained if we replace 
the Poisson—Boltzmann equation by the linear Debye—Hiickel equation. If we 
substitute (25) into (24) and then expand each side of (24) in powers of v*, we 
obtain identical coefficients of the first three powers of v®. 


§6. ALTERNATIVE DERIVATION OF SECOND FORM 
FOR FREE ENERGY 
e free energy of our system, consisting of one particle immersed in a ver 
The f gy of ystem, ting of particl d y 
large volume of electrolyte, is assumed to have the form 


F(n) = > N,€,—n&, + f,(n) + z Nixs(1) — Aya (1) 


+kT | (2 n,n n, HORE (iar) Ley (26) 
s 


Jv \i=1 
In this section all quantities refer to the actual physical state A=1. The 
quantities €,(2=1,...,.s) are functions of volume (or pressure) and temperature 
only and represent the part of the chemical potential of the ions due to the 
short-range van der Waals interaction with the surrounding solvent molecules.* 

fm) is the non-electrical part of the free energy associated with the surface 
phase of m adsorbed ions of type 1. ,(1) is the electrical self-energy of an z ion 
in the dispersion medium. ‘The last two terms on the right of (26) can be 
interpreted as follows. Consider an element of volume dv in the dispersion 
medium, containing n,dv ions of typez (z=1,...,5). The free energy attributed 
to these ions includes an osmotic term k7n,;lnn;dv and an electrostatic term 
nex dv which arises from the Coulomb interaction of the ions. Thus in the 
chemical potential expression of an z ion in the volume element dv we have the 
contribution e,s+kT Inn, which has the same value throughout the whole of 
the dispersion medium. If we compare the element dv with another dv , the 
density of 7 ions and potential in the latter being n° and 4 =0 respectively, then 
the condition of equal chemical potentials obviously yields the Boltzmann 
distribution law (5). We now integrate over the dispersion medium and _ halve 
the resulting electrostatic contribution since we have been counting the Coulomb 
interaction energy of each ion pair twice by this method. The corresponding 
electrostatic contribution (again halved) from the surface ions is finally added. 
If we use the second form in (3) for E(1,) then alternatively we can assign the 
usual form for the field energy to each volume element dv. 

If we substitute (5) for n, (after putting A=1), we derive that the ‘osmotic 
term’ on the right of (26) becomes 

8 § f 
ot i> | n, (In n?—e,h/kT)do=KT | x N, In n,9—n In nm? =| ppdv, ..(27) 

i=1/V i=1 V 

s 
since p= & ne, Applying the Boltzmann distribution law, we have 
i=] 
nev | 1+ . [fe (—ex%/kT)—1}dv | = o a5 teste 


*jInstead of 7; In n;, one should write n; (In n;—1), but the additional term is absorbed into &. 
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and thus, if the volume V is large, a very good approximation is 


1 — In{(V,-2)/V}, 7=1, 
In no+ 5 \, fexp (—ep/AT) — 1}do= hy, ae eRe SC) 
On substituting this approximate form for In n° into (27) and making use of (3), 
we find that (26) can be expressed in the form 


Fn) SF (n)* z Neuf) Aw (DEC Lee 

very nearly, where F(z) is given by the right-hand side of (14) and 
Fy(n)= E NAE,+kT In (N/V)}+ h(n), eee (30) 
mhere He Eis LE oak T ein eee (31) 


The quantity F,(n) is the free energy of our system in the completely discharged 
state. To derive the expression (31) for f(z) we have expanded F(z) in powers 
of n/N, and neglected terms of order (n/N,)?. 

The physical significance of the preceding method of evaluating the free 
energy can be described as follows. We imagine that the ions are all discharged 
and that those in the dispersion medium are uniformly distributed. Then these 
discharged ions are redistributed so as to be in their final average positions, 
corresponding to the fully charged state. The second last term on the right 
of (26) (the osmotic term) represents the work done against thermal motion in 
this redistribution. The last term in (26) (the electrostatic term) is the work 
done in charging all the ions against the potential ys, assuming the ions are fixed 
in their average final positions. This method of charging is not an isothermal, 
reversible, process, and the approximation introduced thereby is implied in the 
Poisson—Boltzmann equation itself. 
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The Conversion and Auger Electrons of Radium D 


By D. Keeur Tl snp W.. D. BRODIE 
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ABSTRACT. The conversion and Auger electrons of Radium D have been studied in a 
lens B-ray spectrometer with a post-focusing electron accelerator. The number of conver- 
sion electrons due to the 46-7 kev. y-ray is found to be 59:2+ 5 per hundred disintegrations 
but reasons are given which suggest that this value be regarded as a lower limit. The Auger 
electrons have been compared with those resulting from the Th B-ThC transition. 


Sf INTRODUCTION 

ECAUSE of the softness of the RaD radiations doubt remains regarding 
B its mode of disintegration. For many years attempts have been made to 
detect the primary f-radiation, but the problem is a difficult one since the 
end point is believed to be in the 15kev. region. The expansion chamber 
results of Richardson and Leigh-Smith (1937) are the most reliable to date, 
and these workers estimate that 60° of the B-rays have energies less than 4 kev. 
It has been suggested that the spectrum is complex, there being a second 
component with an end point of some 60 kev. (Kinsey 1948, Cranberg 1950), 

but this has not been observed. 

The presence of a 46-7 kev. y-ray has long been known both in emission and 
through the natural photoelectrons from the L, M and N shells of the residual 
atom. ‘The unconverted y-radiation is known, with fair certainty, to have an 
intensity of 3-5+0-4°% (Bramson 1930, Gray 1932, Von Droste 1933, Stahel 
1935), but the values given by various workers for the intensity of the internal 
conversion lines range from 15°% to 80° (Kikuchi 1927, Stahel 1931, Lee and 
Libby 1939, Ouang, Surugue and Tsien 1943). 

Besides the 46-7 kev. y-ray, other weaker y-rays have been discovered more 
recently (see Feather 1949) and, in all, seven y-rays are assumed to be known. 


. Their energies and intensities are listed in Table 1. Internal conversion electrons 


of one of these y-rays have been reported by Cranberg (1950). 


Table 1. The y-Rays of RaD 


Quantum energy (kev.) 46:7 42:6 Sars, Bil SaeWe3 23a ee. ierilse Mes) 9/2342 Coz! 
Emission intensity (%) 3:5 0:2+0:1 0:2+0-1 0-44+0-2 ~1 ~0 ~10 


The present investigation of the electron and f-ray spectrum of RaD was 
made with a 6-ray lens spectrometer which incorporates a post-focusing electron 
accelerator, enabling electron spectra to be studied conveniently down to 
energies as low as 4kev. A thin-windowed Geiger—Miiller counter is used as 
detector. Details of the spectrometer have already been published (Butt 1949, 
1950). 
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§2. SOURCE PREPARATION 


RaD was extracted from a supply of old glass radon needles. ‘The needles 
were boiled in chromic acid, rinsed in distilled water and finely crushed under 
about 10 ml. of concentrated HNO,. This solution was boiled and centrifuged, 
and the concentrated HNO; containing the RaD was decanted. The solution ~ 
was then evaporated to dryness and taken up again with about 2 ml. of 0:1N HNO. 
The RaD was collected on a platinum foil 0-3 cm? in area by an electrolysis in 
which the platinum foil was positive with respect to the platinum crucible 
. containing the solution. A current of 10~*amp. was passed for 10 hours. By 
this means only RaD was collected on the anode, since it alone of the elements 
present forms a complex negative ion. In particular this method of collection 
has the advantage of rejecting any traces of mercury that might have been included 
in the radon tubes. The RaD was finally transferred to a piece of gold leaf 
0:22 mg/cm? in thickness and 0-6cm. in diameter, backed, for extra strength, 
by a collodion film of weight negligible compared with that of the gold. This 
transfer was accomplished electrolytically. The platinum was mounted 
horizontally and a drop of 0-1N HNO, spread over its upper surface. ‘The 
gold leaf was then lowered until contact with the HNO, was made. A current 
of 10-tamp. was passed for one hour through the cell so formed, the gold leaf 
being the anode. From this, after washing with distilled water, the surplus 
liquid was removed with the corner of a piece of filter paper. ‘The source so 
formed was a disc of about 45mm. diameter and was a practically invisible 
deposit on the gold leaf. The strength of the source was approximately 10 uc. 


§3. THE INTERNAL CONVERSION ELECTRONS OF THE 46-7 kev- 
TRANSITION 


Figures 1 and 2 show the spectrum obtained after equilibriun between the 
RaD and RaE had been established. With maximum current in the spectro- 
meter coil it was impossible to focus electrons with momenta greater than 2,500 
gauss cm. so the upper part of the RaE spectrum was taken from the results 
of Neary (1940). This part is drawn with a dotted line in Figure 2. 

The two most prominent lines, e and f (Hp=600 and 720 gauss cm. 
respectively) are formed by internal conversion electrons due to the 46-7 kev. 
transition. Line e corresponds to conversion in the L,, L,, and Ly; shells of the 
RaE prodvct nucleus and line f to conversion in the M and N shells. Of necessity 
the resolution was too low for the lines due to the individual shells to be resolved. 

Since, for all practical purposes, the only mode of disintegration of RaE is to 
the ground state of RaF by f-emission, a comparison of the areas of lines e and f 
with the area of the RaE f-spectrum gives the fraction of the decays involving 
excited atomic levels of the RaE atom due to the internal conversion of the 
46-7 kev. y-ray. Before this comparison can be made, however, it is necessary 
to estimate the effect on the spectrum of back-scattering from the source support. 
This work falls into two parts. 


(i) Correction of the RaE B-spectrum. When compared with that of Neary 
our RaE spectrum shows deviations due‘to back-scattering at energies below 
180 kev. (Hp =1,550 gauss cm.). Accordingly the RaE spectrum was corrected 
below 180 kev. using Neary’s results, and the area under the corrected curve was 
taken to represent the intensity of the B-rays of RaE. 


The Conversion and Auger Electrons of Radium D 793 


(ii) Corrections to be applied to the area of the conversion lines e and f. ‘The 
intensity due to RaE in the region of the conversion lines was found by examining 
in the spectrometer the B-radiation from a source of pure RaE deposited electro- 
lytically on a gold foil similar to that used forthe RaD. The source was prepared 
by the method of Broda and Feather (1947). Fitting of the resulting curve to 
the previous one gave the line marked ‘RaE’ in Figure 1. This line was used 
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ds a base line when the intensity of the conversion electrons was being considered. 
To correct the area of the conversion lines for back-scattering use was made of 
the formula given by Hamilton and Gross (1950). ‘These authors state that the 
back-scattering ratio (the number of scattered electrons per primary) reaches an 
asymptotic value called the back-scattering constant S for a thickness t) of 
scattering material at an energy V, given by V,=800[(Z?/A)t)]"?, where V.. is 
in kev., Z and A are the atomic number and mass number of the scatterer and 
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t) its thickness in gm/cm®, The value of S has little energy dependence. 
Applying this equation to the source support used in the present investigation we 
obtain V,=67kev. This value appears to be a reasonable one for saturation 
back-scattering since our RaE spectrum begins to depart from that of Neary at 
about 180kev. It is nevertheless possible that the source backing was actually 
thinner than 0:22 mg/cm? owing to the solution of some of the gold in the 
0-INHNO, during the electrolysis of the RaD. However this may be, we 
assume that back-scattering has reached saturation over the whole range of energy 
(45 kev. and less) with which we are at present concerned. 

To find the value of S the following experiment was carried out. ‘Thorium B 
was deposited by recoil on 0-2 mg/cm? aluminium leaf and the spectrum in the 
neighbourhood of the F line (148kev.) was plotted. A fairly reliable estimate 
of the area of the ThB F line free from back-scattering was obtained from this 
curve. A platinum tipped brass rod was then placed immediately behind and 
in contact with the aluminium leaf and the spectrum plotted again. This gave 
the area of the F line of ThB under conditions of saturation back-scattering. 
Finally, the experiment was repeated under the original conditions to make sure 
that, after correction for the decay of the source, the third set of counts corresponded 
with the first set. A comparison of the areas of the F line of ThB in the two 
cases gave the value 0:55 +0-05 for the number of back-scattered electrons per 
primary under the solid angle conditions obtaining in our spectrometer. We 
regard this as an upper limit for our Ra D correction since, as indicated above, the 
source mount was probably thinner than 0-22 mg/cm”. 


Table 2. The Relative Intensities of the Conversion Lines of the 46-7 kev. y-Ray 
Ellis Curtis Cranberg Butt and Brodie 


ies 1-0 1-0 1-0 

ibe 0-04 0-06 0-09 +0-015 1-0 
[ies 0-01 — 0-019 + 0-004 

M 0-4 055 0-29 +0-02 } 

N 0-2 0-2 0-085+0-01 0-266 + 0-03 
fe) = 0-01 0-021+0-007 


Applying the appropriate corrections to the areas of the Ra E spectrum and of 
lines e and f, the number of conversion electrons due to the 46-7 kev. y-ray was 
calculated as 59:2+5 per hundred disintegrations, of which 46-7 +4 (i.e. 79°) 
were found to be due to conversions in the L shells. Since we have probably 
over-corrected for back scattering these values should be regarded as lower 
limits. 

Taking the intensity of the unconverted 46-7 kev. y-ray as 3-5 per hundred 
disintegrations, the ratio of the number of L-conversion electrons to the number of 
unconverted y-rays is 13-5. 

An indirect method of finding the number of conversion electrons depends 
upon the determination of (i) the number of L x-rays produced as a result of 
internal conversion in the L shells, (11) the fluorescence yield of the L shells and 
(iii) the fraction of the internal conversion electrons which originate from the 
L shells. Stahel’s (1935) result for (1) is 25-1 per hundred disintegrations, and 
Kinsey (1948) has obtained 0-47 for the fluorescence yield. These figures, with 
the present result of 79°% for (iii) give 67-5 conversion electrons per hundred 
disintegrations. ‘This method eliminates the fitting of our incomplete RaE 
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spectrum to Neary’s, but ignores any ionization of the L shells which may be 
produced by conversion of the recently discovered low energy y-rays. 

Since there is still some uncertainty regarding the relative intensities of the 
conversion lines, Table 2 has been compiled giving the results up to date. 


$4, THE VERY LOW_ENERGY REGION OF THE SPECTRUM 

Lines 5, c and d of Figure 1 are L-Auger lines associated with the ionization 
of the L shells of the RaE product nucleus by the various y-rays present. Their 
energies are compared with those of the L-Auger lines of ThC, an isotope of 
RaE, in Table 3. The agreement is good. The ThC Auger lines were obtained 
in a previous investigation with this spectrometer (Butt 1950). Line @ has an 
energy of 6-85kev. and can be identified either as an Auger line, as shown in 
Table 3, or as the L internal conversion line of the 23-2 kev. y-ray. ‘This internal 


Table 3. The L-Auger Lines of ThC and RaE 


a b c d 
ONG ae 7:10 8-11 10-39 12-26 12-62 
Hp (gauss cm.) 284-2 303°5 343-5 373°0 378°5 
RaE kev. (6°85) 7:93 10:20 12-70 
Hp (gauss cm.) (279) 300 340 379 


conversion should produce a fairly prominent line, since the unconverted 


23-2 kev. radiation, as seen from Table 1, has an intensity of approximately 1% ; 


we therefore suggest that line a is due to this internal conversion. 

If we assume that the Auger lines just discussed are produced by ionization 
of the L shells as represented by line e only, then, using Kinsey’s value of 0-47 
for the fluorescence yield the Auger lines should have a combined intensity of 
0:53 x 46-7=25 per hundred disintegrations. ‘The combined intensity of all 
four lines a, b, cand dis roughly 24 per hundred disintegrations. ‘This agreement 
is accidental because an appreciable fraction of the area under these lines may well 
be due to the f-radiation of RaD together with any conversion electrons of the 
23-2kev. y-ray. In view of these considerations it seems probable that some 
source self-absorption was taking place at these energies. 


§5. CONCLUSIONS 
It is clear that the disintegration of Ra D via the 46-7 kev. transition can, at 
most, account for about 70% of the disintegrations. Feather (1949) has suggested 
an alternative mode of disintegration involving some of the low energy y-rays 
which might account for the remaining percentage (Figure 3). So far, however, 
there is no experimental evidence to confirm this scheme. 


Parity Spin 


™RaD = OO 
/3 ~ Wkey. 
pRaE > =] 
31-3 kev, 
eRe = 2 
210 
amr = 0 


Figure 3. 
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As regards spin and parity relations, since the even—even nuclei RaD and 
RaF have zero spin and even parity and since the RaE f-spectrum is almost 
certainly second forbidden (Feather and Richardson 1948) we conclude that the 
ground state of the RaE nucleus has spin 2 and even parity (this assumes the 
Gamow-Teller selection rules). Again, assuming the end point of the RaD 
B-spectrum is about 15kev. the transition is an allowed one (Feather and 
Richardson 1948), and the RaE* nucleus has spin 1 and even parity.. From this. 
we conclude that the 46-7 kev. radiation carries away one unit of angular momen- 
tum, and that, since there is no parity change, magnetic dipole radiation will 
probably be of comparable importance with electric quadrupole. The high 
value of 13-5 for the L internal conversion coefficient points to the presence of 
a large amount of magnetic dipole radiation. 

As has been indicated, we have obtained no convincing evidence regarding 
the RaD continuous spectrum. Its end point could easily be masked by the 
presence of the intense Auger lines and its expected rapid rise partially annulled 
by self-absorption in the source. Efforts to obtain a cleaner source by activation 
with radon gas have so far failed. 

It will have been noted that no internal conversion line of any of the five 
recently reported y-rays has been detected, with the possible exception of the 
6:8kev. line. This can be explained on the grounds of the weakness of the 
radiations and the fact that the resolving power of the instrument was too low. 
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ABSTRACT. A proton recoil chamber with a thick polythene layer was used. to deter- 
mine the angular distribtuion of the neutrons from thick targets of various light elements 
bombarded with 8 Mev. deuterons. In all cases a pronounced forward maximum in intensity 
was found, the half width of the distribution decreasing with increasing atomic number. 
The absolute intensity of neutrons in the forward direction decreases as the atomic number 
increases but the dependence on atomic number is not a smooth function. Most elements 
show the presence of multiple forward peaks. In a few elements the effect of reducing the 
incident deuteron energy was investigated. ; 


§1. INTRODUCTION 


OR low bombarding energies the (d, n) reaction has usually been considered 
} to proceed in two successive steps, in accord with the compound nucleus 

theory of nuclear reactions (Bohr 1936). Measurements made on the 
angular distribution of the neutrons emitted from light elements bombarded with 
deuterons of energy approximately 6 Mev. are in agreement with predictions from 
this theory (Pool 1938). 

At higher energies however, because of the low binding energy of the deuteron 
and its consequent extended structure, other types of interaction have been 
predicted (Oppenheimer 1935, Dancoff 1947a). In one possible interaction the 
deuteron grazing the edge of the nucleus, or coming within the range of the specifi- 
cally nuclear forces, may cause the proton to be “stripped off’ and captured, while 
the neutron continues in approximately the same direction as the incident deuteron 
with approximately the same velocity. Calculations have shown that at 200 Mev. 
this stripping process is the most important reaction (Serber 1947, Dancoff 1947 b). 
Experiments with incident deuterons of energy 200 Mev. (Helmholz, McMillan 
and Sewell 1947) and of 15 Mev. (Roberts and Abelson 1947) showed the existence 
of this stripping process, characterized by the pronounced forward maximum in 
intensity of the neutrons. 

The suggestion has been made that even at relatively low energies the major 
part of the (d, n) cross section could be attributed to the stripping process (Peaslee 
1948, Butler 1950a,b). Experiments in the energy region 15-20 Mev. showed 
approximate agreement with the angular distribution predicted by the stripping 
theory, but disagreement with the predicted neutron energy distribution and the 
variation of neutron yield with the mass number of the target element (Ammiraju 
1949, Falk, Creutz and Seitz 1948, 1949). As no measurements have been made 
on the (d, n) reaction in the region of 8 Mev., which is above the potential barrier for 
most light and medium nuclei, and since this may be the transition region for the 
onset of stripping, the present series of measurements was made. 
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§2. APPARATUS 

Deuterons were accelerated by the 37-inch cyclotron in this laboratory, and 
had an initial energy of approximately 8mev. ‘The deuteron beam entered a 
subsidiary target chamber through a 5-foot extension tube connected to the main 
target chamber. Adjustable slits in this extension tube served to limit the angular 
spread of the beam, and achieved some small measure of energy resolution in the 
cyclotron stray field. ‘The main features of this subsidiary target chamber were 
an insulated target box of thin brass, in which the various targets were mounted, 
and an insulated flap of lead which could be rotated through a vacuum-tight joint, 
so that it either intercepted the beam before it struck the target or lay completely 
out of the path of the beam. Since both the target box and the lead flap were 
insulated from earth, it was possible by connecting them to a beam-integrating 
circuit to measure the total charge which either had received during a particular 
experiment. All parts of the target chamber which could be bombarded by stray 
beams were covered with thick lead foil, to reduce the neutron background. 

Most of the experiments cited above were performed with radioactive threshold 
detectors to measure the neutron angular distribution. ‘The present measure- 
ments were made using a proton recoil chamber, with a thick polythene layer. 
Although the energy resolution of the thick layer proton recoil chamber is very 
poor, it has been shown that the efficiency is high, and increases approximately 
linearly with incident neutron energy (Allen and Wilkinson 1948). This form of 
neutron detector will therefore give more weight to the higher energy neutron 
groups. By varying the bias on the recording instrument itis possible to use the 
chamber as a threshold detector with variable threshold (Barschall and Bethe 1947). 
The proton recoils were detected in a parallel plate ionization chamber with a grid 
(Bunemann, Cranshaw and Harvey 1949). Argon ata pressure of 10 atmospheres 
was used as the chamber filling. . 

At this pressure the dimensions of the chamber were such that a proton with 
energy greater than 8 Mev. passed through the grid. A higher pressure would have 
been desirable to contain all the recoils within the region between the polythene 
layer and the grid, but the intense y-ray background from the cyclotron tank 
caused too much piling up of pulses. Although it was not possible to run the 
chamber with voltage saturation of the pulses, the change in pulse amplitude for 
any reasonable voltage variation of the regulated high voltage supply, due to mains 
variation, was negligible. An amplifier with a rise time of 4 wsec. was used, as the 
chamber operated under electron collection conditions. The differentiating time 
constant was made approximately 10 ysec., in order to minimize the piling up effect 
of y-ray and neutron recoils formed in the high-pressure argon. Thepre-amplifier 
was mounted directly on the recoil chamber, and the pulses fed back through a 
cable to the main amplifier and recording gear, located remotely. The pulses 
were recorded by a 20-channel pulse analyser. The use of a pulse analyser enabled 
the data obtained to be treated so that the effect of varying the threshold energy 
could be investigated. ‘The recoil chamber was mounted on a light framework, 
which moved around the target on an arc of a circle of one metre radius, in sucha 
way that the centre of the chamber was at the same level as the beam. The 
experimental arrangement was such that the neutrons could be detected at angles 
up to 75° with respect tothe cyclotronbeam. The monitoring beam integrator and 
the analyser recording counters could be simultaneously switched on. 'To keep 
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down the background of stray neutrons from the cyclotron tank itself, a thick wall 
of parafhn was used. All the targets were ‘thick’, i.e. thicker than the deuteron. 
range in the target material. 


§3. PROCEDURE 


The current to the target was approximately 0-1 ja. throughout the experi- 
ments, and at each angle the neutron counts were recorded for two runs each of 
33 u-coulombs. Since, as was to be expected theoretically and was actually found 
in preliminary experiments, the neutron angular distribution was symmetrical 
about the direction of the incident deuteron beam, measurements were made on 
one side of this direction only. Despite the paraffin shielding there was still a 
considerable background of stray neutrons from the cyclotron tank, and so the 
following procedure was adopted to enable corrections to be made for this effect. 
Measurements had shown that the neutron yield from lead at 8 Mev. is very small, 
probably the major contribution being due to light element impurities, and an 
oxygen surface film. At each angle therefore the neutron yield was measured both 
from the target and from the lead flap, with equal irradiations. The whole of the 
neutron yield when the lead flap was being irradiated was taken to be a measure 
of the background effect, and was subtracted from the yield obtained when the 
target proper was irradiated. ‘This background correction varied from approxi- 
mately 24°% of the yield in the forward direction in the best cases (Li, Be, B) to 
approximately 15% in the worst case (O,). Of this correction only about 30% 
was attributable directly to the yield from the lead flap. The energy of the 
deuteron beam was varied by putting lead foils of known thickness in front of the 
lead flap, so that the beam energy was lowered both when measuring the yield from 
the target, and while measuring the background effect. 


§4. RESULTS 

The measurements of the angular distribution are presented graphically in 
Figures 1 to 12. The relative intensities of recoil protons, and consequently of 
the emitted neutrons, are plotted as a function of the angle of emission of the 
neutrons with respect to the incident deuteron beam direction, the laboratory 
coordinate system. All of the distribution curves are characterized by the narrow 
forward peak of intensity, similar to that found in the measurements at higher 
energies. Since the amplitudes of the recoils were recorded ona 20-channel pulse 
it was possible to count all proton recoils with energies greater than a certain ‘ bias 
energy’, obtainable from a calibration of the detecting equipment with «-particles 
of known energy. 

All the measurements, except those of Figures 8 and 10, were taken with a bias 
corresponding to 1-5 Mev., this being the lowest bias which could be used without 
y-ray pile-up making a significant contribution to the counting rate. ‘The measure- 
ments of Figures 8 and 10 were recorded at a bias energy of 2Mev. Figures 1, 2, 
8 and 10 show the presence of subsidiary maxima, such as have been previously 
reported (Allen, Nechaj, Sun and Jennings 1949). From Figures 7, 8, 9 and 10 
it is seen that in magnesium and aluminium the subsidiary maxima appear only at 
the higher bias. Almost all these elements which do not show subsidiary maxima 
at the lowest bias show the feature at high bias values. In certain elements 
(Figures 2, 8 and 10) the effect of lowering the beam energy by about 1 Mev. is 
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shown. It can be seen that, although the yield is considerably reduced by this 
decrease in incident deuteron energy, the general form of the angular distribution 
function is unchanged. Interpretation of the results in any detail is difficult 
compared with recent (d, p) experiments (Holt and Young 1950, Rotblat 1950, 
Burrows, Gibson and Rotblat 1950), since neutron groups arising from different 
states are not resolved, and so the angular distribution measured is the cumulative 
effect of the distributions from many states. ‘The presence of the forward peak 
however tends to add weight to the suggestions that even at low energies the major 
part of the (d, n) cross section can be attributed to a stripping process. 
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Deuteron Bombardment of Neon 


By R. MIDDLETON anp.C. T. TAI 
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Communicated by H. W. B. Skinner; MS. received 22nd March 1951 


ABSTRACT. Neon gas of normal isotopic constitution was bombarded by 7:8 mev. 
deuterons from the Liverpool University cyclotron. Photographic plates were used as 
detectors, enabling the protons from the reactions ?°Ne(d, p) and ??Ne(d, p) and alpha- 


particles from *°Ne(d, ~) and ?*Ne(d, x) to be studied, as well as a group of inelastically 
scattercd deuterons. 


SF UN TRO DCCA TON 
HE reactions *°Ne(d,p) and ?*Ne(d,p) have recently been studied by 
| several authors (Elder, Motz and Davison 1947, Zucker and Watson 1950) 
but all in cases low bombarding energies have been used, and few proton 
groups have been observed. The reactions 2° Ne(d,«) and 22Ne(d,«) have not 
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been previously investigated and are capable of giving information about the low- 
lying states of $F. By using the gas scattering technique of Chadwick, May, 
Pickavance and Powell (1944) it is possible to study the spectra and angular 
distributions of the charged particles emitted by 7°Ne and ??Ne under deuteron 
bombardment. This paper describes measurements of the energy of the particles 
emitted by 2°Ne and 22Ne when bombarded with 8 Mev. deuterons. ‘The recent 
theory of Butler (1950) for stripping reactions makes the angular distributions 
for (d, p) reactions of use in determining the spins and parities of the states involved 
which will be discussed in a later paper. 


§2. EXPERIMENTAL. TECHNIQUE 


The experimental method used has recently been described by Rotblat et al. 
(to be published) who give details of the scattering geometry. _Spectroscopically 
pure neon of normal isotopic constitution was used at pressures between 10 and 
20cm.ofmercury. Ilford C2 and E1 emulsions 100 » thick were used as detectors, 
the plates being so located as to enable measurements to be made on the scattered 
particles over an angular range of 15° to 165° with respect to the deuteron beam, 
and at distances of 5 to 20cm. from the scattering volume. Exposures of about 
1 micro-coulomb of deuterons were found to give a suitable track density. 

Observations of proton and alpha-particle tracks were made at various angles, 
using standard microscope technique. At each angle the horizontal projections 
of all the tracks were measured within an area of about 10mm? containing some 
3,000 tracks. A collective correction to determine the true range was applied to 
the track lengths at each angle, corresponding to the angle of incidence of the 
particles in the emulsion. ‘The particles were identified by general appearances 
(e.g. grain density and end-scattering) and by comparison of the measurements 
with the calculated O-values corresponding to the different angles. 


§3. RESULTS 
(1) Proton Groups 


Range measurements have been made for ten angles, in the region from 30° to 
150°. Figure 1 shows a typical histogram for the angle 110°. The mean ranges 
of particle groups and their half widths were deduced from the histogram by 
fitting gaussian curves to the individual groups. In the case of a weak group the 
results of several different angles were superposed to obtain the Q-value corre- 
sponding to that group. By comparing two such superposed histograms at 
fairly widely separated angles it was possible in most cases to detect an incorrect 
assignment of the neon isotope responsible for the group of particles. 

The energies of the groups were obtained by using the range-energy relation- 
ship of Rotblat (to be published). A correction was applied to the range of a 
group according to the distance travelled by the particles through neon gas. 
For this purpose we assumed the value of the atomic stopping power of neon to be 
1-250 for protons and deuterons, and 1-235 for alpha-particles (Livingston and 
Bethe 1937). Since the correction for most of the proton and deuteron groups 
is less than 0-1 Mev., any uncertainty in the value of the atomic stopping power 
has a negligible effect on the calculated Q-values. 

The energy of the incident deuteron beam at the centre of the scattering volume 
was determined by measurement of the mean range of the elastically scattered 
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deuterons at various angles. The contribution to the elastically scattered peaks 
due to the presence of 10° of 22Ne was subtracted on the assumption that the 
isotopes have equal differential scattering cross sections. At certain angles, see 
Figure 1, there is a proton group merged in the deuteron peak which was either 
subtracted or, if sufficiently small compared with the elastic peak, neglected. 
Since the position of the centre of the scattering volume changes slightly with the 
angle, there is also a corresponding change in the incident deuteron energy due to a 
slightly different path length through the neon gas. In Table 1 are listed the 
measured values of the incident deuteron energy as a function of angle, and also 
the values corrected to the 90° position. From the average of the latter values were 
deduced the deuteron energies at various angles which were later used for the 
calculation of Q-values. 


D ineiastic 


t 


D elastic 


200 300 400 


Figure 1. Histogram of proton and deuteron groups at 110°. 
(1 scale division=1-7,) 


With statistically well-established groups, it is possible to assign a particular 
proton group to *°Ne or **Ne by comparing the Q-values calculated from the range 
at various angles. If the reaction has been assigned to the incorrect isotope a 
systematic deviation occurs in the Q-values of more than 0-1 to 0-15 Mev. in the 
complete angle range. In the case of such groups, this deviation is quite apparent, 
and immediate differentiation can be made between the isotopes possibly respon- 
sible for the reaction. ‘The reaction responsible for groups statistically not well 
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established, having a probable error of the order of + 0-07 Mev., was determined 
by a method of superposition of histograms. This method of superposition has 
yielded fairly definite assignment of most groups, only one or two remaining 
uncertain. 


Table 1. Incident Deuteron Energy as deduced from the Range 
of the elastically scattered Deuterons 
Angle 30° 374° 45° 524° 60° 
Ey (mev.) observed Dass Weds —— ASO esis He 7D) 
ae HOS) > aeOH0S, | seW0s se O08 ae Mei) 
E> (Mev.) corrected to scattering position 90° 7:800 7-785 7-823 7:809 7-786 
Ey (Mev.) used in calculating Q-values ceases Fees “te 7fOSHONS 
10 Ol =O O15 ete: Otte 0) Ol eect Oni 
Angle 824° 100° 110° iIBS* 150° 
Ey (Mev.) observed 72834 7-800)» 7-190) 7-795 eee oS 
= OF03" =E10-037 9 st 0-03) a0: OS Ons 
Ey (Mev.) corrected to scattering position 90° 7-835 eS 7:798 7°802 7-809 
E> (mMev.) used in calculating QO-values 7-800 7°800 (2/97 OS 7°785 
S(O aeRO ar Oil sasOeil se Wer 
Table 2 
(a) *°Ne (d, p) #4Ne 
Proton group No. 12, 125 Ba 12, |Pe 1B, 
QO-values (Mev.) 4°54 4-21 2°86 ilo7S 0-81 —0-17 
+ 0-04 + 0:03 + 0-04 + 0-02 +0:04 +0-02 
Energy level of 24Ne 0) 0-33 1:68 PIS) 3°73 4-71 
+ 0:05 + 0-07 +0-05 +0:06 +0-05 
Level (Z & W) 0 0-34 1-78 2°84 3-6 — 
Level (F & T) 0 0-38 1275 2:76 — = 
+ 0-02 + 0:04 +0-05 
Proton group No. RP, 125 Bo Pio ea Bae 
O-values (Mev.) —0-90 —1:20 (—2-12 --2-76 —3:74  —4-37 
+ 0-02 +0-02 +0:08) +0:04 +0:04 +0-04 
Energy level of ?4Ne 5-44 5-74 (6:66 7-30 8-28 8-91 
OS OOS + 0-09) + 0-06 +0:06 +0:06 
Level (Z & W) — oo _ — — — 
Level (F & T) — — — — — — 
(0) **Ne (d, p) Ne 
Proton group No. 1 2 
O-value (mev.) 1:17 (—0-04 
+0-05 + 0-08) 
Q-value (Z & W) ENS — 


The values in parentheses are uncertain as to the responsible isotope. 


with the more likely isotope. 


They are listed 


(Z & W) Zucker and Watson 1950; (F & T) French and Thomson 1951. 


Som 
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In Table 2 are listed the results obtained for the various proton groups 
with their probable errors, together with the recent results of other workers. 
Errors, other than statistical errors, may be attributed to: (i) an uncertainty in the 
incident deuteron energy of +0-01 Mev., (ii) an error in the range-energy relation 
of +0-02 mev. (since our Q-values are averaged over different angles, different 
portions of the range energy relation are used, thus minimizing the error), (ili) an 
uncertainty in the position of the photographic plates, which will affect the O-values 
by about +0-002Mev. In the most unfavourable case, the sum of these errors 
would introduce an error of + 0-025 mev. in the O-values of Table 2. 


(11) Inelastically Scattered Deuterons 


In addition to the observed proton groups an inelastic deuteron group due to 
20Ne has been observed. The deuterons were identified by the large systematic 
error occurring in the O-values when assigned to a (d, p) reaction. ‘The inelastic- 
ally scattered deuterons correspond with an energy level of 1-66 + 0-02 Mev. in ?°Ne. 
This agrees with the level of 1-64+0-05 Mev. from observation of y-rays by 
Jelley (1950) and is probably the same level (1-4 Mev.) as was observed from the 
inelastic scattering of protons by ?°Ne (Powell et al. 1940). 


(111) Alpha- Particle Groups 
The alpha-particle tracks were distinguished from proton and deuteron tracks 
by virtue of their larger grain density and smaller scattering. A further check 
was given by the requirement of consistency in the resulting O-values at various 
angles. ‘The ranges of the tracks were measured to within }y by using a 96 x 
objective in conjuction with a 10x or 15x eyepiece. Figure 2 illustrates two 


a, 
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Figure 2 (a). Histogram of alpha-particles at 35°. 
1 scale division=0°83 ,) 


Figure 2 (b). Histogram of alpha-particles at 1323°. 
(1 scale division=1 ,) 


typical histograms (35° and 1323°). ‘The alpha-particle tracks were measured for 
ten angles and a similar analysis applied to the histograms as for the case of the 
protons. In the case of the alpha-particles the assignment of the responsible neon 
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isotope is easier, since the corresponding difference in the Q-values for different 
angles is larger. The contribution of ?2Ne was found to be very weak, and at 
many angles unobservable, as was also the case with the protons. 

The correction due to the stopping power of neon is larger than that for protons 
and deuterons. An intrinsic check on the adopted value of the atomic stopping 
power of neon, namely 1-235, was obtained by comparing the calculated energies 
for alpha-particles belonging to the same group, but having traversed different 
path lengths in the neon. The Q-values obtained are listed in Table 3. None of 
the groups listed from the ?°Ne(d,«) reactions has previously been observed. 


Table 3 
{a) ?°Ne (d, «) 18F 
Alpha-particle group No. Ly dy on Ola Os 
Q-value (Mev.) 2:78 hoy} 0-95 0:58 0-17 
+ 0-02 +0-02 O02 + 0-05 + 0-04 
Energy level 18F Ground 1:05 1-83 2°20 Deon 
+ 0-03 +0-05 +0:05 +0:05 
Alpha-particle group No. Og * Oly Xs hg X49 
Q-value (Mev.) —0-45 —1-14 — 1-64 —2:23 — 2°83 
+ 0-07 + 0-03 +0:10 +0-08 +0-10 
Energy level 18F B23 3°92 4-42 5:01 5-61 
+ 0:07 + 0-04 +0-10 +0:08 + 0:10 


* Group may be a doublet. 
(b) ??Ne (d, a) 2°F 


Alpha-particle group No. 1 2 
Q-value (Mev.) 2:62 2:05 
+0-10 + 0:07 
Energy level 2°F Ground 0-57 
+0-13 


The two groups from 7*Ne (d, a) ?°F correspond to the first two proton groups observed 
by Allen and Rall (1950) from the reaction !°F (d, p) 2°F. 


By assuming that the observed alpha-particles of longest range correspond to 
the transition to the ground state of 18F, and utilizing the known masses of ?°Ne 
(Mattauch and Flammersfeld 1949), the mass of 18F was calculated to be 
18-00672 +0:00006. ‘This is to be compared sae the value 18-0065 + 0-0006 
given by Bethe (1947). 

Finally it is interesting to observe that the level spacing of 18F bears a close 
resemblance to that of 17O (Rotblat 1950). From our determinations it would 
appear that the level spacings of !8F above 2-60 mev. have re the same value, as 
for 1’O above 3-88 Mev., namely 0-60 mev. 
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Further Results on the Fluctuation Problem in Electron-Photon 
Cascade Shower Theory and the Probability Distribution 
Function 


By H. MESSEL 
School of Theoretical Physics, Dublin Institute for Advanced Studies 


Communicated by W. Heitler; MS. received 1st May 1951 


ABSTRACT. Results are presented both for the fluctuation and distribution of electrons 
in electron—photon cascade shower theory. ‘The results were obtained in approximation A 
and are for the case of a single primary electron giving rise to the shower. 


§1. INTRODUCTION 
' HE so-called fluctuation problem in electron—photon cascade shower 
| theory may be roughly divided into two parts: 

(a) If we consider an electron (or photon) of energy E, impinging on a 
layer of material of thickness t, what is the probability that exactly N electrons 
(or photons) with energies greater than E=«E, will emerge? 

(b) If we know the average number of particles above an energy E=c«£, at 
any given depth ¢, due to a primary electron or photon of energy Ey, what © 
fluctuations from this average number may we expect? 

In order to answer (a) knowledge of the actual distribution function for 
electron—photon cascades is required. Until recently it was also believed that 
(b) could not be answered without solving (a) first (Arley 1949). ‘This, however, 
has not proved to be the case. 

Janossy (1950) has shown that a cascade could be represented by the following 
set of equations, when ionization loss by electrons is neglected. 


(%) 1 
LINED = | £E%/(1—€1), 1EM(ey 1) — GM, Nw) de, i= 1,2 
0 


Past oy (1) 
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where the index 7 refers to the nature of the primary: 7=1, electron primary, 
i=2, photon primary, and G™ is the generating function 


GE, Oy Ut) =D) OO le myers 1) eeas O (2) 


Ny, No=0 
with ®(e, ,, m,t) the probability that at a depth ¢ a primary of energy E, has 
given rise to m, electrons and nm, photons with energies greater than H=«E, and 
an arbitrary number of electrons and photons with energies less than FE. Note 
that in equation (1) we have omitted the variables U, and U, for convenience only. 

In principle equation (1) should be capable of being evaluated by numerical 
methods, since the initial conditions are known, and thus, the answer to (a) 
obtained. However, this has proved to be a very difficult task and to date no 
satisfactory solution has been obtained. 

Bhabha (1950), Bhabha and Ramakrishnan (1950) and Scott (private 
communication*) have also given the solution of the fluctuation problem a good 
deal of attention. ‘They set up the fundamental equations for the distribution 
function in a manner different from that of Janossy, but made no essential progress 
in their solution. 

Recently Janossy and Messel (1950, hereafter referred to as I) Goole a step 
forward. ‘They gave analytical solutions and fairly extensive numerical results 
for the fluctuations in number about the average. The correlation between 
electrons and photons in a shower has also been worked out, and numerical 
results have been given. ‘That work has confirmed quantitatively what had long 
been suspected, viz. that on either side of the cascade maximum the fluctuations 
deviate wildly from those which would be expected if the distribution in number 
were Poissonian. It should be pointed out that the solutions obtained in I were 
in approximation A, 1.e. when ionization loss by electrons was neglected. ‘The 
computations even in this approximation were extremely heavy and one wonders 
what hope there is of obtaining numerical results when ionization loss is included 
(approximation B). In fact the obtaining of an analytical solution in approxi- 
mation B for the second moments of the distribution is most forbidding—especially 
when one considers that we have not as yet obtained a completely satisfactory 
method for obtaining values for the average numbers in approximation B! 

The results as presented in I were of value for interpreting the physical 
behaviour of the fluctuations, but of little direct use to the experimentalist. 
We shall present in a 'T'able below, numerical values of the relative fluctuation, 
which may be utilized more directly. : 

In view of the fact that there is little hope of immediately obtaining a solution 
for the electron—photon distribution function, we have decided to make use of 
an assumption first introduced by Arley (1949). He assumed that the distribution 
function could be well approximated by the Polya distribution (discussed below), 
which contains, essentially, two parameters, the first and second moments. 
Since we now have the values of the second moments it is directly possible, using 
the above assumption, to evaluate the distribution function. This is done 
below. (Note that Arley found values for the distribution function also; 
however, 1n his case he did not have the values of the second moments and had to 
approximate them.) We shall only present results for the case of a single 
primary electron producing electrons, since this is the case for which the second 
moments have heen evaluated. 


* The author is indebted to Dr. W. T. Scott for communicating his results prior to publication. 
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S25 DHE POL YyA ‘DISTRIBUTION 
The Polya distribution function expressing the probability of finding exactly 
NN particles with energies greater than <Hy at a depth ¢ is given by 


~ 1.0. +8)...1+(N=1)d} (_ N_\* 
poe Seen a) P(e, 0, t) 
(20, NeOand V=1,2 0.) 5 6 ane (3) 
with PU) = (le ON eee eae ade (3 a) 


In our case N expresses the average number of electrons with energies greater 
than <H, at a depth ¢ cascade units due to a single incident primary electron of 
energy E). ‘The parameter b is expressed by 


fa eeeNC ONY Nee ae Ae (3b) 
One obtains the Poisson distribution function by simply allowing the parameter 
4 to approach zero, thus oe A 
P(e, 0, t)>exp(—N), Pe, N 1) ae eee (4) 
6-0 : 


Note that for a Poisson distribution ie Plate fluctuation, defined as 
5 =(N?—N’)"/N is simply equal to Nos “The significance of the parameter b 
appearing in the Polya distribution is now obvious. It expresses the difference 
between the square of the relative fluctuation for the actual distribution and that 
for a Poissonian one. 

If one knows both N and 6 then the distribution function P(e, N,t), given 
by (3), may be immediately evaluated. ‘The value of N may be obtained from 
the extensive tables given by Janossy and Messel (1951) and the value of 6 found 
using the results given in I. 


63. RESULES 
In Table 1 we have tabulated N, b, —1/b,1 +5N and the relative fluctuation 
p= (N?— N'Y" /N as functions of the depth ¢ measured in cascade units and the 


. energy ratio expressed by —loge. Atacascade maximum and its neighbourhood 
~ the value of the parameter 5 takes on negative values. In these cases we have 


omitted to give the value of 6. It may be pointed out here that the fact that b 
takes on negative values proves by itself that the true distribution function is 
not exactly a Polya distribution, for if it were, then in all cases 6 would necessarily 
be greater than or equal to zero. 

In Table 2 we have evaluated the Polya distribution function using Table 1 
for the case —loge=2-0. The values for a Poisson distribution are also given 
for comparison. For the cases in which )<0 we could use only the Poisson 
distribution function. ‘This, however, should not bring in serious errors for, as 
shown in I, the fluctuations at the cascade maximum where )<0 are pearly 
Poissonian. 

The value of P(e, N,t) for —loge =3-0, 4-0 and 6-0 and any N and ¢ may 
also easily be obtained by using ‘Table 1 and formula (3). 

Figure 1 gives a plot of ithe relative fluctuation 5 against the depth ¢ for 
—loge=2-0, 3-0, 4-0 and 6:0. We have also plotted the relative fluctuation 
assuming the distribution to be Poissonian, for comparison. Figure 2 gives the 
probability distribution curve for both a Polya and a Poisson distribution, for 
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Table 2. The Probability Distribution Function —log P(e, N, t) 


ving the negative logarithm of the probability of a shower containing exactly N electrons with energies greater 
1/100 of the primary energy, at a depth t cascade units. The shower was initiated by a single primary 
sctron. We have also given the results for a Poisson distribution, for comparison. (The top line of each group: 
is the Polya and the lower line the Poisson distribution.) 


t 2 3 4 5 6 i 8 9 10 12 14 16 18 20 


, . ner 1:513 0-890 0569 0-345 0-126 0-0432 0-0139 0-0043 0-0012 
0 2545 3-249 3-227 2-734 2-083 4.478 0.995 0-644 0-405 0-149 0-0517 0-0169 0-0052 0-017” 


; ang 1:017 0620 0-504 0-509 0-734 1-121 1:590 2-106 2-601 
1 1-777 2-375 2-356 1-935 1-402 5.946 0-635 0-473 0-435 0-613 0-976 1-426 1-926 2-400 


Pees 5 ee 108120 0-609, 0°66 5.20.832..1°31 te 1-841 3550 Roe 8907 
2 1310 1-802 1-786 1-437 1-022 9.745 0.576 0-603 0-766 1-377 2-202 3-136-4-148 5-099 


Men idared Gop aie) o.n10.0:772.0°130, 0-94011-228°1-879.2-487 3-009 13-5322 3166s 
ee 0208140581592 1-115 0-819 9661 0-693 0-909 1-27412-318 3-603 5-022 6-546 70974 


ee man 0846 0-953 1-284 1-654 2-441 3-100 3-619 4135 4-420 
4 0854 1-134 1-123 0-918 0-740 9.537 9.935 1-340 1-906 3-383 5-130 7-033 9-069 10-974 


. neg 1:006 1-235 1-676 2-109 3-001 3-694 4-204 4-711 4-923 
> 0°785 0-958 0-951 0-818 0-758 9.298 4.974 1-868 2-635 4-546 6-753 9-141 11-688 14-071 


é ; ; gee 1236 1-565 2:103 2-582 3-560 4-277 4-774 5-270 5-409 
6 0-795 0-863 0-858 0-797 0-855 4.444 1.692 2-475 3-443 5-787 8-456 11-328 14-387 17-247 


ye oar tentg 1°523 1-934 2558 3-068 4-116 4-851 5-338 5-818 5-882 
7 0872 0-834 0-832 0-843 1°019 4.457 9.177 3.149 4-318 7-096 10-225 13-582 17-153 20-490 


ng, 1:858 -2:335 3-034 3-564 4-672 5-418 5-888 
8 1-007 0:863 0-864 0-947 1-241 4 858 5.790 3-880 5-252 8-46212-053 15-894 


‘erg 2235 2°762 3°526 4-068 5-227 5-982 
9 1-193 0-943 0-948 1-103 1:514 5.555 3.315 4-664 6-236 9-880 13-932 


: : ; ; : 2-648 3-211 4:034 4:578 5-781 
10 1-425 1-069 1-077 1-304 1-833 2-718 3-954 5-492 7-26611°344 


; ; : : : 3-094 3-679 4-554 5-094 
11 1:699 1-236 1-247 1-546 2-194 3-228 4-636 6-363 8-338 


: : ; ; : 3:568 4:164 5:084 5-614 
2 2010 1-442 1-455 1-826 2°592 3-775 5-355 7-271 9-447 


ons 4:069 4:663 5-623 
Seamer? 65001-6810 17698 2-441, 3-025 oe tng og o1g 


; 499 47593. 5-175 
iA 27734 1-954 1-973 248653400. 7'0o oa 


Peete den: rong acoge 8 13925-6098 
BS 2243142 2-256 2-278 2-865 3-985 252 2044 


: F 508 2°205 
16 3:578 2:586 2-612 3-271 4-508 6-287 
17 4-041 2-942 2-971 3°702 5-058 
18 4528 3-323 3-355 4-158 5-632 


19 5°039 3:728 3°763 4-638 6-230 


20 5°572 4:155 4-193 5-140 6-850 
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the case of —loge =2-0 and t=8 cascade lengths. Finally in Figure 3 we have 
plotted the logarithm of the ratio of the probability of finding N particles with 
energies greater than «EZ, at a depth ¢ for a Polya and a Poisson distribution. 
Curves are given for t=7, 8, 10 and 14 cascade lengths. 


§4. DISCUSSION OF RESULTS 
The outstanding feature of the curves in Figure 1 is the extreme rapidity 
with which the relative fluctuation 5 increases once the cascade maximum has 
been passed. We note in comparison the rather slow increase (as N—"?) of 8 for 


1-4 -Loge 
2-0 
4 2-0 
Ss 
vo 
Sy 
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— ee YP 
Figure 1. A plot of the relative fluctuation 6=(N?—N ) /N for the case of a single primary 
electron of energy Ey) producing electrons at a depth ¢ cascade units with energies greater 
than H=cE,. (Log (86+1) is plotted for convenience.) The full curves refer to the true 
relative fluctuation, the dotted curves to the relative fluctuations if the distribution were 
Poissonian. The minimum of each curve corresponds to the cascade maximum. 
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Figure 2. The probability distribution function Figure 3. The logarithm of the ratio of 
P(e, N, t) for the case of —log «=2-0 and t=8 finding N particles with energies greater 
cascade units. For comparison the curve for than #= cE, at a depth ¢ cascade units, 
the corresponding Poisson distribution is also given. for our distribution function (Polya) 


and a Poisson distribution. In all cases 
—log «=2.-0. 
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a Poisson distribution. The above behaviour brings out forcibly the fact that 
the probability of finding two particles at a great depth is of the same order of 
magnitude as the probability of finding one particle. This behaviour had 
already been pointed out by Janossy and Messel. 

Both Figures 2 and 3 and Table 2 bring out the above feature more directly : 
that is, the probability distribution function for electron—photon cascades beyond 
the cascade maximum is heavily skewed in comparison with the Poisson 
distribution. For example, the probability of finding two electrons with energies 
greater than one hundredth of the primary electron energy at a depth of 18 cascade 
units is approximately 1-3 in 1,000 whereas for a Poisson distribution the 
probability is of the order of 7 in 100,000. ‘The average number of electrons 
expected under these conditions is only 0-012. For larger values of N and greater 
t the difference is further accentuated. 

The amount of confidence one can have, from a quantitative point of view, 
in our results depends solely on how well the Polya distribution function approxi- 
mates to the true electron—photon distribution function. We recall the well- 
known theorem in statistics that in principle, if all the moments of the distribution 
are known, then the distribution function itself is fixed. Note that the Poisson 
distribution utilizes the first moments only. The Polya distribution contains 
both the first and second moments. From a qualitative standpoint our results 
do throw a good deal of light on the behaviour of electron showers and should be 
of value for the interpretation of these showers, observed in cloud chambers and 
photographic plate work. 

The values of the relative fluctuation 6 may He treated with confidence in 
approximation A. 
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ABSTRACT. A calculation has been made of the ferromagnetic properties of a system of 
electrons using plane wave functions and Bloch exchange integrals. Our results show that 
below a critical temperature T, the spontaneous magnetization decreases slowly with in- 
creasing temperature. At 7, the magnetization drops suddenly to zero as if a change of 
phase were taking place and there is therefore a latent heat of demagnetization at this 
temperature. Above T, the model is paramagnetic. ‘The reasons for the differences 
between the results and those obtained by experiment and on Stoner’s theory of collective 
electron ferromagnetism are discussed, and the possibility that they may be due to the 
neglect of positional correlations among the electrons is considered. 


§1. INTRODUCTION 

EFORE Stoner (1938, 1939a) introduced his theory of collective electron 
ferromagnetism there was no general scheme sufficiently flexible to account 
for the behaviour of the ferromagnetic metals and their alloys in any but 
qualitative terms. For example, the Heisenberg (1928) theory, although it 
explained the existence of ferromagnetism and the existence of the Curie point, 
could not (without ad hoc assumptions) be made to account for the non-integral 
numbers of spins per atom which the observed values of saturation moment at 
absolute zero required. Nor could it explain the departures of the susceptibility 
above the Curie point from the Curie-Weiss law. Stoner, using Fermi—Dirac 
statistics, calculated the properties of an assembly of electrons in a parabolic 
band when their exchange interactions could be represented by aterm of the form 
— Nk6'C?/2 in the total free energy; ¢ is the relative magnetization of the electron 
assembly, (0 <¢ <1), N the total number of electrons, k is Boltzmann’s constant 
and 9’ is a constant which is supposed to be independent of temperature. With 
the further assumption that the holes at the top of a nearly filled band were 
magnetically and thermally equivalent to electrons in a parabolic band, this 
theory was able to explain the above facts which found no explanation on the 
Heisenberg theory. ‘The scheme has recently been successfully developed by 
Wohlfarth (1949 a, b) in a study of alloys of the ferromagnetic metals. But although 
successful in its application the basis of the theory is not above discussion (Van 
Vleck 1945). In particular, there is at the moment no rigorous theoretical 
foundation for the representation of the total exchange energy by an expression 
of the form used by Stoner; this can only be regarded as justified pragmatically 

by the success of the theory. 

In the present paper we shall adopt a different approach; we shall calculate 
the ferromagnetic properties of an assembly of free electrons systematically using 
the exchange energies given by the quantum mechanics of this system. Our 
results show that below a critical temperature 7, the magnetization ¢ decreases 
slowly from unity with increasing temperature. At T,, € drops suddenly to zero 
as if a change of phase were taking place; there is therefore a latent heat of 
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demagnetization at T=T,. Real ferromagnetics do not behave in this way, the 
magnetization falling steadily to zero at the Curie temperature. (This is not 
true for some alloys of manganese with the Group V elements phosphorus, 
arsenic and bismuth, which do behave in the same way as our hypothetical 
model; it is likely, however, that in this case the behaviour is due to other effects.) 
Our results are thus very different from those of Stoner’s theory. However, 
using the same model as we have done, Wohlfarth (1949c) has attempted to 
justify Stoner’s assumption that the total exchange energy is — NkR9@’C?/2 but 
his considerations were confined to the behaviour at absolute zero and thus 
failed to disclose a strong temperature dependence of the exchange energy. 

Now this free electron model is relevant to a discussion of the ferromagnetism 
in iron, cobalt and nickel because of the equivalence to it of the holes in the nearly 
filled 3d bands of these metals. It is not directly applicable in the way it would 
be if the magnetically effective electrons were those in the almost empty 4s bands. 
However, the model is only ferromagnetic under certain conditions (discussed 
in §3) and when these are not satisfied we should expect to be able to apply it 
to the nearly free valence electrons in, say, sodium. (Calculations by Von der 
Lage and Bethe (1947) have shown how closely the wave functions and kinetic 
energies of the 3s electrons in sodium approximate to those for perfectly free 
electrons.) ‘This has been done by Wohlfarth (1950) who has shown that the 
free electron model only gives a specific heat of the same order as the observed 
electronic specific heat provided we make a correction to the exchange energy 
for the correlation effect: otherwise it is very much less than is observed. It is 
possible in our case also that inclusion on the correlation correction would give 
better agreement with experiment and, in particular, make the change from the 
ferromagnetic to the paramagnetic state less abrupt. Unfortunately it is extremely 
difficult to estimate the changes which the inclusion of a correction for electron 
correlation will bring about. Further consideration is being given to this 
problem. 


§2. THE ONE-ELECTRON APPROXIMATION AND CALCULATION 
OF THE FREE ENERGY 


It is well known (Seitz 1940) that on the basis of the one-electron approximation 
in the theory of metals, the energy of an assembly of electrons is 


U= aed | bat) det E | ye) ae 


~ 8ar2m 


+E [fly aCey fe loa(t) Bees de 


MAN 


2 
— © | [penne dt dtp. oo. (2.1) 
AOI EE AW 
Here %,,(r) is the space dependent part of the wave function for state n of an 
electron, and V(r) is the potential energy of an electron at r due to the ion cores 
of the lattice. ‘The summations in the first two terms are over all occupied 
electron states n, and in the third term over all pairs of occupied states. These 
three terms are respectively the kinetic energy of the electrons and the Coulomb 
energy of the electrons with the lattice and among themselves. In the fourth 
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term, which is the exchange energy, the summation is over all pairs of states 
occupied by electrons with mutually parallel spins. 

Now for perfectly free electrons the wave functions % are of the form 
= V—1? exp(zk.r), where V is the volume of the metal and k is the momentum 
vector. The integral in the first term in equation (2.1) is then just k?. The 
integrals in the second and third terms are independent of the states for which 
they are taken and thus these two terms in U are independent of any change in 
the occupation of states. That is, we can ignore the Coulomb energy in con- 
sidering changes in magnetization or changes brought about by temperature. 
The exchange integral between an electron in state k and another (with parallel 
spin in state k’ was first evaluated by Bloch (1928) who obtained 


us RD: 
aa 


The absence of an exchange integral between two states of opposite spin 
allows us to consider our assembly of N electrons as two independent sub- 
assemblies, one of , parallel spins, the other of m, anti-parallel spins. The 
relative magnetization ¢ is defined by 


C= Gr nN bile fan (2.3) 


Hence, if we can find the free energy of the parallel spins F, say, and thus also 
of the anti-parallel spins /,, the free energy of the whole assembly is (2) = F, + Fy. 
The equilibrium magnetization is that for which F(Z) has its least value. The 
problem, then, is the calculation of the free energy for an assembly of electrons 
with individual kinetic energies h?k?/8x?m and exchange energies given by (2.2), 
where only one electron can be assigned to each allowed value of k. These 
allowed values are determined by the Born cyclic condition, which requires all 
components of k to be integral multiples of 27/V1°. If we assign to each state k 
a distribution number ¢(k) which can be zero or unity only, then the energy of 
any particular configuration is 


J kee 


See apa ene ES 
U= Ba Mik) 58 2.4) 
“eee Sii(k) =n ea «a eee (2.5) 


The summations are over all allowed k vectors; for a cubical volume V these 
form a simple cubic lattice of points in k-space. In principle we can now obtain 
the free energy F from the partition function Z calculated from (2.4) and the 


equation 
Zag URE yg OR 
(p(k)) 


where the summation is over all sets of ¢(k) consistent with ¢(k)=0, 1 and, the 
equation (2.5). However, it appears that the present resources of statistical 
mechanics are not adequate for evaluating a partition function of this type which 
is equivalent to solving the general order—disorder problem where each atom 
interacts with every other one. We therefore look for an approximate method. 


The average distribution $(k) will be a function of the temperature and of the 
magnitude of k only: we write it as f(k). If the exchange energy in (2.4) were 
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absent f(k) would be just the Fermi—Dirac function. We take as the internal 


te ae jp okatle R*f(k) — 1. y te ARR) 
82m Zen V |k—k’pP’ 


This is not quite accurate since the average of 4(k)d(k’) is not just the product 
of the averages of 4(k) and d(k’). To get the free energy we need now an expression 
for the entropy. As the main approximation we assume that the entropy S 
is the same function of f as it is for the Fermi—Dirac assembly where no exchange 
interactions are present. This is not to assume that the entropy is unaltered 
because f will be a different function of k when we include the exchange 
interactions S=—k»,[flnf+(1—f)ln(1—f)]. This assumption is due to 
Koppe (1947). Hence for the free energy we have 


tne fy) 
22, V |k-kP? 


# i finf+ i —f)In(1—f)]. 


F=E-TS= 


Let us in the usual way replace these sums by integrals in k-space introducing 
the density of states factor k? ae (one electron to each state k), then 


= es in fk) dk — aa if fF) fhe) Fhe In| ae 2| dk, dhs 
i: BD oes i Ti ple feeds lanl Weep be dee 4 a eee (2.7) 
“0 
vith m= aa) ARRAR (2.8) 
nh ra 6 Ce eS oe re le eee. 8 ane: . 

The form of f(k) is found by putting the variation d/=0. In this way we get the 
equation he re bee! 1—f(k) 

ee ag) ee y , ' = 

os ab Rf(R’)\n || dk’ — n= kT In( ay ). 


The parameter 7, which is equal to the chemical potential, is determined by 
(2.8). This is the equation Bardeen (1936) assumed for the distribution function, 
since the second term on the left-hand side is just the exchange energy of an 
electron in state k with all the others in the assembly. However we have little 
hope of solving this equation for f(k) or of evaluating the integrals in (2.7): we 
therefore apply the approximate method used by Koppe (1947). This consists 
in assuming a form for /(k) containing two parameters: one of these is eliminated 
through the condition that the total number of electrons is constant (2.8), the 
other by making F a minimum, a condition necessary for thermodynamic 
consistency. 

We shall use two forms for f(A), one for ‘low’ temperatures when the electron 
gas is degenerate, the other for ‘high’ temperatures. 


(i) Low T aie 


a—l 
poe ock< fra-Bh, = 


Very low temperatures correspond to «>1. 


<k< 


p= Uy Re 


TI!R 


818 A. B. Lidiard 


(ii) High Temperatures. 
f=X- phe Ok =cip.  f=O0 Ikea) 5 eee (2.10) 


It might be thought that (2.9), and (2.10) would give a very crude approximation 
indeed. But it appears that the free energy F is rather insensitive to the form used 
for f. A calculation, in the same manner as below, was made for the electron 
gas without exchange interactions and the results compared with the accurate 
functions tabulated by Stoner (1939b). ‘The differences were no more than a 
few per cent, whilst the internal energy at high temperatures was given accurately 
as 3n,kT/2. Curves of the free energy, internal energy and entropy are given 
in Figure 2, and details of the calculation are given in Appendix A. 


i) “i OP we 


a 06 08 10 12 
k BS kT/e, 
Figure 1. Approximate forms for the Figure 2. Thermodynamic functions for a 
function f(k): (i) @21; (ii) a<1. perfect Fermi—Dirac gas. Curves marked 


a are obtained using the approximate 
forms for f(k) and those marked } from 
the accurate values obtained by Stoner: 
(i) internal energy ; (ii) entropy ; (iii) free 
energy. 


The form of the integrals in (2.7) is such that it is more convenient to work 
with the parameters «/8 and 1/«=x than with « and B themselves. Thus (2.9) 


gives f(k) for x <1 and (2.10) applies to x>1. By (2.8) we get as the equations 
for eliminating «/8 


ENG 3x yg 
m= gag) (1- F+e8- 5), x<il. 


Se oN 1 
~ Dane 3) x” iat 


Using equations (2.9) and (2.10) to evaluate the integrals in (2.7), we find that 
F, has the form BRT 


2 


3 3 
Fy, = 21,1) A(x) — qm te, B(x) —ny XE) peers (2.11) 


5 
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2 2/8 1/3 
where le)= Ee (=) and le,=é (=) Se Ae), Bx), and C(x) care: the 


8m \a a7V 
following functions of x, in which P(«) =(1—3x/2 +x? — x3/4)-l. 
= Pols Sen\6 x 
ee lee Cm 9) OS, SLOW ras ayitay' Nels ae (2.12) 
= 2. eal we AN 
B(x) = = Se * 1 +(1—«)*] In (1 — 4x) —(1—x)* In(1—x)} 
oe (8 — 16x + 14x? — 6x3 + x*) 
a GOO, Oi x 
US (ee eee rae 
+P (5 5 + Gs) in(5*), YES O NA 
41/3 
me (28 — 16 In 2)x-2/8, jo 
C(x) = P(«— x + 1x3), <1, | 
2 5 ZO a ie ee (2.33) 
te tae Te SAS pull 
= 3inx 5 (x 1) in(—*) + 9 ae 3D eee 
The parameter x is determined 2 putting 0F'/dx =0, which gives 
7 ep AaX)— 7 6,5 (4) eR C (x)= 0p" to 2 (2.14) 


where the primes indicate differentiation with respect to x. At least in principle 
we have thus solved the main problem, which was the calculation of the free 
energy of either of the sub-assemblies of electrons. As we have already remarked, 
we may consider them separately because the exchange energy is zero between 
two electrons of opposite spin. ‘There are m, electrons in the sub-assembly of 
parallel spins and m, in the sub-assembly of anti-parallel spins. Then by the 
definition (2.3) of ¢ we have m,=3N(1+¢) and n,=4M(1—¢). Thus it follows 
from equation (2.11) that the free energy of the whole assembly of electrons, 
which is simply the sum of the free energies of the two sub-assemblies, is given by 


Fee = 7g (+ DAG) + (1-0 5*4(e9)} 


= 32 (1+ D(a) +(1- DB (a) 


2 i= (le Cm Can te. (2.15) 
Phere ze) anh orn (Fy. Se aeh (2.16) 


€ 9 is just the Fermi energy, whilst «; measures the strength of the exchange 
interactions. ‘The parameters x, and x, are given as functions of € and RT/eq 
by the equations 

(1+2)°*4'(e,)~ 7 2(1 + Be.) — 5 C=O, 


3(1- 6)2.A’(x) — 720- £)¥8B’ (x2) — 5 Cl) 0. 
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Just as RT/e, gives a convenient measure of temperature for the assembly, so 
<,/€y gives the strength of the exchange interaction. We shall work entirely in 
terms of these reduced quantities £, RT/e) and €;/e9. However, before considering 
the numerical analysis of these equations, let us examine the conditions under 
which the model is ferromagnetic. 


§3. THE EXISTENCE OF FERROMAGNETISM: 


If we reduce the temperature T to zero, then equations (2.17) give for all 
¢, x, =x,=0; A(0) and B(0) are both unity, so that 


Ne ns Ne = 1d +o + 1-9") - sal oe 


where £(£) is simply the internal energy. The equilibrium value of ¢ is clearly 
that which gives E(€) its least value. The curves of [E(¢)—E(0)]/Neo against ¢ 
for particular values of €;/e) (Figure 3) show that equilibrium is given by 


t=1) for fey - 42" + 1)/5, 1 
C=) for “ejeg = 42h d)/5: f 


The model is thus only ferromagnetic if €;/e9>1-81. No intermediate values 
of ¢ are possible as they are in Stoner’s theory. This sudden change in ¢ as 


ee (eres (ee ee el sn (Pe a (EU fe SC ne |e | eee 
“0065 9,027.03 04 05 06 07 08 09 10 OM Or 02 03 04 US GE W080 r0 
¢ 
Figure 3. The internal energy as a function Figure 4. Curves of free energy against 
of € at T=0. Numbers on the curves relative magnetization for e€j/e9=2-0. 
are the values of €;/€o- The numbers on the curves are the values 
of RT/é. 


RT /e, is increased, was noticed by Wohlfarth (1949c). Equation (3.1) corrects 
an error he made in locating the transition point. The lowest value of «;/e) for 
which E(¢) <E(0) for all € is €;/ey=2: in fact d?E/dl? at ¢ =0 is less than, equal to, 
or greater than zero according as ¢,/e) is greater than, equal to, or less than two. 


§4. MAGNETIZATION-TEMPERATURE CURVES 
Temperature effects are in principle deducible from the equations (2.15) and 
(2.17), which by elimination of x, and x, give Fas a function of and T. However, 
owing to the very complicated nature of the functions A, B and C and their 
derivatives, this elimination can only be done numerically. This we have done 
for €;/e9=2 and for RT/«,=0-2, 0-4 and 0-6. These three curves are sufficient 
to demonstrate the magnetic behaviour of the model. ‘The method was to 
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tabulate the functions A, B and C and their derivatives, and then, using equations 
(2.17) to prepare tables of (1+2)-¥* as a function of x for the three chosen 
temperatures. Inverse interpolation in the tables of (1+) 1% against x then 
gave x, and x, corresponding to any particular ¢. ‘These could then be sub- 
stituted in (2.15) to give F(¢)/Ne, for this value of €. 

The results of the calculation are shown as curves of [F'(¢) — F(0)]/Neo against 
Cin Figure 4. In general, we see that dF(¢)/d¢ =0 has three roots (one of which 
is always ¢ =0), so that in considering the variation of equilibrium magnetization 
with temperature, it would not be sufficient merely to use this equation: that is 
why we have obtained complete curves of the free energy. It is also to be seen 
from Figure 4 that the equilibrium value of ¢ decreases slowly from unity (at 
T=0) with increasing RT/e) until a Curie temperature 7’, is reached at which 
F(¢)=F(0): for ¢;/e9=2 we find RT,/e9=0-52. Above this temperature the 
equilibrium value of ¢ is zero, and so at T= T, there is a latent heat of demagnetiza- 
tion equal to E(0)—E(£): for ¢,/e9=2 this is about 0-05 Ney. Rough values of 
the equilibrium magnetization against RT/«) have been obtained by eliminating 
x, and x, from (2.17) and the equation dF(¢)/df =0: these are shown in Figure 5. 


cab 


(ith | i ‘ees alis 
0 01 0:2 0:3 04 05 06 
RT/e, 


Figure 5. Curve of equilibrium magnetization ¢ against temperature for €;/e7=2-0. 


With these results for the special case ¢,/ey=2, and Figure 3, we can see 
qualitatively how the magnetization will vary with temperature for other values 
of this parameter. Firstly, as «,/ey and the corresponding Curie temperature 
decrease, the (approximately) rectangular area enclosed by the (¢, 7) curve, will 
become thinner, the value of ¢ for T=0 being always unity. This will continue 
until ¢;/ey9 =4(2"3 + 1)/5 when the rectangle is vanishingly thin: we have already 
seen that below this value the model is not ferromagnetic. Clearly if we now 
increase e€,/€y above 2, this effect will be reversed, and the area enclosed by the 
magnetization curve will become broader. We should also expect the value of 
¢ at T=T, to decrease, and it is possible that, by increasing the exchange inter- 
action sufficiently, we might get a completely second order transition to the 
paramagnetic state, without the abrupt change in magnetization found here. 
However, we cannot be certain of this without repeating the calculations for other 
values of e;/€9, and in view of the general divergence of our results from experiment, 
it was felt that the labour would not be justified. ‘Thus, even if we suppose that 
for a sufficiently high value of ¢,/e) the model would reproduce the behaviour of, 
say, nickel, dilution of nickel with copper must bring the value down into the 
region we have examined, whereas none of the ferromagnetic nickel-copper 
alloys show any behaviour such as we have found. | 
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§5. DISCUSSION 

The general conclusion to be drawn from our calculations is, therefore, that 
although the model reproduces several of the characteristic features of ferro- 
magnetism, it cannot be used to provide an exact detailed description. ‘This 
arises because it predicts a sudden drop in magnetization at the Curie point 
instead of the steady fall to zero found experimentally. Also the specific heat 
curve of, say, iron (see Awbery and Griffiths 1940) has a fairly broad (finite) 
peak at the Curie point whereas our model yields a latent heat, which would be 
observed as an extremely narrow tall peak. We can conveniently discuss why the 
model fails to predict the correct behaviour by comparing it with Stoner’s theory 
which agrees well with experiment. The basic assumption of Stoner’s theory 
is that the total exchange energy is a certain function of the magnetization, but 
not explicitly a function of the temperature. Clearly a sufficient condition for 
this result to be true is that the exchange interaction between two electrons, one 
in state k, the other in state k’, is independent of k and k’. It is interesting to 
notice, however, that the condition is also necessary; this result is demonstrated 
in Appendix B. This suggests that the difference between our results and those 
of Stoner is due to the strong dependence of the exchange integral on k and k’, 
namely as |k—k’|-*. For a satisfactory theory it would be necessary to derive 
an exchange integral independent of, or at any rate much less dependent on k 
and k’ than (2.2). Any correction to the exchange integral must however be 
sought outside the scheme of the one-electron approximation. 

A very likely reason for the unsatisfactory nature of our model is the use of 
the one-electron approximation, which does not take account of positional 
correlations between the electrons other than those implied by the Pauli principle. 
Now we have already mentioned in §1 the difficulties which arise in connection 
with the specific heat of sodium when the Bloch exchange terms are included. 
Wohlfarth (1950) has shown how these may be removed by replacing the Coulomb 
potential integral by a screened potential 


Virg)=ere exp (=A) ie eee eee (535) 


while the one-electron character of the wave function is retained. (‘This 
suggestion was first made by Landsberg (1949) to explain the width of the tail 
of the soft x-ray emission curve for sodium.) Modification of the Coulomb 
potential in this way is regarded as an empirical way of taking into account the 
positional correlations of the electrons. We may therefore enquire whether 
a similar modification of the potential would not remove the difference between 


our results and experiment. ‘The potential (5.1) gives as the modified exchange 
integral 
47e? 1 


ae aaa | P+? ele sieireke (5.2) 


which has a weaker dependence on k and k’ than the Bloch integral (obtained 
by putting A=0). It is clear then that by choosing suitable values of A, the 
dependence of (5.2) on k and k’ could be made sufficiently weak to give a second 
order transition from ferromagnetism to paramagnetism in the way predicted by 
Stoner’s theory. It does not seem likely, however, that calculations with (5.2) 
would be of much significance since (5.1) is merely an empirical expression 
without a theoretical basis. In addition the formula is open to criticism on 
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qualitative grounds. ‘The correlations we are considering are of short range and 
arise from the terms in the total Hamiltonian representing electron repulsions. 
The ‘correlations’ taken into account by the screened potential (5.1) would seem 
to represent long-range effects, since the term exp (—Ar 2) effectively reduces the 
product of wave functions in the exchange integral only at large separations 74p. 
It is difficult to see how these long-range effects are to be interpreted (see, however, 
a note by Bohm and Pines (1950)). In view of this, further consideration is 
being given to the problem of positional correlations. 
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AP PEN ea. 


We give here the details of the calculation of the thermodynamic functions for 
a perfect Fermi—Dirac gas using the approximate forms (2.9) and (2.10) for the 
distribution function. When there are no exchange interactions there can be no 
ferromagnetism: we therefore obtain an approximate expression for the free 
energy of the perfect Fermi—Dirac gas by putting x;=x,, €=0 and ¢«;=0 in 
equation (2.15). ‘This yields 

1 eS) 3 kT 

Ne, => 5 A(x) — 2 G C(x). @ el ejelcle (A 1) 


In this equation 3.A(x)/5 is the reduced internal energy E(x)/Neg, and 3C(x)/2 is 
the reduced entropy S(x)/Nk. ‘The equilibrium value of x is that for which 
F(x) isa minimum; this gives the equation 


RT Jey =2A'(x)/5C’(x). ere (ND) 


From the tabulated values of A(x), C(«) and their derivatives it was therefore 
simple to obtain F/Ne) and RT/e, as functions of x. Figure 2 was then con- 
structed from these results. We give below the forms taken by the functions 
at very high and very low temperatures. 
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(a) Very low temperatures. From equations (2.12) and (2.13) we find, for 
x<1, A(w)=1+5x*/12 and C(x)=x. Therefore, by (A2), x=3R7/eq so that the 
internal energy F is given by 


wena lt p ety | 
Ne, 5 A ak 


whereas the accurate expression is 


ee SHB rece we 
Nea 5 (PAu, 


(b) Very high temperatures. The limiting forms of A(«) and C(x) for large 
values of x are A(x) =2.4?%x?/3/3 and C(x) =25/18+(2/3)Inx. Therefore 
x23 = 15RT.4-2/?/4e, which gives for the internal energy and entropy respectively 


E=3NkT/2 and 
on 25 5 /15\22 
wee] tat La(Z) 


ay ae die) J OVAUE 
pons li ~ aad 
In this equation-we have substituted for ¢y from equation (2.16). We note that 
whereas the expression for the internal energy is accurate, the entropy is in error 
by a constant amount, the accurate expression being 
S ae (2rmkT)3? 2V ie 2) 
Nk N h Jigs 


the difference in the constant terms is, however, only 0-104. 


as a 390} 


APPEND IEXSE 


We shall find below the conditions under which the.sum of the exchange and 
Coulomb energies for the whole assembly of electrons depends only on the 
magnetization and not on the temperature. The individual Coulomb and exchange 
integrals which appear in equation (2.1), we denote by C,,, and Jy, respectively. 
Let the distribution function for the sub-assembly of parallel spins be f,(k) and 
let it be f,(k) for the anti-parallel spins. ‘The total exchange and Coulomb ‘ 
energy is then 


G= Seats cc + 3 ee Z folk folk ict’ 
i > fulPfUR Cae + E22 fafa 
= (A(R) Geese tet alt ns PRE Aes: (B 1) 


k k’ 

In this equation we have introduced implicitly a self-exchange integral, so 
that the inner summations in the first four terms may be over all allowed k’ vectors, 
not only over all k’~k. Clearly, if we make J, = — Cy, the self-exchange terms 
in the first and second summations will cancel the self Coulomb terms in the third 
and fourth. ‘The last summation in (B1) represents the Coulomb energy of the 
parallel spins with the anti-parallel spins; terms in C,, will therefore occur here 
in any case. 


: 
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Now if there are n, parallel spins and n, anti-parallel spins, f,(k) and f,(k) 


_ must satisfy the conditions 4, f,(k) =m, and X,f,(k) =m. It is given that G is 
independent of temperature, and therefore the variation 5G is zero, for all 
) variations in f, and f, 


6G =0= X4f,(k) UR Sia + Ldf2(k) Ufeo(h)S ict 
+ Vf (Rk) UHR) Cw + & Ofo() Ufo(k') Crac 
+E 8A Efk Cue + D5falh) Dhl’) Coes 


subject to the conditions 


Ba Stas (B 2), Beek) =O. And pet (B 3) 
Introducing the undetermined Lagrange multipliers A, and A,, we have therefore 
Oe NTP ie 1h Cea Salk Cu} 
PROB) at DL falh Tite + Sal ® Cae $F (F')Crae]} =0. 


As this must be true for all arbitrary variations 6f, and 5f, we must have 
eI a + Cy) rth eho Eee iat yd | Wteliaraacc (B4) 


& Sal) Fra + Cys) + © fi(F Cre ere Geer (B5) 


These two equations show that if the energy G is not an explicit function of 
temperature, then the total Coulomb and exchange energy of a single electron 
with the remaining electrons of the assembly is likewise independent of 
temperature. Equation (B4) expresses this for an electron with parallel spin 
(in the f, distribution) and (B5) for an electron with anti-parallel spin. Now 
take (B4) and consider further variations df, and 6f,, obtaining 


a OF (RS scree + Chet) + = Ofo(R’) Cun: = 0. 


Combining this with the conditions (B2) and (B3) by the aid of further 
undetermined Lagrange multipliers j, and j,, we have 


2 AiR (Siac eC us 9) = Ofo(R’)( Cray +Jo) = 9. 


As this also must be true for all arbitrary variations 6f, and 6f, we must have 
Cae = —Jo and Jig + Cy = —J, which can be put’in the form Cy,=C and 
Jy, =J. If we now substitute these results into equation (B1), we find that 
the total Coulomb energy is simply N?C/2 (where N=n,+ n,) and is therefore 
independent of magnetization. Likewise the total exchange energy is 
1N?J +iN?*JC? where, as before, ¢ is the reduced magnetization. We have thus 
shown that if the sum of the exchange and Coulomb energies is to be a function 
only of the magnetization, then, within the scheme of the one-electron 
approximation the total Coulomb energy is a constant, whilst the total exchange 
energy can only have the form assumed by Stoner. Asa corollary it follows that 
if the exchange energy is assumed to be some other function of ¢, it must contain 
parameters which depend on the temperature. 
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The Velocity of First Sound in Liquid Helium 
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ABSTRACT. ‘The velocity u, of ordinary sound in liquid helium has been measured in 
the temperature range from 1:2° K. to 4:2°x. Its value extrapolated to 0° kK. is compared 
with the velocity uw, of second sound at 0:1° K. and a slight departure from the Landau 
relation u,=1w,/ 1/3 is indicated. The behaviour of u, in the neighbourhood of the A-point 
has been studied in detail. 


Silky JUNAURONDNOKC TEI KON) 


STUDY of ordinary sound in liquid helium is of interest for several reasons. 
A First, a determination of the velocity at low temperatures would make it 
possible to test more critically the relation u,=u,/\/3 between the 
velocities of first and second sound at absolute zero predicted by Landau (1941). 
Measurements of second sound velocity at temperatures down to 0-1°K. have 
recently been completed (Atkins and Osborne 1950), and it was desired to extend 
the existing first sound measurements to lower temperatures in order to make 
possible a more accurate extrapolation to absolute zero. Second, the behaviour of 
the velocity in the immediate vicinity of the A-point is of interest. The Ehrenfest 
relations for a phase transition of the second order predict a discontinuity of 24%. 
at the A-point, using accepted values of the thermodynamic parameters involved. 
Previous investigators have found no evidence for such a discontinuity under the 
saturated vapour pressure (Findlay, Pitt, Grayson Smith and Wilhelm 1938), 
although it has been observed at higher pressures (Findlay et al. 1939). A know- 
ledge of the exact shape of the curve in this region might lead to a better under- 
standing of the mechanism of the d-transition. Third, the question of dispersion 
has never been adequately investigated. It is probable that there is no large 
dispersion in the region between 1 and 15 Mc/s. (Pellam and Squire 1947), but a 
small effect of a few per cent is not excluded in this range and nothing is known 
about the behaviour at other frequencies. Finally, the attenuation is of interest, 
particularly near the A-point and in the low-temperature region, as it may 
provide information about the internal structure of the liquid. 

Measurements of velocity were made by Findlay et al. (1938), using a 
continuous-wave resonance method at a frequency of 1-338 Mc/s., and of both 
velocity and attenuation by Pellam and Squire (1947), using a pulse technique at 
15 Mc/s. The present work is in agreement with these previous results in the 
temperature range which they covered (from the normal boiling point down to 
1:6°x.). In addition velocity and attenuation measurements have been extended 
to 1:2°xK., and a much more detailed investigation of the behaviour in the neigh- 
bourhood of the A-point has been made. 
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§2. DESCRIPTION OF APPARATUS 


A pulsed ultrasonic method similar to that employed by Pellam and Squire was 
-used. Pulses of approximately ten microseconds duration from a 14 Mc/s. 
| oscillator were fed into a quartz crystal immersed in the liquid helium. The 
_ resulting ultrasonic pulses were reflected from a polished brass reflector, and the 
echoes received by the same crystal. After amplification and demodulation the 
echoes were displayed on a Cossor Model 1035 double-beam oscilloscope, together 
with timing pips provided by a standard radar unit (a ‘ Gee’ unit). 

Figure 1 shows the etalon with which the measurements were made. An 
X-cut quartz crystal G, 1:9cm. in diameter, was mounted in a specially designed 


A 
A’ 
Figure 1. 
Etalon. A: Polished reflector. F: Brass electrode. 
B: Aligning screws. G: Quartz crystal. 
C: Beryllium—copper spiral springs. H: Coaxial cable. 
D: Beryllium-copper springs. I: Dystrene insulator. 
E: Brass case. 


holder. ‘The top surfaces of the electrode F and brass case E were machined 
accurately plane, and the crystal was held against these surfaces by means of three 
beryllium copper springs D. The top face of the crystal was silvered, and contact 
with this face was made through these springs. The assembly was so designed 
that thermal contraction would not alter the alignment of the surfaces against 
which the crystal was pressed, thus making certain that the central electrode 
would remain in contact with the crystal face and that no stresses would act upon 
the quartz when the apparatus was cooled. The electrode was made in the shape 
of acone to eliminate spurious echoes caused by the sound travelling back from the 
crystal. Such echoes would have to undergo many reflections before ever 
reaching the crystal again, and would be highly attenuated. The electrode 
was connected to the oscillator by a coaxial cable H extending to the cryostat head. 
The crystal holder assembly was attached to the reflector A by three aligning screws 
B, and was held in position by the beryllium copper springs C. By immersing the 
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apparatus in water and adjusting these screws to maximize the echoes, it was 
possible accurately to align the reflector and transducer, and it was then found that 
no appreciable misalignment occurred when the apparatus was cooled to liquid 
helium temperatures. The reflector was of brass, and was machined to be as 
plane and smooth as possible. Two alternative reflectors (A’ and A”) were 
provided in order to obtain longer or shorter path lengths. By interchanging 
them with A and varying the position of the aligning screws, it was possible to 
obtain a number of different path lengths. : 

All the electronic equipment was synchronized by means of the ‘ Gee’ unit, 
which provided trigger pulses at a repetition rate of 500 per second and calibrator 
pips spaced at intervals of 62, 66%, and 3334 microseconds. By the use of separate 
variable delay circuits to trigger the modulator and the oscilloscope time base, 
it was possible to move the pulse with respect to the calibrator pips and also to 
examine any part of the time base on an expanded scale. ‘The small correction 
which had to be applied for the non-linearity of this expanded scale was determined 
by photographing a pulse of fixed duration at anumber of positions along it. The 
capacitance of the crystal and holder was tuned out by a small variable inductance in 
parallel with it. ‘The output of the crystal was fed through a pre-amplifier and 
then through a 14 Mc/s. amplifier with a bandwidth of 4Mc/s. ‘The amplifier 
output was demodulated and fed on to the vertical plates of the oscilloscope through 
a low-gain video amplifier passing frequencies up to 2 Mc/s. 

The ultimate accuracy of the time measurement was limited by the uncertainty 
in the positions of the feet of the echo and the input pulse. Because of the poor 
impedance match from the crystal into the liquid, the transducer had a small 
bandwidth and the echo was distorted so that its rise time was about 2 microseconds. 
The apparent position of the foot of the echo therefore depended on the signal-to- 
noise ratio, and this dependence was investigated by varying the echo height at 
constant temperature by means of an attenuator. When the echo height was 
approximately equal to the noise level this pulse amplitude correction amounted to 
0-7 microsecond, but after it had been applied the position of the echo could be 
estimated to within 0-5 microsecond visually and 0-1 microsecond from enlarged 
photographs. ‘The uncertainty in the position of the foot of the input pulse arose 
because this pulse was many orders of magnitude larger than the echo, so that any 
small build-up region at its foot could have been unduly emphasized and could 
conceivably have introduced an error as large as 0-5 microsecond. The length of 
the etalon could not be measured to better than 0-05 mm., or $%. Acorrection of 
4° was applied for contraction of this length upon cooling, but this was only an 
approximate estimate. Asa consequence of all these sources of error it is believed 
that relative values of the velocities are accurate to within 4%, but that the 
absolute values may be in error by as muchas 1%, 

Temperature measurements were made with a mercury manometer at 
pressures above 4cm. Hg, and with a butyl phthalate manometer at lower pressures 
(the accuracy being of the order of mm. in both cases). A differential butyl 
phthalate manometer was also provided to register small changes of pressure, and 
proved to be of particular value near the A-point. By comparing the oil and mercury 
manometers at one or two temperatures the conversion factor relating the density 
of oil to the density of mercury was obtained. Above the A-point a correction was 
applied to allow for the temperature gradient in the helium bath due to the 
hydrostatic pressure head. This will be discussed in more detail in §5. The 
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A-point, as determined by the sudden cessation of boiling, occurred at a pressure of 
3-77 cm. Hg, which corresponds to a temperature of 2:179°k. on the 1949 scale 
(van Dijk and Shoenberg 1949). Since the accepted value (2:186°k.) of the 
A-transition temperature is probably not very accurate, the observed deviation from 
it is not very significant. In any case, although this deviation may possibly be 
accounted for by errors in manometry which have been overlooked, temperature 
differences measured from the d-point are certainly correct to within 0-001°K. 
In order to obtain very accurate temperature control near the A-point, a constantan 
heater was wound round the etalon, and the power fed into it varied in such a way 
that a balance was maintained between the amount of liquid helium evaporating 
and the speed of the pumping system. In this way the temperature could be kept 
constant to better than 10~¢ degree. 


5. DHE VELOCIRY-LEMPERATURE CURVE 

The velocity as a function of temperature is shown in Figure2. The circles and 
the full curve represent the present results; the broken curve is the one given by 
Findlay et al. as the best fit to the 1-338 Mc/s. data. For the sake of clarity the 
results of Pellam and Squire are not shown, but within the limits of their experi- 
mental error they are consistent with both sets of data. The difference between 
the two curves in the neighbourhood of the A-point may indicate the presence of 
some dispersion in this region, particularly in view of the excellent agreement 
above 3°xk. It should be emphasized, however, that the disagreement is very 
little larger than the sum of the two experimental errors. 
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Figure 2. Figure 3. 
Velocity of sound in liquid helium. Velocity of sound in liquid helium (low 
—~-— Findlay, Pitt, Grayson Smith and Wilhelm. temperature range). 
—o-o-— Present results. Path length=4-310 cm. The vertical line indicates the size of the 
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In Figure 2 there is a faint suggestion of a maximum near 1:3° K.; however it 
was believed that the apparent dip at the lowest temperatures was due to the small 
size of the echo there and the resulting difficulty in estimating the position of its 
foot.* Measurements were therefore made with a shorter etalon of length 4mm., 
in order to have larger echoes. Figure 3 shows the results of these measurements, 
which were made photographically on the second echo, giving an effective path 
length of 1-6cm., and a relative accuracy of within+0-2%. Since the error in 


* The attenuation rises very rapidly below 1:4° x. This will be discussed in greater detail in a 
later paper. 
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measuring the length of the etalon was more important in this case, the absolute 
value has been adjusted to agree with the more precise absolute determinations 
made with the longer etalon. From this curve it may be seen that there is no evi- 
dence for a maximum in the velocity at low temperatures. The value extrapolated 
to absolute zero is 237 + 2m. sect. 

Figure 4 shows the adiabatic compressibility 8,4 calculated from the usual 
formula B,4=1/pu,°, where p is the density and w, the velocity of sound. This 
calculation is based on the assumption that there is no dispersion between zero 
frequency and 14Mc/s. Values above 2:5° x. have been given by Findlay et al. 
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Figure 4. 
Adiabatic compressibility of liquid helium at 14 Me/s. 


$4. VELOCITIES OF FIRST AND SECOND SOUND NEAR O°K. 

If the above extrapolation of the velocity of first sound w, from 1:2°k. to 0°K. is 
justifiable, then Landau’s prediction is that the velocity of second sound at 0°K. 
should be uy =u;/4/3 =137 +1m.sec™, which is to be compared with the value 
Uy =152+5m.sec at 6-1°K. given by Atkins and Osborne (1950). The second 
sound measurements were only exploratory and not particularly accurate, and 
final judgment should be reserved until second sound has been more fully 
investigated below 1°k., but there is obviously some evidence for the existence of a 
discrepancy of the order of 10%. ‘This discrepancy is not large enough to 
raise doubts about the fundamental principles of Landau’s theory, but it does 
stimulate a search for reasons why Landau’s relation might not be accurately 
obeyed. The relation arises because, at sufficiently low temperatures, the 
important elementary excitations in Landau’s theory are phonons, or standing 
sound waves, and the relationship between their energy « and momentum / is 
€=u,p. The thermodynamic functions such as entropy and internal energy can 
be derived from a knowledge of the elementary excitations, as in the Debye theory 
of the specific heats of solids, and the velocity of second sound can then be expressed 
in terms of the thermodynamic functions, giving the Landau relation u, =u,/1/3. 
It is important to notice, however, that at a temperature T the average frequency of 
the phonons is of the order of RT/h, which is still as high as 2 x 10®c/s. at 0-1°K. 
The discrepancy therefore suggests that the velocity of the phonons (i.e. the 
velocity of first sound) at a frequency of 2 x 10°®c/s. is greater than the velocity 
measured at 1-4 x 10%c/s. Dispersion always occurs when the absorption is so 
large that the wave is appreciably attenuated in a distance comparable with the 
wavelength, and present knowledge of the attenuation indicates that this may well 
be the case at a frequency of 2 x 10° c/s. and a temperature of 0-1°K. This will be 
discussed in greater detail in a later paper describing some results on the 
attenuation of first sound. 
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§5. BEHAVIOUR AT THE A-POINT 

The question of the behaviour of the velocity of sound near the A-point is an 
interesting one. If the transition is assumed to be second order and the best 
possible values of the various quantities are taken to be those given by Keesom 
(1942), the Ehrenfest relations predict a discontinuity of approximately 24% in the 
velocity of sound. The experimental evidence on the magnitudes of the dis- 
continuities in specific heat and thermal expansion is so inadequate, however, that 
this prediction is very approximate, and it is possible that a discontinuity much 
larger than this or even zero might still be consistent with the Ehrenfest relations. 

The measurement of velocity near the A-point is made difficult by the extremely 
high attenuation in this region (Pellam and Squire 1947). Itis therefore necessary 
to extrapolate the curve over an interval of a few thousandths of a degree, and since 
the velocity depends strongly on the temperature in the vicinity of the A-point, 
considerable uncertainty is introduced in this way. The use of a shorter etalon 
reduces the attenuation and makes it possible to approach closer to the A-point, but 
the shorter time delays involved make the measurements less accurate. It is there- 
fore necessary to compromise between these two sources of error. The procedure 
adopted was to approach as closely as possible using a long etalon for maximum 
accuracy, and then to make supplementary measurements with a much shorter 
etalon. The points obtained with the shorter etalon were then adjusted in 
absolute value to lie evenly about the more accurate curve in the region of overlap, 
and were sufficiently accurate to determine the shape of the curve very near the 
A-point. 

Another factor which made it difficult to come very near the A-point was the 
existence of a temperature gradient in the liquid above the A-point due to the 
hydrostatic pressure head. ‘This effect is mentioned by Schmidt and Keesom 
(1937) and has been investigated more recently by Lock, Pippard and Shoenberg 
(1951), but its existence does not seem to be generally recognized. ‘The thermal 
conductivity in liquid helium I is so poor that whenever there is a heat flow upwards 
there is a tendency for a large temperature gradient to be set up; however, the 
temperature at any point cannot rise above the temperature corresponding to a 
vapour pressure equal to the pressure at the surface plus the hydrostatic pressure of 
the column of liquid above the point, otherwise the liquid would vaporize. ‘The 
temperature at a point below the surface is therefore greater than the temperature 
at the surface by an amount determined by the hydrostatic pressure and the slope of 
the vapour pressure curve. ‘This amounts to somewhat over 0-001 deg/cm. near 
theA-point. When, however, the surface of the liquid reaches the A-point, it is no 
longer possible for this temperature gradient to be maintained because of the high 
thermal conductivity, and the entire volume becomes superfluid at once. This 
leads to a sudden jump in temperature, which may be as great as one or two 
hundredths of a degree near the bottom of the flask. The exact mechanism by 
which the transition takes place is obscure, and needs further study. The present 
work afforded an opportunity to demonstrate the existence of this hydrostatic effect 
in avery convincing fashion. Measurements of velocity were made both above and 
below the A-point, first with a pressure head of about 20 cm. and then with a much 
lower pressure head of approximately 4cm. The data from this experiment are 
shown in Figure 5, where velocity is plotted against surface temperature determined 
manometrically. ‘The agreement of the two sets of points below the A-point 
indicates that the etalon temperature was then the same as the surface temperature, 
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and hence no hydrostatic temperature gradient was present. Above the A-point, 
however, the two curves are clearly separated, the upper curve corresponding to 
the higher pressure head. If this curve is shifted to the right by 0-019°k., the 
temperature difference resulting from a pressure head of 16 cm., it agrees well with 
the lower one. Because of the existence of this effect, measurements above the 
A-point were made with the liquid level only a few millimetres above the reflector, 
and the theoretical correction was then applied, measuring the pressure head from 
the middle of the etalon path. In this way it was possible to approach within 
0-002° k. of the A-point from above. 
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Figure 6 shows the results of measurements near the A-point, after all corrections 
have been applied. For all the points on this curve the echo was large and the 
pulse amplitude correction discussed in § 2 was completely negligible. One of 
the most striking features of this curve is the sharp dip in the velocity within the 
last 0-02° k. above the A-point, and a corresponding tendency for the velocity to 
decrease rapidly as the A-point is approached from below. It is obviously 
impossible to say much about the discontinuity at the A-point; in fact it is 
conceivable that both branches drop to very much lower values than have yet been 
observed. It should be remembered, however, that the Ehrenfest relations 


predict a discontinuity in the isothermal compressibility 8, which is related to 
the velocity of sound by the formula 


ae! OV mee 
p= — (55) = ae 
where y, the ratio of the specific heats, can be expressed in terms of the coefficient of 
thermal expansions: y—1= Tuj’«?/C,,ify=1. yis a rapid function of temper- 
ature in the neighbourhood of the A-point and may be entirely responsible for the 
peculiar behaviour of w, in this region. Calculations of 8 from u, using the rather 
inaccurate values of y at present available suggest that 8 increases monotonically 


over the entire temperature range and has a discontinuity of at least 1°/, at the 


\-point, but a firm decision on this issue must await more accurate measurements of 
the coefficient of thermal expansion «. 
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The Flow of Liquid Helium II through Wide Capillaries 


By’ K. R. ATKINS 
Royal Society Mond Laboratory, Cambridge 


Communicated by D. Shoenberg;, MS. received 29th May 1951 


ABSTRACT. A study has been made of the isothermal, pressure-induced flow of liquid 
helium II through glass capillaries of bore between 2:6 x 10-3 and4x10-%cm. The critical 
velocity is found to be less than 1 cm. sec! except, possibly, in the case of the finest capillary. 
The theory that there is a force of mutual friction between the normal and superfluid 
components predicts the correct order of magnitude for the flow velocities, but there are 
discrepancies which suggest the existence of some other source of friction. 


§1. INTRODUCTION 


HEN liquid helium II flows through very narrow channels, or through 

\ x / the helium film, at velocities less than a certain critical value v,, the 
frictional forces opposing the flow are very small and may even vanish 

(Daunt and Mendelssohn 1946, Atkins 1950 b, Bowers and Mendelssohn 1950). 
Above this critical velocity, however, the frictional forces increase rapidly and 
appreciable pressure gradients are necessary to produce even a small increase 
in velocity. Similar phenomena are observed when the flow is produced 
thermally (Kapitza 1942). Velocities up to v, are produced by temperature 
differences too small to be measured, but as soon as u, is exceeded the required 
temperature difference begins to increase rapidly. However, in channels of 
width greater than 10-?cm. the situation is entirely different (Allen and Misener 
1939). In such channels there is, as yet, no experimental evidence for the 
existence of frictionless flow and the critical velocity, if it exists, must be very 
small. The velocity of flow varies in a markedly non-linear manner with the 
pressure gradient and the frictional forces obviously have a complicated and 
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unprecedented character. Although superfluidity is more clearly manifested 
in the frictionless region below the critical velocity, the manner in which it 
breaks down in the imperfect region above the critical velocity may provide 
valuable clues concerning its nature. 

One aspect of the problem, thermal conduction in capillaries, has already 
been extensively studied and, apart from minor discrepancies, is satisfactorily 
explained in terms of the two-fluid theory if one assumes that there is a mutual 
force between the two components varying approximately as the cube of their 
relative velocity (Gorter and Mellink 1949). The present paper extends the 
experimental evidence by a detailed study of the case of isothermal, 
pressure-induced flow through glass capillaries. Still further evidence is 
provided by the experiments of Hollis-Hallett (1950) on the damping of the 
motion of a disc performing large-amplitude torsional oscillations in the liquid. 
In a later paper these various lines of evidence will be compared in an attempt 
to deduce the nature of the frictional forces and to decide whether the 
Gorter—Mellink theory of mutual friction is adequate to explain all the results. 


§2. DESIGN OF THE-EXPERIMENTS 

An experimental survey of the whole field of flow phenomena in capillaries 
has already been made by Allen and Misener (1939) and it is clear from this 
work that the situation is very complex. ‘The present research was therefore 
mainly concerned with a more detailed investigation at the lowest temperature 
easily obtainable (about 1:2°x.), because at this temperature the normal 
component was only about 3% of the whole, so that one was studying the flow 
of practically pure superfluid and the.results could be more easily interpreted. 
Even if the flow of the normal component had been restricted only by its own 
viscosity, of the order of 40 micropoises (Andronikashvili 1948), its contribution 
to the total flow could not have amounted to more than 1°% except in the case 
of the widest capillary when it might have amounted to 5% at the highest pressure 
gradient. ‘To obtain some idea of the variation of various quantities with 
temperature, the measurements were repeated at 1:5°k., which was the highest 
temperature at which it was reasonably certain that the flow was not being 
influenced by the thermal effects discussed below. 

Figure 1 shows the three types of flow vessel used. The rate of flow of liquid 
through the capillary C was measured by taking cathetometer readings of the 
meniscus in the monitoring capillary M which had an internal diameter of about 
0:-5mm. ‘This procedure was preferable to making direct observations of the 
meniscus in the main reservoir R, because the smaller meniscus in M was more 
sharply defined and also the surface tension rise in M kept its meniscus clear of 
the bath meniscus, enabling reliable observations to be made at small pressure 
heads. Capillaries longer than 10cm. were accommodated either in the form of a 
spiral (Figure 1 (a)) or as six adjacent straight portions joined together by much 
wider tubing (Figure 1(0)). The relative advantages of these two arrangements 
will be discussed in §4(i). The constriction at F reduced transfer of liquid 
through the film, but, because of certain considerations discussed below in 
connection with the thermo-mechanical effect, it was not permissible to make this 
constriction so narrow that the film flow correction could be entirely neglected. 
Therefore, in order that this correction might be reliably estimated, it was 
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necessary to avoid the increased rates of flow resulting when the film is formed on 

a surface contaminated by solid deposits (Bowers and Mendelssohn 1949), and 
the experiments were therefore performed inside a cryostat specially designed to 
avoid such contamination (Atkins 1950a). The correction was then based on 
_ the assumption of a rate of film flow equal to 8 x 10-°cm? sec"! cm™, which was 
confirmed in the case of two of the vessels by studying the flow when the end of 
the capillary C was sealed off. Because the experiment was performed inside a 
_ closed chamber containing a constant amount of liquid, there was no correction 
__ to the pressure head arising from evaporation effects, and the small change in the 
position of the outer level accompanying a change in the inner level could be 
calculated accurately from purely geometrical considerations. 
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(a) (b) (C) 
Figure 1. The vessels. 


Thermal effects associated with the flow are sufficiently important to justify 
a detailed calculation of their magnitude. ‘The mechano-caloric effect results 
in a heat input to the reservoir R of 7r?ovQ* erg sec~1, where 7 is the radius of the 
capillary C, p is the density of the liquid, v is the mean velocity of flow and O* is 
the amount of cold transported per unit mass of liquid, which is unknown but is 
certainly not greater than T'S. ‘This heat is dissipated by a distillation process 
in which liquid evaporates from the reservoir and is condensed again in the bath. 
Evaporation from the reservoir surface can take place only if there is a pressure 
discontinuity Ap, between liquid and vapour and simple kinetic theory 
considerations show that the rate of evaporation is 


1 /3M\12 ae 
i(Rr) IND OH GC CEE Noa te sree t ce (1) 


M is the molecular weight, R is the gas constant, T the temperature, A the 

cross-sectional area of the reservoir, and we have assumed that the vaporization 

coefficient is unity. Equating the heat input to the required latent heat 
(L ergs per gm.) we obtain 

Anr’upQO* (RT \12 

Ap = 4rrr*vpQ* , 
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There is also a pressure drop Ap, across the constriction F arising from the 
viscous flow of gas through it. If this constriction has radius a and an effective 
length A and the gas has a density p, and a viscosity 7, then, using Poiseuille’s 
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formula to relate the volume of gas evaporating per second to the pressure drop 
across the constriction, 
1/3M\12A,_. za’ 
Z (ar) ae Sn) Ap.» i 
pea 20M S00 eA ea iene a 
Aps ~ zat RT Pe, 
As the cross-sectional area of the bath is so large, the pressure discontinuity 


there may be neglected, and so the temperature difference needed between the 


inner and outer levels to produce the pressure difference between them is 
AT = (Ap, + Ape)/(dp/4T )o, 
ee 
~ (dp/dT),LA\3M Cae SAR MERA Prt ee 
(dp/dT), being the gradient of the vapour pressure curve. Because of the 


thermo-mechanical effect this temperature difference is equivalent to a 
hydrostatic pressure head 


me 
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The uncertainty in the relevant value of O* means that this correction cannot 
be reliably estimated, and so the apparatus must be designed in such a way that 
it is negligible. This was achieved by using as large a reservoir as possible 
(large A) and choosing the diameter of the constriction F so that the second term 
was no larger than the first. ‘The investigation was in any case restricted to 
temperatures below 1:5°k. where the thermo-mechanical effect is small. Under 
these circumstances it can be shown that, even if O* is given its maximum 
value T'S, the correction was entirely negligible, except at pressure heads less 
than 10-*cm. of He, when the situation was, however, also complicated by the 
presence of oscillations (see Table 1). 

An alternative approach to the problem would be to produce the flow 
thermally by deliberate supply of heat and to measure both the temperature 
difference and the hydrostatic pressure head, but this method was not adopted 
for two reasons. It is not yet certain that a temperature difference is the exact 
equivalent of a pressure head, and, in any case, the conversion factor between 
the two is not yet known. Moreover, if this conversion factor should turn out 
to be the one applicable to the thermo-mechanical effect in narrow slits (see 
equations (12) and (13)), then the smallest measurable temperature difference 
would correspond to a pressure head of about 5 x 10-%cm. of helium, whereas 
the true hydrostatic pressure head can be measured to an accuracy of 5 x 10-3 cm. 
of helium and it is thought that the thermo-mechanical effect can be made as small 
as 4x 10-*cm. of helium by suitable design of the apparatus. 

In the case of the very finest capillary of bore 2-62 x 10-cm., any permissible 
outlet at the top of the reservoir would have introduced an intolerably large film 
flow correction. ‘Therefore, as shown in Figure 1 (c), the top of the reservoir 
was sealed off by a large thin-walled glass bulb B with a surface area of about 
20cm?. ‘Thermal contact between the inside and outside of the reservoir was 
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achieved by thermal conduction through the glass wall of this bulb. The figure 
quoted for the thermo-mechanical effect in Table 1 is based on the thermal 
conductivity of glass at 1-3°K. as measured by Keesom (1945). 

To minimize effects due to heat entering the cryostat from outside, the vessels 
were illuminated by sodium light from which heat radiation had been removed 
by a filter of Chance ON19 glass and the efficacy of this procedure was proved 
by the fact that rates of emptying were equal to rates of filling to within the 
experimental error (see Figure 2). 


§3. EXPERIMENTAL DETAILS 
Table 1 contains a complete description of the various capillaries and their 
associated reservoirs. The internal diameter of a capillary was measured at 
various points along its length with a travelling microscope, and deviations from 
the mean are given in case these should prove to be important when the flow is 
better understood. The film flow correction is calculated on the basis of a rate 


: Table 1 
Capillary I II Ill Ges 
Mean int. diam. (cm. x 10?) 4-40 2:03 0-815 0-262 
Variation of int. diam. (cm. x 107?) 0-10 0-10 0-013 0-01 
Cross section of reservoir (cm?) B25) 0-892 0-201 0:00571 
Capillary lengths used (cm.) 48:6, 8:04,1:75 46-6, 7-90 8:03, 0:83 7-76 
Form of longest capillary ; Zigzag Zigzag Straight Straight 
Diam. of film constriction (cm.) 0-2 0-1 0-09 0) 
Film flow corr. to velocity (cm/sec.) 0-03 0-08 0:55 0 
Thermo-mech. corr. (cm x 107-4 of He) 3:6 Des 2:0 ~2 


of film flow equal to 8 x 10->cm? sec-!cm~!. The thermo-mechanical correction 
is an upper limit obtained by putting Q*=T7'S in equation (5) and refers to a 
temperature of 1:2°K. The value at 1-5°K. is approximately four times as large. 
This correction is proportional to the mean velocity of flow through the capillary 
and the value quoted is for a mean velocity of 1 cm. sec}. 
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Figure 2._ Mean velocity as a function of pressure head. 
Bore of capillary=0-203 mm. 'Temperature=1:-22° k. 


The experimental procedure was to raise or lower the vessel until there was a 
level difference of about 1 cm. and then to time the meniscus in the monitoring 
capillary (M in Figure 1) as it traversed the scale in the eyepiece of a cathetometer. 
This scale had previously been calibrated against a scale engraved on metal and 
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found to be linear to within 1°%%. The two principal sources of random error were 
probably fluctuations in the bath temperature and optical flaws in the walls of 
the Dewar flasks, and it was necessary to repeat the procedure several times in 
order to obtain sufficient accuracy. Figure 2 shows the size of the random error 
and illustrates how an accurate curve could be obtained by smoothing out the 

data from six runs. . 


§4. CORRECTIONS 
(i) End Corrections 


When an ideal classica. liquid with zero viscosity enters a narrow pipe in 
which the velocity is v, there is a pressure drop of 4pv" at the entrance, arising 
from the Bernoulli effect, but there is also an equal and opposite pressure rise 
at the exit, so that the total end effect is zero. In the case of a viscous liquid, 
however, the effect of the frictional forces is to produce a non-zero total pressure 
drop, giving rise to the well-known Hagenbach correction (see, for example, 
Prandtl and Tietjens 1934). When the flow of liquid helium IT is subject to 
frictional forces, as it is in wide capillaries, one might expect to find a similar 
end effect which, if comparable with the Hagenbach effect, would involve an ~ 
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Figure 4. End effects. 
Temperature=1-2° k. 


Figure 3. Effect of varying the length 
of capillary I. 


important correction. Evidence for the existence of such an effect is provided 
by curves such as those in Figure 3. In this figure, x is the level difference and 
J the capillary length, so that in the absence of an end effect pgx// would be the 
pressure gradient, grad p. The mean velocity v is defined as the rate of flow in 
cm’ sec”* divided by the mean cross-sectional area of the capillary. ‘The mean 
velocity should be a unique function of grad p, independent of capillary length, 
but is obviously not a unique function of pgx//, the discrepancy going in such a 
direction as to indicate that pgw// is greater than grad p because of an end effect. 
To investigate the matter further’ and to confirm that this was really an end 
effect rather than, say, a thermo-mechanical effect, the vessels shown in the inset 
to Figure 4 were used. A capillary of internal diameter 0-48 40-01 mm. and 
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length 8-1cm. was attached to a reservoir of cross section 3-86 cm? with a film 
constriction of diameter 2mm. ‘The flow measurements were repeated with a 
bulb inserted in the middle of the capillary and the capillary length only slightly 
different (7-3.cm.). The bulb introduced extra end corrections and the pressure 
head needed to produce a given velocity of flow was considerably larger. Denoting 
the end correction by <(v)cm. of helium and assuming that in the capillary 
grad p=/(v) then, in the absence of the bulb, 


PRIS OeneO) ee ae (6) 

With the bulb : PE vg Glo Le OD) pan ke ws Paes (7) 
Eliminating f(v) and remembering that «x, and x, then refer to the same 7, 

pele — wa = 2 / bed leo). gate Pe (8) 


In Figure 4 «=(x2/l,—.%x,/1,)/(2/l,—1/l,) has been plotted against v?/g to give a 
straight line of slope 1-05 + 0-05, so the end correction is 1:05pv which includes both 
the entrance and exit of the capillary. 

A correction of this magnitude was applied to the results at all temperatures and 
forallcapillary sizes. Figure 3 proves that the correction is independent of temper- 
ature, because at high velocities where most of the pressure head is due to end 
effects, the curves at 1-2°K. and 1:5°k. coincide. ‘There is some doubt whether 
the correction is the same in the finest capillaries in which the flow may be tending 
towards the ideal ‘ superfluid’ character it manifests in the film, but the velocities 
are then small and the correction not very important. In addition to the end 
correction there is probably also a ‘ constriction correction’ to be applied at every 
point where the bore is slightly constricted or whenever the bore tapers. It is 
easily seen that very slight constrictions at very few places could produce an effect 
comparable with the end correction and it is unfortunately impossible to make any 
estimate of the effect. For this reason the zigzag type of capillary shown in Figure 
1 (d) is to be preferred to the spiral shown in 1 (a). When such a spiral was cut up 
and calibrated it was found that its bore had been appreciably distorted and this 
probably explains why such spirals always required anomalously high pressure 
heads to produce a given velocity. In the case of the zigzag capillary the end 
correction was multiplied by six, but could be estimated with tolerable accuracy, 
and it will be seen in Figures 5 and 6 that, after the application of the correction, 
the curves for different capillary lengths coincided. However, one must not 
conclude from this that the ‘constriction correction’ is unimportant, as this 
correction may be proportional to the length of the capillary. Finally, it should be 
mentioned that, although this end correction is important in the present experi- 
ments, it is completely negligible in experiments on thermal conduction through 
capillaries of comparable size. 


(ii) Acceleration Correction 


As the pressure head decreases the velocity continually decreases and the 
liquid always has a finite acceleration. ‘The frictional forces are therefore not 
quite balanced out by the pressure head and a small correction should be applied 
for this. ‘The required correction is easily shown to be 
mre Ov 
mead (9) 
Under all conditions in the present experiments this correction was smaller than 
the observational error and could be completely neglected. 
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§5. PRESENTATION AND DISCUSSION OF RESULTS 


The experimental results are collected together in Figures 5 and 6, and the 


region of small pressure gradients is shown in greater detail in Figures 7 and 8. 
The previous results of Allen and Misener (1939) show a similar behaviour, but 
their velocities are always rather lower than the present ones. ‘This may be a 
consequence of the fact that their capillaries were spiralled and therefore more 
susceptible to constriction effects. ‘The use of spirals in the present research 
always led to lower velocities and all results obtained with them have been rejected. 
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Figure 7. Variation of velocity in the region Figure 8. Variation of velocity in the region 
of small pressure gradients at 1:22° k. of small pressure gradients at 1:52° K. 
Capillary 1@ /=48-6 cm. Capillary III © /=8-03 cm. 


Capillary II -e- /=46-6 cm. Capillary [IV @ /=7-76 cm. 


It is important to notice that for velocities below 5 cm. sect when the end correction 
is too small to influence the matter, there is good agreement between the curves for 
different capillary lengths. ‘This could not be the case if the thermo-mechanical 
effect were important. 


In an earlier paper (Atkins 1950 b) it has been suggested that flow results can be 
conveniently discussed on the basis of the empirical formula 


U=U;(r)+v_(gradp)+vu;(gradp,r). «=» saeeee (10) 


The critical velocity v, 1s such that the frictional forces are negligibly small until v, 
is exceeded. wv, is believed to decrease as 7 increases. In Figures 5 to 8 v, would 
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appear as an intercept on the velocity axis. In the film v, is about 40 cm. sec“, but 
it is clear from Figures 5 to 8 that when 7 is greater than 4 x 10-cm. v, cannot be 
greater than Icm.sec"!. It is noticeable, however, that the curve for the 
0-0262 mm. capillary runs parallel to the curve for the 0:0815 mm. capillary but is 


_ displaced upwards from it by about 2 cm. sec~!, suggesting that v, is greater for the 


| 


finer capillary by this amount. It would be unwise to place too much reliance on 


the precise values of the velocities obtained with the 0-0262 mm. capillary, as this 


capillary was very difficult to clean and to keep free from blockages, but it is 
difficult to escape the conclusion that the velocity was greater in this very fine 


capillary, particularly as precautions were taken to eliminate film flow effects (see 
§2). The product of critical velocity and radius seems to be at least 


2x 10-%cm?sec~1, as compared with a film transfer of 8 x 10->cm? sect. A more 


detailed and careful research into capillaries of this size would probably be very 
fruitful, but it would be necessary to select very uniform capillaries and to develop a 


—technique for cleaning them and measuring their bores accurately. 


0 5 10 5 0 05 10 15 20 
Grad p (dyne cm~3) . Log v 
Figure 9. Difference curves. Figure 10. A as a function of velocity. 
Capillary 1 @M-227K. 11-527 K. Flow through a 0:0815 mm. capillary 
Capillary, I (a 1222° 1K he 522 Ke Salis 22 skews @ li Drake 
Capillary, BV /\ 1-222 x. @ 1-52? x. Thermal conduction in a 0-095 mm. capillary 


ate dlagiane 1:61° kK. 


The component v, of the velocity can be eliminated by subtracting curves for 
different radi: 


u(r’) — u(r) =0,(7') — v,(r) + v,(7’, grad p)—v,(r, grad p). ...... (11) 


This has been done in Figure 9 by subtracting the curve for the capillary of bore 
0-815 x 10-? cm. (capillary III) from all the others. It is immediately obvious that 
these difference curves are straight lines to within the experimental accuracy and 
that the curves for the two temperatures coincide. The component v, may there- 
fore be represented approximately over the whole of the temperature range by the 
empirical formula v;=«(r) grad p. 'To evaluate vy exactly would require a much 
more extensive study of the exact form of v, and vs, but it is fairly clear that for the 
capillary of bore 0-815 x 10-*cm. (capillary III) v, is small (Figure 7) and v, is also 
small (Figure 9) so that the curve for this capillary may be taken to indicate the 
nature of v,. Itis probable that v, varies approximately as (grad p)", but it is not 
possible to verify this relation more precisely. 

The theory of Gorter and Mellink (1949) makes quite definite predictions 
concerning the velocity of flow in wide capillaries and these predictions will now be 
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compared with the experimental results. The fundamental equations suggested 
by Gorter and Mellink are 


nF > ral gradp+p,SgradT+Apyp.(Yp—%)*, vee (12) 
an = ae 2 47,V2 13 
Po ape > Se ee eae ee TS Apep( 0, = 0s) aan VeUE seeeee (13) 


Apart from the basic assumption that the two-fluid theory is still applicable at 
velocities exceeding the critical velocity, the novel feature of these equations is the 
term Ap,p,(v,—v;)* which represents a force of mutual friction between the two 
components. For isothermal flow through capillaries these equations reduce to 


(puiplgradp = Apis, 0), eS ee (14) 
(Pu/p) grad p as paps ae D.) ty Vans TODO OD (15) 
Adding, rad pan-uss = ah et Ih SG! See (16) 


This is the usual equation for viscous flow and its solution is v, =(7’grad p/8,), 
whence, from equation (14), 


—grad *) 138 7? grad p 
vs= (—>——]) +— >) Re ee 17 
( Appn 870 oe, 
ga Coe eto. Tt! ee a eee (18) 
p p 

ps (—grad ey" ry? grad p 
= = (ot nnn wes 19 
p ( Appn 87 a? 


The first term obviously represents v, and the second term has the correct form to 
represent v3, but the Gorter—Mellink theory does not take critical velocities into 
account and so v, is naturally absent. For the 0-815 x10-?cm. capillary is 
approximately equal to v, and, by equating it to (p,/p)[(grad p)/App,]"?, A has been 
deduced and plotted against vin Figure 10. ‘Thesame graph shows the values of A 
deduced by Gorter and Mellink from thermal conduction experiments. There is 


Table 2 

Diameter of capillaries Viscosity (poise) from 

Temp. Difference Oscillating 
(CCK) r’ (cm.) r (cm.) curves discs 

PS I. 4-40x10-2. 0-815 «107? 53 39 
(TSO? tx. 4-40x10- 0-815 x 107? 53 1s 
TaD aikKcs AMO Oasis KO 32 39 
1S Osi AMD SUNOS Oettiliss s¢ MO = 32 13 


a satisfactory order of magnitude agreement, but exact agreement is not to be 
expected in view of the possible errors of an experimental nature in both sets of 
experiments and the difficulty in eliminating contributions from v, and v3. 
Significant discrepancies arise, however, when vs, is compared with 
(7? grad p/8,). It is immediately obvious that, whereas vg is practically indepen- 
dent of temperature in the range 1-2°K. to 1:5°K., 7, is known to vary rapidly over 
this temperature range. In ‘lable 2 values of 7, deduced from the slopes of the 
difference curves in Figure 9 are compared with the values measured in oscillating 
disc experiments (Andronikashvili 1948). The disagreement appears too large to 
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be explained away by experimental errors. For example, if one assumes that the 
discrepancy is entirely due to the constriction correction discussed in § 4, this 
correction can be evaluated by assuming that equation (19) has to apply at 1-2°x., 
and one then finds that if a similar correction is applied at 1-5°K. there is still a 
large discrepancy at this temperature. 

The possibility that the flow became turbulent should be carefully considered. 
Turbulence in the normal component alone can be completely discounted, as 
the Reynolds’ number v7p,/7, never rose above 750. Ifthe superfluid component 
is considered to be a liquid with zero viscosity, then turbulence is very likely 
to be set up in it, but if the fundamental equation describing its behaviour is 
taken to be curl v,=0, then it would be very unlikely to become turbulent. The 
question of turbulence in the superfluid component is therefore best left in 
abeyance until there is strong experimental evidence in its favour. 


§6. OSCILLATIONS 

The type of oscillation observed in very narrow channels (Allen and Misener 
1939) and in the film (Atkins 1950 a) was also observed in these wide capillaries. 
In the region of small pressure heads the frictional forces are small and vary 
approximately as v’, so that when v is less than about 1 cm.sec™! the frictional 
forces are not measurable as a pressure head. Consequently the equilibrium 
position is approached with a velocity of about 1 cm.sec~4, the inertia of the 
liquid is sufficient to produce an overshoot and subsequent oscillations, and the 
frictional forces are then so small that the damping of the oscillations is only 
slight. The period of the oscillations can be shown to be 


au 
tT) =20 (2 5) OMT ha OL ee (20) 


pry 
In the present case there are two complications arising because (a) the 
thermo-mechanical effect cannot be ignored, and (4) the film is in parallel with 
the capillary and alters the period by a few per cent. In Table 3 the observed 
periods are compared with the periods calculated from equation (20). Results 


Table 3 
Capillary I IV V VI 
Mean int. diam. (cm. x 10~?) 4-40) 0-262 3:29 2-10 
Capillary length (cm.) 48-6 7°76 76°8 93-6 
Cross-sectional area of reservoir (cm?) 3°50 0-00571 3250) 0-704 
Form of capillary Zigzag Straight Spiral Spiral 
Observed 12 Se 60-0 L557, 12052 92°8 
period (sec.) \.1°5° k. 64:3 19-2 — 100-2 
Calculated f{1:2°x. 59-5 18-4 99-7 86:3 
period (sec.) |. 1:5° K. 62-9 19-6 —= 92-0 


for the spirals are included because end effects and constriction effects are not 
likely to be important at the small velocities involved in the oscillations. Both 
the thermo-mechanical effect and the film effect would decrease the period so 
that where the observed period is greater than the calculated period, as in some 
of the spiralled capillaries, this suggests that the mean radius was slightly less 
than the measured value because of the distortion produced when the capillary 
was bent. Equation (20) is based on the assumption that the velocity is constant 
over the cross section of the capillary, and, if this were not so, 7 would have to be 
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replaced by an effective radius. The approximate agreement between observed 
and calculated periods excludes such possibilities as that the flow is mainly 
confined to the walls of the capillary. In a more accurate experiment specifically 
designed to investigate oscillations and to eliminate the disturbing effects 
mentioned above, it might be possible to decide whether the velocity really is 
constant over the cross section of the capillary. 


§7. CONCLUSIONS 

In capillaries of diameter greater than 10-? cm. the critical velocity, if it exists, 
is less than 1cm.sec~!, but there is evidence that it is rather greater than this 
in a capillary of diameter 0:26 x10-?cm. The major part of the flow in these 
wide capillaries is opposed by frictional forces and is not truly ‘superfluid’. 
The Gorter—Mellink theory of mutual friction is almost adequate to explain the 
results, but there are discrepancies, and it seems that the results are explicable 
only if there is some other type of frictional force in addition to the mutual friction. 
This conclusion is considerably strengthened by the oscillating disc experiments 
of Hollis-Hallett (1950), and will be discussed in further detail in a later paper. 
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X-Ray Asterisms from Deformed Crystals 


It has been suggested (Barrett 1940) that x-ray asterisms from deformed metal crystals 
originate from deformation bands, i.e. narrow regions wherein variations in lattice orientation 
occur. Recent work with aluminium has shown that the term deformation band embraces 
at least two different phenomena (Honeycombe 1951). A distinction has been made between 
bands of lattice curvature which have been called kink bands and bands in which local 
orientation differences have arisen because of slip on a second set of planes. The latter 
have been called bands of secondary slip. Evidence obtained by optical (Cahn 1950, 1951) 
and x-ray micrography (Honeycombe 1950) has suggested that the kink bands are responsible 
for the major part of the asterism smear in x-ray Laue photographs of aluminium single 
crystals deformed in tension. 
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Ordinary x-ray methods, using beams of 0:5 to 1mm. diameter, cannot give direct 
evidence since the kink bands occur on a microscopic scale (Figure 1). They are usually 
about 0-01 mm. wide and between 0-02 and 0:2 mm. apart, depending on such variables 
as the extent of deformation, the orientation of the crystal and the purity of the metal. 
It appeared, however, that application of the x-ray microbeam technique (Hirsch and 
Kellar 1951, Gay et al. 1951) would provide direct evidence on this problem. The specimen 
examined was aluminium of 99-99°% purity and of very large grain size (1-2 cm.). It was 
electropolished in the annealed condition, then strained in tension by about 7:5% elongation. 
One grain which had developed widely spaced kink bands was selected for detailed exami- 
nation using an x-ray microbeam of 100 diameter. ‘To enable accurate location of the 
beam on chosen microscopic areas, a special specimen holder was used, details of which 
will be published elsewhere. 

It was found that if the beam was directed on a region between two kink bands, the 
resulting Laue spots showed little asterism (Figure 2). If the width of the kink band was 
completely covered by the incident beam, then marked asterisms appeared in the x-ray 


Figure 2. Section of amicrobeam Laue 
photograph | taken with the x-ray 
beam incident between the kink 
bands. 


Figure 3. A similar section of a micro- 


Figure 1. Micrograph of surface of beam Laue photograph taken with 
electropolished deformed specimen the x-ray beam straddling a kink 
of aluminium, showing kink bands band. The reflections are the 
(magnification 75). same as those shown in Figure 2. 


photographs (Figure 3). For the cases examined, the range of orientation within the 
band was approximately three degrees. No marked structure was observed in the asterism 
streaks. In previous macrobeam experiments on lightly deformed aluminium crystals, 
the asterisms usually contained intensity maxima. These have been attributed to 
reflections from the slightly disorientated lamellae between the kink bands, the latter 
being responsible for the diffuse part of the asterism. This view is confirmed by the fact 
that in the present experiments, in which only one lamella is irradiated, the asterism streaks 
are continuous. 

It is recognized that asterisms also arise from macroscopically bent crystals, but these 
experiments show that for aluminium crystals deformed by small amounts in tension the 
asterism is mainly caused by the kink bands, which are regions of curvature on a microscopic 
scale. It is relevant to note that crystals of some other metals, e.g. cadmium, which do not 
develop kink bands when deformed in tension, show very little asterism. 


Cavendish Laboratory, P. Gay. 
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On the Identity of B-Rays with Electrons 


Goldhaber and Goldhaber (1948) have pointed out that if the identity of negative /-rays 
with atomic electrons were not complete, then, in accordance with the Pauli principle, f-rays 
might be expected to undergo capture into bound orbits which are already filled with atomic 
electrons. Such a process must occur as a prelude to the capture of negative mesons by 
nuclei (Fermi, Teller and Weisskopf 1947). If this effect exists one should detect character- 
istic X-radiation from the capturing atom, modified by screening and the possible differences 
in mass or spin between the $-particle and an atomic electron. Goldhaber and Goldhaber 
(1948) gave an experimental upper limit of 3°% to the number of f-rays from a 4C source 
which could have been captured into the K shell of a lead absorber. 

A more suitable source and more efficient detection have allowed us to improve on this 
figure. 

We have examined the bremsstrahlung from the f-rays emitted by a source ot 9H 
absorbed in a thin layer of zirconium (approximately 0:2 mm. thick) in which sensibly 
all the B-rays were stopped. A Nal (TI) crystal (Harshaw Chemical Co.) and a photo- 
multiplier (E.M.I. type VX5032) were used to analyse the x-radiation in the manner des- 
cribed by West, Meyherhof and Hofstadter (1951). 
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Pulse height distributions for bremsstrahlung due to tritium f-rays stopped in zirconium (curve A) 
and L x-rays (mean energy 11°8 kev.) of bismuth (curve B). 


In the Figure is shown the pulse height distribution from the tritium source; on the 
same Figure is plotted, as an indication of the resolution and the calibration, the distribution 
for the bismuth L x-rays. 

Tests showed that over the energy range 10-100 kev. the pulse height varied linearly 
with quantum energy. 

The general form of the spectrum can be explained as bremsstrahlung from the tritium 
B-rays modified by the spread in pulse height and the heavy absorption in the Zr below 
about 6 kev. (corresponding to 12 volt pulses). 

From our results we may set an upper limit to the number of f-rays of 7H which could 
have been captured into a K shell of Zr. Since the mass of a f-particle of 3H can differ 
from that of an electron by at most about 4°, we should expect to find the radiation from 
its capture into a K orbit of Zr close to the normal K line, i.e. at approximately 16 kev. 

The pulse height distribution fails to reveal any trace of an x-ray line in this region; 
we may set an upper limit to the intensity of such a line by attributing to it all pulses of 
height exceeding 23 volts (corresponding to hv=15-8 kev.). The resulting counting rate 
is 10* per second or, when corrected for self absorption, 2:5 x 104 per sec. Since the crystal 
was sufficiently thick to stop all the x-rays reaching it from the source we have assumed 
that the detection efficiency was close to 100%. 

The number of f-rays per second from the source of 3H was estimated as 10! from the 
known volume of tritium occluded in the zirconium. Hence not more than 2°5 in 10® 
f-rays of tritium can have been captured into a K shell of zirconium. 
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In recent notes Charpak and Suzor (1950, 1951) explain certain details of their experi- 
ments on the back scattering of positive and negative f-rays (as compared with internai 
conversion electrons) by the assumption that f-rays are emitted as charged particles differing 
fromelectrons, into which they are subsequently transformed bysome unspecified mechanism. 
Irrespective of whether this suggestion can be substantiated or disproved, it raises the 
question as to how soon after the instant of disintegration a f-particle may be identified 
with an atomic electron. 

For the f-rays of tritium we can set a limit to this interval if we can estimate the time 
for their slowing down and supposed capture into a K orbit in zirconium. ‘The time to 
reduce them to thermal velocities will be of the order 10-1 sec.; the time for capture of 
a meson of electronic mass into a K orbit has not been calculated, but if we assume that 
it is of the order of the lifetime of an x-ray transition (i.e. 107! sec. in zirconium) then the 
predominant factor will be the slowing-down time. 

Thus we may conclude that in the case of the f-rays of 9H not more than a few in a 
million can be other than electrons at a time of the order 10-™ sec. after disintegration. 


The Clarendon Laboratory, W. T. Davies. 
Oxford. M. A. GRACE. 


31st May 1951. 
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Measurement of Fast Neutron Velocities 
by Delayed Coincidence 


In a nuclear reaction in which neutron emission is followed by a gamma-ray from an 
excited nucleus, a study of the angular correlation between these particles may enable one to 
infer properties of compound and residual nuclei. In such an experiment (Thirion and 
Muller 1951) there are in general many neutron groups, and the interpretation of results may 
be difficult. By means of a delayed-coincidence method we have found it possible to 
measure the velocity of a neutron that is followed by a gamma-ray in a nuclear reaction. 
Although no conclusive results have yet been obtained, we feel that our method is of sufficient 
interest to justify a short description. 

Photo-multipliers with liquid scintillators are used for detection of both neutron and 
gamma-ray. The multipliers are E.M.1. type end-window 14 stage tubes, with gains about 
10°. The scintillators consist of terphenyl in toluene (0:25 gm. per 100 cm?) and are con- 
tained in glass flasks coated on the outside with magnesium oxide. One phosphor (the 
‘gamma detector’) is one litre in volume, and situated close to the target. The ‘neutron 
detector’, five litres in volume, is placed one metre from the target. 

The principle is to measure the time of flight of the neutron to the large phosphor. 
This is done by using a fast coincidence arrangement and producing a delay between gamma 
counter and coincidence unit. Thus a2 Mev. neutron that travels one metre from the target 
arrives 5:1 x 10~® sec. later than the gamma-ray following it (corresponding to a length of 
10-2-metres of Uniradio 32 coaxial cable). In order to produce rectangular pulses of short 
duration, the output from each multiplier is fed to the grid of an E.M.I. type E.2133 
secondary-emitting pentode, whose twofold purpose is to transform the high output 
impedance of the multiplier and produce pulses of uniform height. The pulses are ‘ clipped’ 
by means of a shorted delay line (Bell and Petch 1949). A plot of coincidence rate against 
delay time using a ®°Co source showed a half-width of 2:2 x 107° sec., which for a 2 Mev. 
neutron corresponds to an energy spread from 1:4 to 3:2 Mev. 

As an example of what can be obtained under these conditions, the Figure shows the 
delay spectrum for coincidences from the *°Be(d, ny)!°B reaction at 0-48 Mev. bombarding 
energy. The large peak corresponds to gamma-gamma coincidences and to ‘reversed’ 
neutron—gamma coincidences (neutrons recorded by the ‘gamma-counter’ and vice versa). 


ema 
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The small peak corresponds to a mean neutron energy of 2:5 Mev. In order to detect the 
angular correlation between this neutron group and its associated gamma-ray the neutron 
counter was kept fixed and the gamma counter rotated with respect to it, in the plane 
normal to the deuteron beam. This experiment gave isotropy to within 3%. 

Since this work was started the neutron spectrum from this reaction has been studied in 
detail (Ajzenberg 1951), and it seems likely that the ‘ group’ is in fact a mixture of three 
groups, too close in energy to be resolved in our experiment. Hence the lack of correlation 


Coincidence Rate (arbitrary units) 
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Delay spectrum of coincidences from *Be (d, ny) ?°B. Arrows indicate delays corres- 
ponding to known neutron groups (Ajzenberg 1951). 


means little or nothing. In principle the energy resolution can be improved to any desired 
accuracy by moving the neutron counter further from the target, but a limit is set by the 
drop in counting rate; the background rate does not however decrease proportionately, 
since it is due partly to slow neutrons entering the counters from random directions. 

The method described can be used to measure velocities of other types of particles in 
coincidence with gamma-rays, and it is hoped to.apply it to a study of alpha-gamma coin- 
cidences resulting from nuclear reactions. 


Cavendish Laboratory, D. B. JAMEs. 
Cambridge. PB 
15th June 1951. 
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An Investigation of the Reaction °Be (d, p) 


Thin evaporated targets of beryllium have been bombarded with an unanalysed deuteron 
beam accelerated by an electrostatic generator. The deuteron content of the beam was 
determined, by magnetic analysis, to be (92+ 2)%. The particles emitted from the resultant 
reactions were detected in Ilford C2 nuclear emulsions. The voltmeter measuring the 
generator voltage was calibrated by means of the *Li(p, y) reaction, and was known to be 
linear in the voltage range used. 

In some of the experiments the emulsions were placed in a small magnetic spectrograph. 
The intensity of each group detected was plotted as a function of distance along the emulsion 
from a fixed and known reference line. The spectrograph was calibrated by means of the 
direct deuteron beam which was arranged to pass into the spectrograph at known and 
constant energies. 

Although the proton content of the beam was small, sufficient protons were present to 
produce a particle group from the *Be(p, d)*Be reaction. By taking the position of this 

- group as the position of half-maximum intensity, on the high energy side of the group, the 
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Q-value for this reaction was found to be (560+ 13) kev., in agreement with previous values, 
see Hornyak et al. (1950). ‘The position of half-maximum intensity of each group observed 
was therefore taken empirically as the position of that group, following the method 
described by Buechner et al. (1949). 

The group emitted from the reaction *Be(d, p)!°Be* was observed and the Q-value of 
the reaction determined from the measured momentum of the group. This gave 
O=(1,214+ 13) kev. 

The Q-value determined from the range of the particles in the emulsion was 
OQ=1,204 kev.; this is the mean of three measurements each having a standard deviation of 
45 kev. These Q-values are in agreement with that of (1,201+7) kev. obtained by 
Buechner and Strait (1949). 

From a knowledge of the acceptance of the spectrograph, determined from the geometry, 
an estimate was made of the cross section of this reaction at two energies. The results are: 


o (90°) =(0-43 + 0:09) x 10-78 cm? sterad at 301 kev., 
o (90°) =(2:96 + 0°59) x 10-8 cm? sterad~! at’ 440 kev. 


The thickness of each target was calculated from the known weight of beryllium 
evaporated on to each backing. The error in the calculated thickness is estimated at + 20%. 
The pressure at the bottom of the accelerator tube and above the target was 10~° mm: Hg or 
less, the scattering of protons out of the true image was about 2%. The intensity of 
contamination products from the. 1*O(d, p)!7O* reaction was less than 3% of the 
intensity of the °Be(d, p)!®°Be* groups. This was determined by blank target exposures, 
and by evaporating emulsions in an auxiliary camera to all the reaction products,so that 
the group from the 1%O(d, p)!7O reaction also could be observed. The angle of aip of 
the tracks in the emulsion in the spectrograph was about 15°. 

The Q-value for the reaction °Be(d, p)!°Be was determined from the range of the 
particles in emulsions exposed in the auxiliary camera, which is similar to that used by 
Rubin (1947). For this reaction Q=(4:55+ 0-03) mev. ; 

The energy measurements were based on the range-energy relationship due to Rotblat 
(1950). Inallthe measurements made this relationship gave better agreement with Q-values 
than did the relationship due to Lattes and others. 

I am indebted to Mr. D. R. Petrie who designed the Bpettomrapl: 


Research Laboratory, A. J. SALMON. 
Associated Electrical Industries Ltd., 
Aldermaston, Berks., 
2nd July 1951. 
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Note on Experimental Determination of the Half-Life of ‘Be 


«O(y, 4a) photodisintegration stars frequently involve a ground-state ®Be nucleus as 
an intermediary product (Goward and Wilkins 1950). In recent studies of such stars 
with Ilford E1 nuclear emulsions (Millar and Cameron 1951), the observation of single 
tracks produced by *Be nuclei, prior to splitting into two alpha-particles, has been reported 
and a half-life of (5°34 1:1) x 10-1 sec. deduced. This long half-life has been advanced 
(Stearns and McDaniell 1951) as evidence for J=2 for the ground-state of *Be. The 
observations in the present note on 280 similarly formed ®Be nuclei do not support the 
possibility of such a definite deduction of the half-life. 

Examination of the narrow V-shaped alpha-particle tracks from the above-mentioned 
280 1O(y, 4x) stars leads to the following conclusions. Firstly, closely neighbouring 
tracks of this kind can only be resolved sufficiently to allow an upper limit (not an actual 
value) of about 5 x 10~™ sec. to be obtained for the half-life of Be. Secondly, the degree 


850 Letters to. the Editor 


of scattering evident in the tracks is sufficient to produce an occasional illusion of a single 
*Be track. The importance of these factors may be estimated from photographs already 


published (Wilkins and Goward 1950). Processing techniques directed especially towards _ 


studying the ®Be half-life may give a limited improvement in resolution, but the use of 
grain gradation development is more likely to increase the illusion of a single *Be track 
by thickening the unresolved part of the V-tracks. 

The above qualitative remarks may be supported by simple numerical considerations, 
complications such as scattering being neglected. Assume that a separation of about 
0-3 micron (the grain diameter) is sufficient to give clear resolution between tracks, 
whatever their orientation relative to the observer, and further that the *Be nuclei are 
ejected with 2-5 Mev. kinetic energy and disintegrate immediately and isotropically. Then 
about 35% of the resulting V-tracks have the first micron or more unresolved, and an upper 
limit of 9x10-!sec. might be deduced for the half-life. Extrapolation of the tracks 
backwards would reduce this upper limit estimate, but hardly by a factor of more than two, 

The resolution of the method might be improved with considerably higher *Be kinetic 
energies, although the advantage of a higher initial velocity is almost exactly offset by a 
narrower V angle. It is significant, however, that a report on some observations of *Be nuclei 
in cosmic-ray stars by Crussard (1950) makes no mention of single *Be tracks in 16 events 
(see also Perkins 1950). A specific attempt by Hodgson (private communication) to extend 
the present results by measurement on such stars places an upper limit of about 
2 x 10—*4 sec. on the 8Be half-life. 


Atomic Energy Research Establishment, J. J. WILKINS. 
Harwell, Didcot, Berks. . F. K. Gowarb. 
2nd July 1951. 
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On the Supplementary Condition in Quantum Electrodynamics 


It has been shown by Koba, Tati and Tomonaga (1947) and Schwinger (1948) that 
when we pass over from the Heisenberg representation to the interaction representation, 
the supplementary condition for the electromagnetic field has to be modified by adding 
a ‘charge-term’. Since the treatment by these authors is fairly complicated, it is not 
easy to see why the additional charge-term should appear in the interaction representation. 
It will, therefore, be of interest to show that the supplementary condition in the interaction 
representation may be obtained in a very straightforward way by the general method of 
passing over from the Heisenberg representation to the interaction representation. The 
treatment, to be described here, also proves very convenient for obtaining the supplementary 
conditions (or the coordinate conditions) for Einstein’s gravitational field in the interaction 
representation (Gupta, to be published). 

Fermi’s supplementary condition. According to Fermi the supplementary condition 
for the electromagnetic field in the Heisenberg representation is given by 

OAM 
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where we use primed quantities for the Heisenberg representation. The condition (1), 
which is required to hold at all times, is equivalent to the initial conditions 
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In order to pass over to the interaction representation, we can use the following general 
method for any field quantity. First, remove all explicit time derivatives from the field 
quantity in the Heisenberg representation by expressing it in terms of the field variables 
and their canonical conjugates. Then, replace these field variables and their canonical 
conjugates by the corresponding variables for free fields. This gives the required field 
quantity in the interaction representation. 

If we use the Lagrangian density, given by Schwinger (1948), for a system of photons 
and electrons, the canonical conjugates of A, will be given by 
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Using (3), and the field equation []?4,’=—j,’, we get 
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Therefore, corresponding to (2), we have in the interaction representation 
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The above equations may be w-itten in a moe elegant form as 
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the equivalence of (6) and (7) being obvious from the fact that 
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It is interesting to note that it is the second relation in (2), and not the first, which gives 
tise to the charge-term. The initial conditions (7) are, of course, equivalent to requiring 
that the supplementary condition 

OGG) SOsF Ses eee Ay: (9) 
should hold at all times. 

Modified supplementary condition. Fermi’s supplementary condition, discussed above, 
involves the difficulty that there is no normalized state which satisfies it (Ma 1949, 
Belinfante 1949). Therefore, the present writer (Gupta 1950) has given a new treatment 
for the longitudinal part of the radiation field, in which an indefinite metric is used for the 
scalar photons and (1) is replaced by a modified supplementary condition 
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[eA,//ex,]* being the positive frequency part of @A,’/éx,. In the absence of 
interaction (10) may also be written as (04,/‘+)/dx,)¥’=0, where A,/‘~ is the positive 
frequency part of A,’. But, in general one has to use the form (10), because in the presence 
of interaction A,’ cannot be split into the positive and the negative frequency parts. 
However, such a splitting is still possible for 0A,,’/@x,, which always satisfies the wave 
equation for plane waves: []? (0A,’/@x,) = — 0j,//0x,=0. 

If we now pass over to the interaction representation, then, as Bleuler (1950) has already 
done, we ought to replace (9) by 
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where Q“*) is the positive frequency part of Q, which satisfies the wave equation [}?Q=0. 
Thus, we have from (8) : 
34, 4x) : 
BiG) a eee DOG =e eGo) te (12) 
OX ti=t 
A,,*) and D(x) being the positive frequency parts of A, and D(x). 


Cavendish Laboratory, Suraj N. Gupta.* 
Cambridge. 
3rd July 1951. 
* State Scholar of the Government of India. 
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The Band Systems of Thorium Oxide and Hafnium Oxide 


In continuation of the work in this laboratory on the band systems of titanium oxide 
the systems of thorium oxide and hafnium oxide have been investigated. Both spectra 
were developed in arcs between graphite rods, the lower positive of which was filled with 
the oxide in question. Glass and quartz Littrow spectrographs were used to photograph 
the hafnium oxide spectrum, while a 21-ft. grating could also be used for the less expensive 
thorium oxide from A 8500 to A 3800. 

The spectra were very rich in red degraded systems throughout the whole region 
photographed and a number of sequences could be recognized. Extrapolation of the 
vibration frequencies of the other elements of the group of the periodic table suggested the 
corresponding frequencies in thorium oxide and hafnium oxide as about 740 cm~! and 
865 cm™! respectively for the ground states. 

About fifteen systems in ThO and eight in HfO could be distinguished. The energy 
level diagrams of the two molecules showed three important lower levels at values of v: 0, 
2721 and 4177 for ThO and 0, 2084 and (4003) for HfO. These levels are probably *IT 
levels, and the ®IJ,—#II, intervals 2721 and 2084 of ThO and HfO are comparable with 
the °F,—*F; intervals 2869 and 2320 of the arc spectra of thorium and hafnium. It is 
hoped that full details will be published shortly. 


Spectroscopy Laboratory, S. G. KrISHNAMURTY. 
Imperial College, London. 
12th July 1951. 


REVIEWS OF BOOKS 


Weltsystem, Weltather und die Relativitatstheorie, by K. JELLINEK. Pp. xv +450. 
(Basle: Wepf, 1949.) Sfr. 45: 


This is a book the like of which one would not, during the past thirty or forty years, have 
expected ever again to see published: for it is no less than an attempt to retain a universal 
aether and a universal absolute frame of reference while seeking to reconcile with these many 
of the notions and results of relativity theory. 

One must hasten to admit that, in view of recent trends in cosmological theory, the attempt 
must not be condemned out of hand as being entirely misdirected, as it almost certainly 
would have been not many years ago. It is in fact well known that almost all current theories 
of the expanding universe introduce what is usually called ‘ cosmic time’, which may seem 
very like the universal time of pre-relativity physics. This feature has been particularly 
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emphasized during the past couple of years in the writings of Bondi and Gold, and of Hoyle, 
on the possibility of continuous creation; for if newly created matter does appear throughout 
the universe it must apparently serve to define an absolute standard of rest at every point of 
space (i.e. a preferred time-direction at every event in space-time). 

Nevertheless, science is never known to go back completely onits tracks. Ifitreverts toa 
position something like one it has held prior to what has been regarded as a major advance, 
the likeness is found to be superficial or even illusory. The photons of modern physics are 
not really much like Newton’s light corpuscles. Similarly, any preferential frame of reference 
in modern cosmology is conceptually far removed from a frame at rest in a pre-relativity 
aether. ' 

As far as I can understand him, Professor Jellinek wants to keep the nineteenth century 
aether and to make the best job he can of accommodating relativity in it. The job is a 
difficult one, and Professor Jellinek’s intellectual candour is such that he conceals none of the 
difficulties. Indeed, reading his book one repeatedly expects him to produce the standard 
arguments which originally showed how much simpler it is to do away with the aether and to 
accept the orthodox formulation of relativity theory. But he never surrenders. 

Professor Jellinek’s purpose is to provide a serious treatise of not too mathematical a 
character on special and general relativity for the use of experimental physicists and others 
who are not professional mathematicians. He is justified in his belief that no such book 
exists which includes general relativity. The reason is that hitherto the latter has been of 
very little interest to experimental physicists as such. Quite recently it has gained some 
interest for them on account of the possible connections between cosmological theory and the 
theory of the relative abundances of the elements. Owing to Professor Jellinek’s unusual 
presentation, and the extremely involved discussions of particular points to which it leads 
him, his book, unfortunately, seems scarcely to meet the need. 

Here is not the place to attempt a lengthy criticism of Professor Jellinek’s treatment. It 
would be largely a repetition of arguments that were familiar in the early days of relativity 
theory. But it would also demand the admission that much of what he says in criticism of 
other treatments is in itself very pertinent, though I think it can be countered in more simple 
ways than he allows. ‘Two general points may be made. First, a great merit of general 
relativity theory is that it is a ‘differential ’ theory, i.e. it can be formulated in terms of the 
properties of space-time in an observer’s neighbourhood. Even if it ultimately leads to the 
recognition of a preferential universal frame of reference in a cosmological model, it does not, 
as Professor Jellinek seems to assert, require to assume the existence of such a frame at the 
outset. The second point is that it is primarily only in the theoretical models of a ‘smoothed- 
out’ universe that general relativity recognizes a preferential frame of reference. If we 
liken such models to simple surfaces of revolution, as we can, the preferred frame of reference 
corresponds to the meridian sections and circular sections, and it is obvious that these are to 
be recognized as a preferred system. Buton this analogy we should liken the actual universe 
to the surface of, say, a very rough tree-trunk. Though we might again find a preferred 
coordinate system suitable for describing its general features, such a system would play no 
special part in the description of minute local configurations in the surface. Similarly, I 
think that concepts such as ‘ cosmic time ’ have very little to do with applications of relativity 
theory to anything but the very large-scale features of the universe. These concepts do not 
seem to me to offer any justification for thinking that relativity theory requires anything that 
corresponds to a universal aether in Professor Jellinek’s sense. 

Anyone who does want an ‘ aether’ in relativity theory can, however, have it for the 
asking. ‘The theory deals exclusively with the entity it calls ‘space-time’. Where space— 
time possesses such a curvature that the Einstein tensor is non-zero, it asserts that ‘matter’ is 
present. Where the curvature is such that this tensor vanishes, it says that space-time is 
‘empty’. That is, it does not regard matter as something immersed in space-time, but 
as being a property of space-time. Where space-time does not possess this particular 
property, general relativity still treats it as being by no means featureless. Therefore, to 
give it a name like the ‘ aether’ might well be more appropriate than to call it ‘empty’ 
space. This, though not in the context of orthodox general relativity, seems to be the point 
made by Eddington (Fundamental Theory, 1946, p.172). But to retain the term ‘ aether ’ in 
this way would be to run the risk of retaining also those properties of the classical aether 
which relativity has abandoned. It seems to the reviewer that these are just the properties 
that Professor Jellinek does desire to retain. 


854 Reviews of Books 


Nevertheless, it is abundantly evident that Professor Jellinek has pondered long and 
deeply upon the subjects of which he writes, and what he has written commands one’s 
respectful attention. If it cannot be deemed wholly suitable for the requirements of the 
experimental physicist making his first acquaintance with relativity theory, it will probably 
be found of much value to other relativists who may themselves be called upon to supply 
these requirements. It will stimulate them to think of their subject from fresh points of 
view and will be achallenge to them to re-examine the foundations of their own presentations. 

W. H. MCCREA. 


Sourcebook on Atomic Energy, by S. GLassTONE. Pp. 546. (London: Mac- 
millan, 1950.) 24s. 


The book under review is of the type perhaps more common in the United States than in 
Great Britain. It was apparently originally requested as a comprehensive sourcebook on 
atomic energy for the use of textbook authors and editors. It stands somewhere between 
an elementary textbook and an encyclopaedia. As stated in the foreword, Dr. Glasstone has 
written a book which is of far broader usefulness and which provides a source of basic 
atomic energy information for everyone interested in this field. 

Beginning with the earliest theories of the atom and its structure, it describes the growth 
of thought and knowledge in the field—the development of the theories of the nature of 
electricity, energy and matter, the discovery of the phenomena of radioactivity—and so leads 
to problems of modern developments. 'The following chapters deal with the acceleration of 
charged particles, nuclear transmutation and artificial radioactivity, the neutron, nuclear — 
structure and nuclear forces, nuclear fission, utilization of nuclear energy, new elements, the 
uses of isotopes, and finally cosmic rays and mesons. It concludes with a chapter giving an 
outline of radiation protection and ‘ health physics ’. 

It is certainly an amazing concentration of information which, in view of its nature, is 
surprisingly readable, though it is certainly not a book to be read through at one sitting. 
The chapters are divided into a large number of paragraphs each of which is clearly labelled, 
thus making cross reference much easier and more satisfactory. Each of the subjects is 
dealt with largely historically. Only a rudimentary knowledge of science or mathematics is 
required to follow the narrative, which will doubtless be of interest to a very large number of 
readers without much scientific training. 

Atomic energy is of great importance in medicine, biology, astronomy, physics, and other 
branches of knowledge. This miniature encyclopaedia, arranged nevertheless on a logical 
system, will probably be of great interest and use to the type of reader for whom it was 
intended. It will be the basis of a good many popular articles in newspapers even if not, as it 
is surely not intended to be, useful as a basis for further advances. It is a well written and 
a skilfully organized account of the fundamentals of ‘ atomic energy’ and its applications. 

. W. V. MAYNEORD. 


Introductory Nuclear Physics, by Davip Hatiipay. Pp. 558. (New York: 
Wiley and Sons; London: Chapman and Hall, 1950.) 52s. 


Professor Halliday has written a remarkable and interesting book which is intended for 
post-graduate students about to start experimental work in nuclear physics. It is a typically 
American book with plenty of facts well displayed in convenient tables, plenty of references 
to modern work and plenty of diagrams of apparatus which has been in use in the last five 
years. ‘The student is early introduced to the idea that one should consult the literature, 
and above all reminded that the best way to find out whether you understand the meaning 
of an expression is to work out a few elementary numerical examples. These are amply 
provided at the end of each chapter ; on the whole they are well graded and match the easy 
descriptive standard of the text better than is usual in American books, in which a hard 
text is so often combined with trivial examples. 

Clearly an author who sets out to cover at an advanced level the whole of nuclear 
physics in 500 pages has to jettison much traditional material. 

History is the first to go. The opening sentence baldly states that ‘‘ Nuclear physics 
began in 1896 when Henri Becquerel discovered radioactivity’, and leaves it at that. 
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Its further progress is summarized in a table giving the author, date and country of origin 
of the chief subsequent discoveries. The only ghost from the past which is allowed on 
the stage is the familiar diagram of Geiger and Marsden’s apparatus for measuring 
Rutherford scattering. It is, however, rapidly exorcized by a footnote which remarks 
that “‘ with modern techniques this pioneer experiment can be done as a three-hour student 
exercise’. I wonder how many students do it on this side of the Atlantic. 

Secondly, no attempt is made to derive, where possible, well known and important 
results by relatively elementary methods. Although a fairly extensive general knowledge 
of wave mechanics is assumed, the book contains practically no mathematics. 

The author is anardent, logical positivist, and one gains the impression that he thoroughly 
enjoyed throwing overboard the dusty treasure of the past. The merits of this course of 
action, when pursued with Professor Halliday’s vigour and enthusiasm, depend very 
much on the graduate training of the student. The differences in this respect between 
English and American students may somewhat reduce the value of the book to English 
students, who on the whole are probably better grounded in classical physics and less 
immediately aware of the elementary properties of neutrons and high energy particles. 
For example, the book has passages in small print judged too hard for inclusion in a graduate 
course; in what British University would one find students to whom the statement 
that “‘ Mott and Massey give the differential cross section for neutron—proton scattering as 


do=(A?/167”) | (21+ 1)P; cos 0(e2#81—1) |? dQ” 


would be reasonably intelligible and for whom the proof of the formula AE/E=4mM, 
for the head-on collision of elastic spheres, is considered so hard as to be relegated to small 
print? This is only one of many similar anomalies. 

Nevertheless the amount of new material included is stimulating and surprising, a 
there can be few who will not find in these pages some aspect of nuclear physics of which 
they were only dimly aware. 

On the whole the author’s philosophical position is not unduly in evidence, otehs in 
some instances it has led him to omit material which should have been included. | For 
example, the discussion of the specific ionization of charged particles omits much important 
material. The index has no entry under Bragg (in this connection), although the term 
“ Bragg curve ’ is still in common use in American journals; no mention is made of the 
different mechanisms of ionization which account for the falling and rising parts of the 
curve. Furthermore, no mention is made of the convenient and empirical Geiger law, 
although the standard measurements in this field are still tabulated as deviations from this 
law. 

In spite of these minor criticisms the book can be confidently recommended. ‘The 
research student provided with this vade-mecum will find that he can rapidly run down a 
list of references and advanced texts dealing with the problems facing him in the laboratory 
on starting work in nuclear physics. What use he makes of them is up to him. Well, 
so it is; life is like that ! 

I imagine that lecturers who have to give an advanced course in nuclear physics will 
find this volume most useful because it supplies just what is needed in the way of references 
and tabulated data to form the coping-stone of a more traditional treatment. 

The standard of paper, printing and diagrams is excellent, and the price no more than 
holders of a devalued currency can expect. Gp lslg C 


Superfluids, Vol. I: Macroscopic Theory of Superconductivity, by F. LoNpoN. 
Pp. viiit+ 161. (New York: Wiley and Sons; London: Chapman and Hall, 
1950.) 40s. 


In 1937 Professor F. London expounded, in Une Conception Nouvelle de la Supra- 
conductibilité, the views of the nature of the superconducting state which he and his brother 
had introduced some two years before. Now, in the first volume of a projected series under 
the general title Super fluids, he sets out the same views once more, in a considerably expanded 
form. ‘The fact that in fifteen years he has not seen fit to make any substantial alterations 
to the fundamental framework of his theory is to be regarded not as evidence of a conser- 
vative nature, but as testimony to the logical beauty of the theory in its original form. It 
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must be admitted, however, that there would in any case be strong justification for a passively 
conservative attitude, since a rigorous experimental test of the theory has so far not been 
undertaken on account of the technical obstacles to such an investigation. 


The book is divided into five sections. The first, a brief development of the concept of 


the ‘ ideal superconductor ’, with an account of its thermodynamical properties, is followed 
by a long discussion of the electrodynamics of superconductivity. This section, which 
occupies nearly half the book, is perhaps the most satisfactory, and is a real pleasure to read. 
The basic idea, that the supercurrent is linked to the magnetic rather than to the electric 
field, is introduced as a postulate in the form of the well-known ‘ London equations ’, and 
the implications of these equations are then examined in detail. It is proved very elegantly 
that they lead to a unique field distribution if the sources and ‘ fluxoids ’ are prescribed and 
to perfect diamagnetism of macroscopic superconductors of any shape. Among other 
passages of great value may be mentioned especially the lucid discussion of the difference 
between perfect conductivity and superconductivity. 

It is perhaps because the disturbing shadow of experimental facts has not fallen upon this. 
section that it seems more complete than the two which follow. Here the intermediate 
state and the properties of thin films are discussed, and it becomes only too plain that the 
ideal superconductor is too remote a conception. Muchas we may admire the mathematical 
structure which Professor London has erected, we feel that the security of the experimental 


foundations is still open to question, and it is quite certain that the last word on these 


subjects is still to be said. 

Finally Professor London examines the relation between the macroscopic electro- 
dynamics and the quantum mechanics of an assembly of electrons, in order to formulate 
the conditions under which, in his view, the superconducting state may exist. The major 
requirement, that the wave functions of the electrons shall be rigid, in the sense that they are 
not readily perturbed by an applied magnetic field, is one which has not in the past been 
accepted by all who have attempted to construct theories of superconductivity, but it is 
significant that the most recent attempts have been along these lines, and show in other 
respects a considerable improvement on. their forerunners. Though short, this section 
contains much weighty matter, and is certainly one of the most important. 

The general impression which the book leaves is that it is not a balanced textbook on 
superconductivity—it is not even a balanced textbook on the theory of superconductivity ;. 
rather is it a personal account of an important aspect of the theory written by one who has 
contributed perhaps more than anyone else to our understanding of the phenomenon. 
It can hardly be recommended as an introduction to the subject for the general physicist, for 
whom a fuller account of the experimental facts would be desirable, but to the theoretical 
physicist, and above all to the low-temperature physicist, it represents an extremely valuable 
addition to the small body of literature in this field. A. B. PIPPARD. 


Process Heat Transfer, by D.Q. Kern. Pp. xii+871. 1st Edition. (New York 
and London: McGraw-Hill, 1950.) 68s. 


The justification for an addition to the already extensive literature on heat transfer 
must lie in a new or novel approach to the subject, and in this respect Professor D. Q. Kern’s 
book has serious claims to consideration. Whilst we do not consider his aim of providing 
fundamental instruction in heat transfer “‘ using the methods and language of industry ”’ 
as being fully capable of reelization or even as being altogether desirable, the endeavour 
to adhere to it has resulted in a book containing much practical data of direct value to the 
designer. An admirable feature of the book is the linking up of theory and practice by 
means of well chosen practical examples worked out in full and illustrated by clear line 
drawings of plant and equipment. 

According to the author’s definition, Process Heat Transfer deals with the rates of heat 
exchange as they occur in the heat transfer equipment of the engineering and chemical 
processes. A typical problem of process heat transfer is one concerned with the quantities 
of heat to be transferred, the rates at which they may be transfered because of the nature 
of the bodies, the driving potential, the extent and arrangement of the surface separating the 
source and receiver, and the amount of mechanical energy which may be expended to 
facilitate the transfer. 
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The first five chapters of the book deal on orthodox lines with the mechanism of heat 
transfer by conduction, convection and radiation, with fundamental definitions of heat 
transfer coefficients and with the relation between parallel flow and counter flow. 

The main body of the book then treats in detail the more important classes of heat- 
transfer equipment, including all common commercial types of exchangers, heaters, coolers, 
condensers, evaporators, reboilers, and distillationcolumns. In view of recent developments. 
the reader will find a useful summary of the methods of deriving the efficiencies of various 
forms of extended surfaces such as finned tubes. 

Other sections of interest relate to unsteady-state processes as exemplified by the 
regenerator and the reversing exchanger which is now coming into extensive use in gas 
separation plant, and to ‘furnace calculations based upon empirical and semi-empirical 
methods. The final chapter of the book deals with the control of process variables. 

The book is well illustrated and, although the bibliographies are by no means complete, 
adequate references are given to most of the significant work of the past twenty years. 

D. M. N. 
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ABSTRACTS FOR SECTION B 


The Mechanical Properties of Metals, by N. F. Mott. (35th Guthrie Lecture.) 
ABSTRACT. A description is given of edge and screw dislocations in a close-packed 


cubic lattice, which leads up to an account of the theory of Frank and Read of the origin of 


slip bands. Consideration is then given to a crystalline grain containing Frank—Read 
sources; it is shown that one has to expect cross-slip and the formation of deformation 
bands. It is shown moreover that the formation of vacancies by moving dislocations plays. 
an essential role in these processes.’ Some speculations on the stability of deformation bands 
are given. It is suggested that the movement of vacancies (self-diffusion) plays an essential 
part in polygonization, recovery and steady-state creep, enabling dislocations in deformation 
bands to move out of their slip planes and so relieve stress. It is probable moreover that in 
these bands the stresses are several hundred times greater than the applied stress, and that 
this enables diffusion to occur at temperatures at which it would not be possible otherwise.. 
Finally an account is given of low temperature creep of the type which does not involve 
recovery. 


The Deformation and Ageing of Mild Steel: II—Characteristics of the Liiders 
Deformation, by E. O. Hatt. 


ABSTRACT. In this paper, some of the factors which influence the appearance of 
Liiders bands in mild steel are studied. It is shown that the Liiders band is adequately 
described by a uniform shear front, spreading over the specimen. In coarse-grained 
specimens, experiments indicate that this front becomes diffuse; diffuse bands are also 
present in strain-aged material, but here the diffuse fronts become sharper as the ageing 
becomes progressively longer. 


The Deformation and Ageing of Mild Steel: II1I—Discussion of Results, by E. O- 
HALL. 


ABSTRACT. Anattempt is made here to explain the observed phenomena in the yielding 
and ageing of mild steel, described in two previous papers, in the general terms of a grain- 
boundary theory. On this hypothesis, a satisfactory explanation of the variation of the 
lower yield point with grain size may be developed. It is shown that strain-ageing must 
involve two processes: a healing of the grain-boundary films, coupled with a hardening in 
the grains themselves. A discussion of the possible nature of the grain-boundary film is. 
also undertaken. 


Reversible Effects in the Magnetization of Nickel, by R. S. TEBBLE, W. D. CorNER 
and J. E. Woop. 


ABSTRACT. An account is given of investigations on the temperature dependence of the 
contribution from reversible processes to the magnetization of annealed and strained nickel 
wires. This involves the measurement of reversible susceptibility «, over a range of field, 
and of temperature (—50°c. to -+-200°c.), and a rapid automatic method of recording 
(x,, H) curves has been developed. The effect of temperature on x; is discussed in relation 
to Becker’s strain theory and reasons for the discrepancies are suggested. It is shown that 
the contribution of reversible processes to the total change in magnetization, integrated over 
the hysteresis cycle, increases with temperature and strain, from 15°% for an annealed 
specimen at —50° c. to 85% for a strained specimen (17-1 kg. mm~”) at 150° c. The effect 
is not uniform over the whole cycle and in the region of the coercive field a considerable part 
of the change in magnetization must be attributed to irreversible processes. 
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The Surface Tension of Supercooled Phenyl Ether, by C. Dovp. 


ABSTRACT. Determinations of the surface tension of phenyl ether have been made, 
using Jaeger’s method, for the liquid in both the ordinary and the supercooled state. 
No abnormal change with temperature, of the type already found for viscosity and 


dielectric constant, has been detected for surface tension as the liquid enters the supercooled 
region. 


Freezing of Supercooled Water, by A. W. Brewer and H. P. PALmer. 


ABSTRACT. Rau has claimed that, by repeated freezing and thawing, the temperature 
to which a drop of water supercools before freezing can be lowered and, with special 
care, the freezing temperature can be reduced to —72° c. ‘This is claimed to be the temper- 
ature of spontaneous nucleation of ice. At that temperature the ice forms in rectangular 
crystals which also melt at —72° cc. This note describes a careful attempt to repeat Rau’s 
work. Like Cwilong, we have been unable to do so except by deliberate contamination, 
and it is presumed that Rau’s apparatus was faulty and his drops were contaminated 
by his cooling fluid. Information is given concerning the reproducibility of the freezing 
temperature of individual drops when repeatedly frozen and thawed. Rau’s rectangular 
crystals have been observed in strong alcohol solutions but no information concerning their 
nature has been obtained. 


Spontaneous Condensation of Water Vapour in Expansion Chamber Experiments, 
by B. J. Mason. 


ABSTRACT. When clean, saturated air is subjected to a sufficiently large adiabatic 
expansion, small water droplets arise spontaneously from the vapour. An equation is 
derived for the growth rate of a droplet in an atmosphere of given supersaturation and 
temperature. This is used in conjunction with the Becker—Déring theory of nucleation to 
calculate the supersaturation and the concentration of droplets at successive time intervals 
during a cloud-chamber expansion, allowance being made for warming of the air by liber- 
ation of the latent heat of condensation. For large expansion ratios, the supersaturation 
attains a maximum and the droplet concentration a sensibly constant value before the end of 
the expansion. ‘The computed value of this maximum droplet concentration is in acceptable 
agreement with the recent observations of Frey. The experimental results appear incom- 
patible with Tolman’s recent theory concerning the variation of surface tension with droplet 
radius. ‘The predicted rates of droplet growth agree well with recent measurements made 
by Hazen. 


Light Flux received by a Spectrograph from a Spatially Extended Refracting Source, 
by L. A. Woopwarp and J. H. B. Georce. 


ABSTRACT. ‘The geometrical-optical treatment of Nielsen, which referred to a spatially 
extended source of refractive index unity, is extended to sources of any refractive index. 
Provided that the direction of emergence of the light from the source is nearly normal, it is 
shown that a change of refractive index from 1 to n modifies the light flux by the factor 1/n?. 
This result is independent of whether or not a condensing lens is used between source and 


spectrograph. Its relevance in connection with observed intensities in Raman spectra is 
discussed. 


Optical Properties of Selenium, by J. J. Down. 


ABSTRACT. Measurements have been made of the refractive index of amorphous 
selenium in the red and near infra-red spectral region. The absorption coefficient of 
amorphous selenium has been measured in the ultra-violet visible and near infra-red 
spectral regions. The transmission coefficients of several single crystals of selenium have 
been measured and an estimate made of the absorption coefficient in the wavelength range 
0-68 to 2°04. ‘The results show that the absorption edge of the crystalline material occurs 
at about the’ same wavelength as that for the amorphous form. A discrepancy between the 
value of the near infra-red refractive index and the dielectric constant of amorphous selenium 
indicates a further absorption band in the infra-red. . 
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On a New Test Method for Spherical Aberration of Electron Lenses, by 
O. KLEMPERER. 


ABSTRACT. The focus of rays from a lens with spherical aberration appears as a spot 
surrounded by a discrete halo ring if.a diaphragm with fine circular aperture is placed 
across the beam in front of the focus. The diameter of the halo allows an estimate of the 
magnitude of the aberration involved. "The geometry of rays forming the halo is explained 
here by schematic drawings. The practical application of the halo test is illustrated by some 
examples. In particular, the negative spherical aberration produced by an electronic space 
charge in a saddle field lens is demonstrated. 


The Mechanism of Positive Ion Collection by a Spherical Probe in a Dense Gas, by 
RwsbeEs Bovp: 


ABSTRACT. "The Langmuir probe technique is not suitable for measuring ion densities 
above 1 mm. Hg pressure because the probe dimensions approach those of the ionic and 
electronic mean free paths.. In addition to the invalidity of the Langmuir theory for this 
case, such a probe also causes a great disturbance of the discharge. 

In this paper a detailed examination is made of the possibility of using a probe collecting 
positive ions as a means of finding ion densities. 

It has been found possible to calculate the potential distribution, outside a space-charge 
sheath around a spherical probe, if the radius of the sheath is known. From this the current 
of positive ions to the probe may be found. Curves are given to facilitate this calculation. 

It is found that, depending on the ion concentration and pressure, there are two rather 
different sets of conditions around the probe. If the ion density is high (greater than about 
108/mean free path (in cm.)) a thin sheath will form on the probe and a complete solution of 
the problem is possible. Under such conditions the energies of the diffusing ions in the 
neighbourhood of the probe greatly exceed the thermal energies of the gas particles. At the 
lower ion densities normally encountered, a thick ‘ diffusion’ sheath occurs and the problem 
can now only be solved if the sheath thickness is known. The ion energies may or may not 
greatly exceed the thermal energies, depending on conditions. 

These results contradict one of the basic assumptions made by Davydov and Zmanov- 
skaja in their approach to the same problem. They assume a thin sheath with the ions in 
thermal equilibrium with the surrounding gas. As a result their conclusion is in conflict 
with the conclusion of this paper, that in most circumstances itis not possible to use a negative 
probe to measure ion densities unless an additional means of determining the sheath 
thickness is available. 


The Temperature of the Upper Atmosphere, by D. R. Bares. 


ABSTRACT. ‘The properties of a number of simple models of the upper atmosphere 
are tabulated. A study is made of the thermal equilibrium in the region of the F layers. 
The rate at which energy is gained from ionizing photons is estimated. Various loss 
processes are examined. It is found that the most important is conduction, and the next 
most important is probably emission by the magnetic dipole connecting the two low levels of 
the ground term of atomic oxygen. ‘These are so effective indeed, that the energy source 
first considered appears inadequate to maintain the high temperature the upper atmosphere is 
generally supposed to possess. Alternatives are investigated. Sufficient energy might 
conceivably be supplied either by band absorption by nitric oxide, or by incoming interstellar 
matter : but neither source is attractive. Attention is drawn to the possibility that the original 
estimate of the contribution from ionizing photons may be muchtoo small. This estimate is 
based on radio measurements. It is shown that these only give a lower limit to the electron 
production rate : for certain ions may be removed so rapidly by recombination that they 
would escape detection even if their formation rate were greater than that of the dominant 
ions in the layers. The suggestion is tentatively made that heat is supplied to the upper 
atmosphere mainly by non-observed ionization. It is pointed out that owing to the possible 
existence of such ionization, the measurements of the radio scientists cannot be assumed to 
give the intensity of the solar emission beyond the Lyman limit. The temporal variation of 
the temperature is briefly discussed. 
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